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ABSTRACT 


Condensation of R-113 was studied using an evaporator/condenser test platform. 
The condenser section contained four horizontal tubes (nominal outer diameter 15.9 
mm) forming a vertical in-line column with a pitch-to-diameter ratio of 2.25 and a con- 
densing length of 1.2 m. The condenser tubes could be operated either individually (ie. 
as a single tube apparatus) or as a small tube bundle. This allowed investigation of the 
effects of condensate inundation on different types of condenser tubes. Tubes tested were 
smooth copper tubes, copper/nickel KORODENSE roped tubes, KORODENSE tubes 
wrapped with wire, and copper/nickel finned tubes (26 fpi). The outside heat transfer 
coefficient, h,, was calculated by experimentally determining the overall heat transfer 
coefficient, U, , and then using a modified Wilson Plot procedure. Great care was taken 
to ensure the results were not vitiated by the presence of noncondensibles. 

Results obtained with the smooth copper tubes are in agreement with published data 
and verify satisfactory operation of the test platform. Furthermore, problems associated 
with the apparatus encountered by previous workers were successfully overcome. In 
comparison to the top smooth tubes, the copper/nickel top KORODENSE tube yielded 
about a 22% increase in h,. When different diameters of wire were wrapped around the 
copper/nickel KORODENSE tubes, an optimum pitch-to-wire diameter of 7 was found 
yielding almost a 90% increase in A,, compared to the top smooth tube. This marked 
increase is presumably due to surface tension effects thinning the condensate film. The 
copper/nickel finned tubes gave the best enhancement (approximately a 7 fold increase). 
Comparison of the data obtained from the top copper;nickel finned tubes agreed well 
with the model of Beatty and Katz. 

The effect of condensate inundation was to reduce h, for the lower tubes compared 
to the top tube in the bundle. Comparison of h, for the second tube compared to the top 
tube showed that the effects of condensate inundation are reduced most by wrapping the 
KORODENSE tubes with fine wire. 
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NOMENCLATURE 


Standard Symbols 


ee iat sem Se ey (2) Scares 


SN er > 


Tee 


area, m?’ 

fin spacing, m 

parameter defined in Eq (2-11) 

specific heat, J/kg K 

diameter, m 

fin height, m 

gravitational acceleration, m/s? 

heat transfer coefficient, W/m? K 

heat transfer coefficient in flooded portion of finned tube, W/m? K 
heat transfer coeffcient in unflooded region of finned tube, W/m? K 
heat transfer coefficient for fins in unflooded region, W/m? K 
thermal conductivity, W/m K 

length, m 

average condensing length defined by Eq (2-4) 

mass flow rate, kg/s 

number of tubes in a vertical column 

Nusselt number 

Prandtl number 

heat transfer rate, W 

heat flux, W/m? 

radius, m 

fouling resistance, K/'W 

wall resistance, K/W 

Reynolds number 

film Reynolds number 

temperature, K 

fin thickness, m 


overall heat transfer coefficient, W/m? K 


V 


Greek Symbols 


Subscripts 
€ 
eq 
ef 
f 
fs 
ft 


Sat 


velocity, m/s 


tape thickness, m 

inreriieiency 

aspect ratio of condensate surface, Eq (2-13) 
rotation angle normal to fin surface 
dynamic viscosity, kg/m s 

density, kg/m 

surface tension, N/m 


flooding angle defined in Eq (2-7) 


coolant 
equivalent 
effective 

film 

fin sides 

fin tips 

inside or inlet 
outside or outlet 
root 


saturation 





I. INTRODUCTION 


A. BACKGROUND 

In the last few vears, budgetary considerations have led to a reduction in the funding 
level available for new vessel construction for the US Navy. In construction of ships, 
weight can often be equated to cost. This has led the Navy to look at various ways by 
which the weight of auxiliary systems can be reduced. 

Of particular interest to the Navv has been the issue of reduced weight and therefore 
cost of refrigeration systems. Two major components of refrigeration systems are the 
evaporator and condenser. If the weight of these two components can be reduced, then 
the overall weight and presumably cost, can also be reduced. 

Condenser weight can be significantly reduced by the proper selection of materials 
for condenser construction. A major weight savings comes from changing the material 
which is used for construction of the condenser tubes. As an example, the DDG-51 was 
originally intended to have a refrigeration condenser made of titanium finned tubes. This 
would provide a substantial weight savings and a reduced cost (presumably). However, 
the system) was not installed on the current flight of ships. Rather, traditional 
copper/nickel finned tubes were utilized for the refrigeration condenser.! 

In re-designing refrigeration condensers, changes in condenser tube material based 
solely on weight consideration should not be the only consideration. The reason for this 
is that a lightweight material may not have heat transfer characteristics which are com- 
parable to copper or even copper/nickel. Therefore, condenser tube designs have to be 
developed to maintain or even enhance heat transfer or else the condenser will have to 
be enlarged. This enlargement is necessary in order to provide the requisite surface area 
to maintain the required level of heat transfer. 

Copper/nickel finned condenser tubes have been used in refrigeration condensers 
onboard US Navy ships for many years. However, the use of titanium, as was planned 
for the DDG-51 class ships, holds the promise of providing substantial weight savings 
(not to mention improved corrosion resistance on the tube side). Unfortunately, 
titanium has a thermal conductivity which is substantially less than that of copper or 


copper/nickel. The question which must therefore be addressed is whether titanium 


1 According to Mr. John Goodhue, NAVSEA code SEA 56Y15, the reason the units were 
not installed was that they could not be ready to support proposed construction dates for flight one. 


condenser tubes can be constructed in such a manner as to maintain heat transfer (due 
to enhancement) despite the lower value of thermal conductivity. 

A major problem along these lines is how to measure the performance of various 
types of condenser tubes. Typically, this is done with single tubes or small tube bundles. 
These tubes tend to be relatively short. Often, the experimental apparatus bares little 
resemblance to a shipboard condenser. Thus, any results become suspect when they are 
extended to the real situation. 

In order to examine heat transfer characteristics of various designs of condenser 
tubes, the best means Is to construct a test platform composed of an evaporator and 
condenser with a size comparable to that of a shipboard svstem. Unfortunately, this is 
seldom feasible due to costs involved and the size of the plants. However, it 1s possible 
to construct a reasonably large scale test platform in which tube performance can be 
evaluated. This evaluation must include both single tube performance as well as tube 
bundle operation. The reason for this 1s that condensate forming on one tube falls onto 
lower tubes. This effect, known as condensate inundation, can result in decreased heat 
transfer on lower tubes due to a thickening of the condensate layer. Hence, another 
aspect Which must be considered 1s whether condenser tube geometry can be changed in 
such a manner as to reduce the adverse effects of condensate inundation. This could 
help make up for decreased thermal conductivity associated with lightweight materials 
such as titanium. | 

The construction of a suitable condenser/evaporator test platform was begun at the 
Naval Postgraduate School by Zebrowsk: [Ref. 1]. The construction and initial operation 
of the test platform was completed by Mabrey [Ref. 2] in 1988. Mabrey conducted ex- 
periments using R-113 and R-114 on smooth copper tubes as well as copper/nickel 
finned tubes (26 fpi). During the course of his experiments, an orange-peel or dimpled 
appearance was observed on the smooth copper tubes. Additionally, the copper/nickel 
finned tubes showed a buildup of “material” between the fins. This material was believed 
to be ice although this was not proved. The cause of the orange-peel appearance was 


also never identified. 


B. OBJECTIVES. 

This thesis concentrates on the condenser section of the test platform with specific 
interest in the condensation of R-113. R-113 was chosen for this study because of 
convienence , ie., it is a liquid at atmospheric pressure. The condenser section consists 


of four horizontal tubes mounted in a vertical in-line column. The nominal] diameter of 


the tubes is 16 mm with a condensing length of 1.2 m. While the condensing medium 
used in this study was R-113, the condensation of many different types of fluids can be 
evaluated. In addition, the condenser can be used to study heat transfer on almost any 
type of condenser tube subject to the nominal outer diameter and condensing length re- 
strictions indicated. 

The main objective of this thesis was to resolve the contamination problem found 
by Mabrey [Ref. 2], and to make the test platform operational. It was felt that such a 
system, once operational, could be used to study heat transfer using condenser tubes 
made of various materials and geometries. In addition, since the tubes are arranged in 


a vertical in-line columin, the effect of condensate inundation could also be studied. 


C. SPECIFIC OBJECTIVES 
The specific objectives of this study were as follows: 
1. Resolve the contamination problem reported by Mabrey [Ref. 2] 


2. Establish bascline condensation data by examining the condensation of R-113 on 
_ smooth copper tubes. 


3. Examine the heat transfer characteristics of R-113 on various tube geometries 


4. Investigate the effects of condensate inundation 


I]. LITERATURE SURVEY 


A. INTRODUCTION 

Condensation occurs when vapor molecules come into contact with a surface which 
has a temperature below the vapor saturation temperature. This vapor, coming into 
contact with the cooler surface, forms a liquid which can accumulate in two ways. If the 
liquid wets the entire surface, a continuous film of liquid is formed. On the other hand, 
if the liquid fails to wet the entire surface, droplets are formed. These two types of 
condensation are referred to as filmwise and dropwise condensation respectively. 
Dropwise condensation yields larger heat transfer coefficients but is difficult to maintain 
(Ref. 3]. All condensers today are therefore designed to operate with film condensation 
conditions. 

In 1916, Nusselt [Ref: 4] reasoned that as vapor condenses on a cooler surface, such 
as a horizontal tube, a continuous film of condensate will form. This film grows in 
thickness as the condensate flows around the tube under the influence of gravity. Hence, 
the thickness of the condensate 1s not uniform around a horizontal tube. Since the 
condensate acts as a resistance to heat transfer, Nusselt argued that the local heat 
transfer coefficient will be a maximum at the top of the tube where the film thickness is 


least and decreases around the tube as the condensate thickness increases. 


B. CONDENSATION OF REFRIGERANTS ON SMOOTH TUBES 

Goto et al. (Ref. 5] examined film condensation of various refrigerants, including 
R-113, on smooth copper tubes. Their results showed values for the vapor-side heat 
transfer coefficient on the order of 1.2 to 1.5 KW/m? K for a wall-film temperature dif- 
ference of 3 to 25 K. Furthermore, their results were about 10% below that predicted 
by Nusselt theory. 

. Masuda and Rose [Ref. 6,7] reported that condensation of R-113 on smooth copper 
tubes yielded an outside heat transfer coefficient of about 1500 W/m? K for wall-film 
temperature differences of 12 to 20 K. Their results showed good agreement with 
Nusselt theory. 

Marto et al. [Ref. 8] showed that for condensation of R-113 on smooth copper 
tubes, the vapor-side heat transfer coefficient was on the order of 1000 W/m? K for a 
temperature drop across the film of 20 to 25 K. Their values were also in good agreement 
with Nusselt theory. 


C. ENHANCEMENT OF FILM CONDENSATION 

The most commonly used procedure for enhancement of film condensation is 
through the use of integral-fin tubes. With this type of condenser tube, two distinct re- 
gions exist; namely, a flooded and an unflooded region. Typically, as fluid condenses 
on a finned surface, it is pulled from the fin tip to the fin root by surface tension forces. 
Gravity then causes the fluid between fins to drain to the bottom of the tube where it 
drains off. Hence, condensation on a finned tube involves a combination of surface 
tension and gravity forces. It is the added surface tension which tends to enhance heat 
transfer over smooth tubes by thinning the condensate film at the fin tip. 

Since the present study uses R-113 as the working fluid, the remainder of this 
chapter will focus on previous experiments which have used low surface tension fluids 
such as refrigerants. For additional information involving other fluids, the reader is re- 
ferred to recent review articles by Marto [Ref. 9,10], Webb [Ref. 11] and Sukhatme [Ref. 
i). 

1. Condensation on Single Enhanced Tubes: Theoretical Models 

The first theoretical model of condensation on finned tubes was published by 


Beatty and Katz [Ref. 13}. Their theoretical model is given by: 


hex = py e+ hy cs (2 — 1) 
Bey Agr 





Here 7, 1s the fin efficiency, A, is the fin surface area, A, is the unfinned surface area, and 


A,, 1s the effective total surface area. The above equation can be simplified to: 
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where the equivalent diameter is given by: 
1 \« Ar 1, Ao 1 
= 1307 + 2—3 
( Deg Ag DHA” Ay pi4 waa 
In this equation, 
_ 2(D;—D;) 
L= Top) ae (2 — 4) 


and 
Ag = nA + wy (2 — 5) 


The equation for the heat transfer coefficient 1s a Nusselt-type condensation equation 
based on the equivalent diameter for a finned tube. This equation predicted the exper- 
imental results obtained by Beatty and Katz to within + 20%. 

There are several important limitations to the Beatty and Katz model [Ref. 13]. 
It does not distinguish between fin tips and fin flanks, essentially assuming that the fins 
have no thickness. Therefore the model 1s not suitable for trapezoidal fins. A second, and 
perhaps more significant limitation is that the Beatty and Katz model fails to take into 
account surface tension effects. Instead, the model assumes that the flow of the 
condensate film is simply a function of gravity. While this may not pose a Serious limi- 
tation when low surface tension fluids such as refrigerants are examined, this assumption 
can become critical when condensation of high surface tension fluids such as steam are 
investigated. 

Another limitation of the Beatty and Katz model is that it does not take into 
account condensate flooding. Rudy and Webb [Ref. 14] argued that the flooded portion 
of the tube was ineffective in condensation. They therefore proposed that the Beatty and 
Katz model be modified by: 


Py 


haw = higx “Tt (Zao) 


where @, the flooding angle, is defined as: 


- 20 
bs = cos ‘lier et (2 —7) 


where A, 1s the outside radius, w is the fin tip half angle and a is the fluid surface tension. 
It is interesting to note that this modification resulted in an underprediction of Rudy and 
Webb's R-11 condensation data by about 35%. 

Webb et al. [Ref. 15] next argued that surface tension drainage must be consid- 
ered in the model. They derived the equation: 


1 


k3 oh + aa 
fP"e fa) 26 f 1, tl \i4 
y= 0983} LAT Da caae) Ve) 





where ¢ is the fin height, t is the fin thickness and b 1s the width of the space between 
the fins. Substitution into the Beatty and Katz expression yielded agreement with ex- 
perimental data of + 10%. 

Owen et al. [Ref. 16] argued that the Rudy and Webb [Ref. 15] model had to 
be corrected in order to account for heat transfer in the flooded part of the tube. Owen 


et al. proposed the equation: 


b= hax e+ h(1-) (2 — 9) 


where h, is the average outside heat transfer coefficient and h, is the heat transfer coef- 
ficient in the flooded region of the tube. This value was estimated using a one- 
dimensional radial conduction model. This corrected model of Owen et al. resulted in 
agreement of + 30% with earlier experimental data. 

In 1985, Webb et al. [Ref. 17] proposed a model that included heat transfer in 
the flooded region as well as surface tension drainage on the fin flanks and gravity 
drainage on the interfin surface. Their model equation for the average heat transfer co- 


efficient is: 
A, Ar 


In this expression, h, is the heat transfer coefficient for the unflooded root surface be- 
tween the fins and f, is the heat transfer coefficient on the fin flanks in the unflooded 
area, and A, 1s the heat transfer coefficient in the condensate flooded region. The term 
C, is the fraction of finned tube surface flooded by condensate, or, 


bet-—" (2-11) 


where ¢, is the flooding angle given in equation (2-7). 
In the Webb et al. [Ref. 17] expression for the average heat transfer coefficient, 
h,, h, and h, are defined mathematically as follows: 


2 J 
u yil. 
kp Pf 


where Re, is the condensate film Reynolds number. 





hy = 2.149 


ky | ean Et ig Ca 
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where S_ is the length of the convex surface and ©,, is the rotation angle normal to the 
fin surface [Ref. 18]. The term € 1s a parameter that characterizes the aspect ratio of the 
condensate surface. 


Finally, A, 1s defined as: 


This complex model of Webb et al. was able to predict heat transfer coefficients to 
within + 20% of experimental data obtained with R-11. 

Honda and Nozu [Ref: 19] and Honda et al. [Ref. 20] have proposed the most 
comprehensive model to date. This model treats the heat transfer problem as one in- 
volving vapor-to-coolant heat transfer through a finned tube wall. In this model, 
trapezoidal fins were divided into flooded (f) and unflooded (u) regions. Equation (2-7) 
was used to determine the flooding angle. Condensate flow in between the fins was as- 
sumed to be driven by gravity whereas for the fin surface, condensate flow was assumed 
to be driven by both surface tension forces and gravity. The condensate film thickness 
on the fin surface was arbitrarily divided into a thick region, where no heat transfer oc- 
curred, and a thin film region. The equation of the condensate film thickness was solved 
numerically. 

Honda and Nozu [Ref. 19] showed that the average Nusselt number could be 


Written in terms of relevant parameters for the flooded (f) and unflooded (u) regions as: 
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Here, 7,, and T7,, are the dimensionless average wall temperatures at the fin root and @¢, 


is defined as: 


~ 
b= (2 — 16) 


It should be pointed out that T., and a are determined by solving a circumferential wall 
heat conduction equation. This wall temperature variation must be considered since heat 
transfer rates through flooded and unflooded regions are very different. 

Honda and Nozu [Ref. 19] used their model to predict average heat transfer 
coefficients. These predictions were compared to experimental data on eleven different 
fluids and twenty-two tubes. The predicted values matched most experimental values to 
within + 20%. 

2. Condensation on Single Enhanced Tubes: Experimental Results 

The first experiments conducted on enhanced tubes were performed by Katz and 
colleagues [Ref. 21,22, 13 ]. These authors looked at condensation of R-12 and R-22 on 
smooth and finned tubes. Their results showed a significant enhancement in heat transfer 
for finned tubes. 

In 1971, Karkhu and Borovkov [Ref. 23] examined condensation of R-113 on 
five finned tubes. Their results showed that the rate of condensation for these tubes with 
trapezoidal fins could be increased by 50 to 100% over smooth tubes. 

Arai et al. (Ref. 24] and Kisaragi et al. [Ref. 25] studied the condensation of R-12 
and R-113 respectively on finned tubes with spiral grooves cut along the length of the 
tube. Arai et al. showed that grooves with a pitch of 0.7 mm gave the best enhancement 
in heat transfer. Kisaragi and co-workers also compared grooved finned tubes with tubes 
on which the fins had been covered with a porous coating. In comparison to the finned 
tube, a two fold enhancement was seen with the grooved finned tube while a four fold 
enhancement was obtained with the finned porous coated tube. 

~ Carnavos [Ref. 26] performed one of the first detailed studies of condensation 
of R-11 on enhanced copper tubes. In this study, eleven different tubes were tested. The 
tubes included finned tubes as well as finned tubes which had an axial spiral groove cut 
along their length. Various sizes and shapes of fins were included. Enhancements of 4 
to 6 were seen when compared to smooth tubes. 

Rudy [Ref. 27] and Webb et al. [Ref: 15] reported experimental data obtained 
for R-11 condensing on copper tubes with different fin geometries.. Their results showed 
that increasing the fin density from 748 to 1378 fpm resulted in about a 70% increase in 
heat transfer presumably due to an increased surface area. 

Honda et al. [Ref. 28] studied condensation of R-113 on copper finned tubes of 
various fin density. Their results showed a 6 to 9 fold enhancement for R-113 
condensation for finned tubes when compared to smooth tubes. A similar study by 
Kabov [Ref. 29] examined the condensation of R-12 and R-22 on tubes with different fin 


densities His results also showed that the degree of enhancement is determined by fin 
spacing. 

Masuda and Rose [Ref. 6] examined heat transfer for the condensation of R-113 
on low integral-fin copper tubes. Fourteen different fin spacings (0.25 to 20 mm) were 
examined. The optimal fin spacing was 0.5 mm which resulted in a 7.3 fold enhancement 
in the vapor-side heat transfer coefficient. 

Lin and Berghmans [Ref. 30] examined condensation of R-11, R-113, and R-114 
on an enhanced tube. This enhanced tube, known commercially as “EVERFIN’, 1s made 
of trapezoidal fins (1378 fpm) with a spacing between the fins of 0.064 mm.2 This tube 
gave an 8 to 10 fold increase in the vapor-side heat transfer coefficient when compared 
to a smooth tube. 

Condensation of R-11 on a variety of enhanced copper tubes has been reported 
recently by Sukhatme et al. [Ref. 31]. All tubes had trapezoidal shaped integral fins. The 
fin density, fin height and fin-tip half angle were systematically varied in order to deter- 
mine the optimal values which would result in maximum enhancement. The maximum 
enhancement obtained by these authors was about 10.3 relative to smooth tubes. This 
enhancement corresponded to a fin spacing of 0.35 mm, a fin height of 1.22 mm and a 
fin-tip half angle of 10 degrees. In addition, a set of specially enhanced tubes were tested 
by these authors. These tubes had axial grooves cut into them with the trapezoidal fins 
present. An enhancement of 12.3 was obtained with these specially enhanced tubes 
when compared to smooth tubes. 

Marto et al. [Ref. 8] studied condensation of R-113 on a variety of copper 
integral-fin tubes with rectangularly shaped fins. The optimum fin spacing was found 
to lie between 0.25 and 0.5 mm. The vapor-side heat transfer coefficient was enhanced 
on the order of 4 to 7 fold. 

| Michael et al. (Ref. 32] continued the studies of Marto et al. (Ref. 8] by looking 
at the effect of root diameter on the condensation of R-113 on integral-fin tubes. They 
found that the root diameter had only a negligible effect on heat transfer enhancement. 
In addition, these investigators found that the optimal fin spacing was about 0.5 mm. 


They obtained an enhancement of 5 to 6. 


_ 2 While the authors quote a fin spacing of 0.064 m, this value seems too small. The schematic 
diagram of the fins together with other dimensions published by these authors tend to suggest that 
the correct fin spacing is 0.64 mm. This valuc is more in line with values quoted by other authors. 
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In comparing enhancements reported by various authors, a word of caution is 
required. In calculating values for the vapor-side heat transfer coefficients, the surface 
area to be used in this calculation must be defined. This area can be based on outside 
or envelope diameter, root diameter or total fin area. The great majority of authors 
define this area by the outside diameter of the tube as opposed to the root diameter. In 
the various sources cited above, this has been the case. That is, the surface area has been 
defined by the outside tube diameter. 

While it is clear that integral-fin tubes can result in an enhancement of the 
vapor-side heat transfer coefficient, other alternatives exist. As an example, Fujii et al. 
[Ref. 33] studied the condensation of R-11 on smooth tubes which were wrapped with 
fine wire. These authors found an optimal pitch-to-wire diameter ratio of 2 which re- 
sulted in an enhancement of 2.5 in the vapor-side heat transfer coefficient when com- 
pared to copper smooth tubes. These authors argued that the enhancement in the heat 
transfer coefficient was due to a thinning of the condensate film. The reason that the 
condensate film thickness is reduced is that the presence of the wire creates a low pres- 
sure region in the film near the base of the wires. This low pressure region is due to the 
concave shape of the condensate film next to the wire. As a result, condensate film is 
pulled along the tube wall towards the wire. This results in a thinning of the condensate 
film between the wires. 

3. Summary of Single Enhanced Tube Condensation 

It 1s clear from the preceding discussion that theoretical models, despite their 
degree of sophistication, still only match experimental data to within + 20%. This 
suggests that the process of condensation on enhanced tubes is extremely complex and 
far from being fully understood. 

In terms of experimental data involving condensation of refrigerants on en- 
hanced tubes, several comments need to be made. The experiments reported here have 
shown that integral-fin tubes can give enhancements of the vapor-side heat transfer co- 
efficient on the order of 7 to 12. The magnitude of this enhancement depends upon such 
factors as fin shape, fin spacing, and fin height. Trapezoidal shaped fins with a spacing 
of 0.5 mm give the highest enhancements. 

Unfortunately, the above experiments were all performed on enhanced tubes 
made of copper. In Naval refrigeration svstems, where seawater is used as the coolant, 
corrosion is of major concern. This necessitates the need for copper/nickel or titanium 
coldenser tubes in. order to minimize corrosion. However, as can be seen from the 


foregoing summary of the literature, essentially no work has been done on the 
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condensation of refrigerants on copper/nickel or titantum finned tubes. Since both 
copper/nickel and titanium condenser tubes have a lower thermal conductivity than 
copper, one must wonder what effect this will have on the vapor-side heat transfer co- 
efficient. Clearly, experiments using these tubes must be carried out. 

4. Condensation on Enhanced Tube Bundles: Theoretical Studies 

Condensation in tube bundles presents an entirely different situation than 
condensation on single tubes. In a tube bundle, two conflicting factors come into play. 
First, condensate forming on an upper tube drains off this tube onto lower tubes within 
the bundle. This increases the thickness of the condensate layer on lower tubes which 
further increases the resistance of these tubes to heat transfer. This effect is referred to 
as condensate inundation. 

Condensate falling from one tube to another can have a second effect. As 
droplets fall from one tube to the next, ripples and splashing (turbulence) are induced in 
the condensate film. This can actually enhance heat transfer on the lower tubes. 

Vapor velocity can also be a significant factor in tube bundle condensation. As 
the vapor travels into the tube bundle and condenses, the local velocity decreases. This 
decrease in loca! velocity leads to a decrease in the heat transfer coefficient similar to an 
inundation effect since vapor shear effects that thin the condensate film are reduced. 

In analyzing condensation in horizontal tube bundles, Nusselt [Ref. 4] assumed 
that all condensate flowed from one tube to another as a continuous laminar sheet. He 
showed that the ratio of the average heat transfer coefficient for a vertical column of 


horizontal tubes with respect to the top tube could be expressed as: 
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where N is the number of tubes in the vertical column. A local coefficient can also be 


calculated for the Nth tube using Nusselt theory. This coefficient is given by: 

—~=N* —(N-1)* @= 15 
Kern [Ref. 34] proposed a less conservative approach. He argued that droplets 

of condensate falling from one tube to another causes ripples which disturb the 


condensate film, thereby diminishing the effect of condensate inundation. His math- 


ematical model took the form: 
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The local value of h for the Nth tube can be expressed as: 
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In 1972, Eissenberg [Ref. 35] proposed that condensate may not drain only in 
a vertical direction, but rather, mav be diverted sideways, presumably due to local vapor 
flow conditions. This so-called “side drainage” model predicts a less severe effect of 


inundation yields: 
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At the present time, numerous investigations have been carried out in an at- 
tempt to better understand the phenomenon of condensate inundation. However, the 
data are scattered and as a consequence, no successful theoretical model of condensate 
inundation exists [Ref. 3]. However, for design purposes, Butterworth [Ref. 36] has re- 
commended that the Kern model be used. 

5. Condensation in Enhanced Tube Bundles: Experimental Studies 

Young and Wohlenberg [Ref. 37] examined the condensation of R-12 on a 
vertical column of smooth horizontal tubes. These investigators found that the effect of 
the number of tubes in the vertical column is less than that predicted by Nusselt. 

Katz et al. [Ref. 38] investigated the condensation of R-12 on a vertical column 
of six horizontal finned tubes. Their results also showed a less severe inundation effect 
than predicted by the model of Nusselt. 

Smirnov and Lukanov [Ref. 39] used R-11 to study condensation on 20 hori- 
zontal rows of finned tubes arranged in a triangular fashion. Their results showed a 
dramatic reduction in heat transfer for the first five tubes. However, after this, the heat 
transfer coefficient remained almost constant for the rest of the tubes in the column. 

Honda et al. [Ref. 40] have studied the condensation of R-113 on a vertical 
column of horizontal low finned tubes. An inundation tube was used to simulate addi- 
tional rows of tubes. These authors described four modes of condensate flow; namely, 


Groplet mode, column mode, column and sheet mode, and sheet mode. Unfortunately, 
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these authors did not calculate the heat transfer coefficient so that no comparison with 
theoretical models is possible. 

Webb and Murawski [Ref. 41] examined the condensation of R-11 on a column 
of enhanced condenser tubes. The tubes tested were a standard 1024 fpi tube, a 
GEWA-SC tube, a Turbo-C tube and a TRED-D tube. The GEWA-SC tube is a 1024 
fpi tube but the fins have a Y-shaped cross-section. The Turbo_c and TRED-D tubes 
are similar to the 1024 fpi tube except that the fins have a saw-tooth shape. Webb and 
Murawski found that a standard 1024 fpm tube showed virtually no row (inundation) 
effect. However, when the experiments were conducted on Turbo-C, GEWA-SC and 
Tred-D tubes, row effects (ie., inundation effects) were seen. The inundation effect was 
greatest for the Tred-D tube which also had the lowest single tube performance. These 
investigators found values for the exponent in equations (2-17, 2-19, and 2-21) on the 
order of 0.12 to 0.26. Recall that Nusselt’s theory used a value of 0.25. 

More recently, Honda et al. [Ref: 42] studied the condensation of R-113 on an 
in-line horizontal tube bundle. The tubes tested consisted of either flat-sided (rectangu- 
lar) annular fins or three dimensional trapezoidal fins. These authors found that the ef- 
fects of condensate inundation were greatest for the three dimensional fin tubes. 
Unfortunately, these authors also had vapor velocities on the order of 3.4 m/s or greater. 
Therefore, it 1s difficult to compare their results directly to the theoretical models of 
Nusselt, Kern, or Eissenberg. | . 

An interesting facet of the experiments of Honda et al. [Ref. 42] deserves men- 
tion. These investigators compared the effects of condensate inundation on the vapor- 
side heat transfer coefficient for smooth and finned copper tubes. Their results showed 
that despite condensate inundation, the heat transfer coefficient for the finned tubes was 
still substantially above that of the unfinned tubes. In addition, the decrease in the heat 
transfer coefficient with increasing row number was less with the finned tubes suggesting 
that the effects of condensate inundation were less severe with finned tubes than with 
unfinned tubes. Presumably, this result is due to the fact that condensate draining off 
of a finned tube drains in columns rather than continuous sheets. Lower tubes will 
therefore not be completely covered with condensate but instead will have unflooded 
regions. 

Honda and Nozu [Ref. 43] examined the effects of condensate inundation on 
smooth and integral-fin copper tubes using a theoretical model derived recently by 
Honda et al. [Ref. 44] This model was a single tube model which was extended to include 


the effects of condensate inundation by realizing that condensate did not flow off tubes 


as sheets but rather as columns. The model took account of the fact that regions which 
received drainage of condensate had different heat transfer coefficients than regions 
which did not. The results support the more recent experimental results of Hlonda et al. 
[Ref. 42]. The effect of inundation was less severe for the finned tubes than for the 
smooth tubes. The authors suggest that this 1s due to condensate being channeled on the 
finned tubes whereas for the smooth tubes, the entire tube surface becomes covered with 
a condensate film of varying thickness. 
6. Summary of Enhanced Tube Bundle Condensation 

Experimental data has shown that the Nusselt model for the effects of 
condensate inundation may be too conservative. Rather, most experimental data seems 
to fall between the Kern and Eissenberg predictions. 

Other experimental data suggests that the effects of condensate inundation are 
less for finned tubes as opposed to smooth tubes. Presumably, this 1s due to the manner 
in which condensate flows from one finned tube to another. The presence of the seems 
acts to channel the flow of condensate leaving large areas relatively free of flooding from 
inundation. 

As a final comment, the use of wire wrapped tubes provides an interesting al- 
ternative to the more conservative finned tube approach. However, the effects of 


refrigerant condensate inundation on wire wrapped tubes has yet to be examined. 
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I. EXPERIMENTAL APPARATUS 


A. TEST PLATFORM 

Details of the design and construction of the condenser/evaporator test platform are 
provided by Zebrowski [Ref. 1] and Mabrey [Ref. 2]. Therefore, only a brief summary 
will be presented here. 

The test platform is shown scheinatically in Figure 1 while a photograph of the unit 
is shown in Figure 2. The unit is made from rolled stainless steel plate (6.35 mm thick) 
and designed to withstand an absolute pressure of 308 kPa (approximately 3 bar). 

1. Condenser Section 

The top section is the condenser. It is 1.3 min length with an internal diameter 
of 0.61 m. The end plates are made of stainless steel plate (6.35 mm thickness) with a 
diameter of 0.71 m. They are connected to the condenser shell using a bolted flange as- 
sembly. A rubber gasket, 3.22 mm thick, provides a vacuum tight seal. The condenser 
unit has five viewports which allow the condensation process to be viewed. 

Four test condenser tubes are contained in the condenser forming a Vertical in- 
line column with a pitch-to-diameter ratio of 2.25. These are held in place at each end 
by a nylon block bolted to the end plate. A rubber gasket provides a vacuum tight seal. 
A Stainless steel plate is fitted over the ends of the condenser tubes and bolted to the 
nylon block. O-rings are placed in between the nylon block and steel plate in order to 
provide a vacuum tight seal. 

In order to provide a means of pressure control during testing of the condenser 
tubes, the test platform was fitted with five auxiliary condenser tubes. These auxiliary 
condenser tubes are constructed of 9.53 mm diameter copper tube, wound into coils. 
Each coil is supported by a stainless steel rod welded to one of the end plates. The aux- 
iliary condenser coils enter and exit through this same condenser end plate. 

A stainless steel shroud assembly, shown in Figure 3, is installed in the 
condenser unit. This shroud surrounds the test condenser tubes as well as the auxiliary 
condenser coils. A window, installed in the side of the shroud, allows visualization of the 
test condenser tubes through the viewports. The purpose of this shroud is to funnel va- 
por from the evaporator below, around the outside of the shroud and then vertically 
down past the condenser tubes. The velocity of the vapor past the condenser tubes was 
0.1 m/s. 
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Figure 1. Schematic diagram of test platform. 
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Figure 3. Schematic diagram of condenser showing shroud. 
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2. Evaporator Section 

The flooded evaporator is formed from rolled stainless steel and has an internal 
diameter of 0.61 m with a width of 0.279 m. Two viewports allow visualization of the 
boiling process. This unit is connected to the condenser by a collar assembly 0.28 m in 
diameter and 0.203 m in length. The volume of liquid contained in the evaporator is 
approximately 0.06 m3, just sufficient to cover the boiler tubes. 

Vapor is produced in the evaporator using a combination of three groups of 
heaters. The first group is the main bundle and is composed of thirty-five tubes in a 
triangular pitch (pitch to diameter ratio of 1.2). Fifteen of these have heaters with a 
nominal rating of 1 kW each. Five of these heated tubes are instrumented with 
thermocouples to enable bundle boiling data to be obtained. In addition, unheated 
dummy tubes are used to guide the two phase flow within the bundle. 

The second group of tubes 1s composed of five simulation heaters, each with a 
nominal rating of 4 kW. These are placed below the main bundle in order to simulate 
further rows of tubes. 

‘The final group consists of four auxiliary heaters, each with a nominal rating 
of 4 kW. Two of these are placed on each side of the evaporator in order to keep the 
liquid pool at saturation conditions; this avoids any subcooling of the liquid. Power to 
all three sets of heaters is provided by three separate STACO 240V, 23.5kVA rheostat 
contollers. . | 

For further information on the design and construction of the test platform, the 
rcader 1s referred to Zebrowski [Ref. 1], Mabrey [Ref: 2], Murphy [Ref. 45], and Chilman 
[Ref. 46]. 

3. Coolant System 

The coolant is a 54/46 by weight mixture of ethylene glycol and water. The 
coolant system consists of two main coolant pumps. These pumps are 0.5 hp constant 
speed pumps which take a suction from the inain coolant sump. One pump is used to 
supply five rotameters. A ball valve, located at the entrance to each rotameter, is used 
to control the flow rate. Flow to each test condenser tube is supplied by an individual 
rotameter. The fifth rotameter is used to supply one of the auxiliary condenser coils. A 
needle valve is used to provide fine control of flow through the coil. This coil is normally 
used to control saturation pressure during experimental runs. 

The second coolant pump is used to supply coolant flow to the remaining four 
auxiliary condenser coils. Flow through each coil is controlled by a needle valve located 


at the entrance to each coil. It should be noted that the four auxiliary cooling coils are 


seldom needed during the condenser experiments; their primary purpose is for use during 
evaporator experiments. (see Chilman [Ref. 46)). 

In addition to supplving the four auxiliary condenser coils, the second coolant 
pump also provides coolant flow to a refrigerant storage tank and to a cold trap located 
at the top of the condenser. The storage tank can be used to store refrigerant during 
condenser or evaporator tube changeout. The cold trap is placed in-line between the 
condenser and a vacuum pump. Vapor from the condenser is drawn through the cold 
trap which is cooled by the ethylene glycol/water mixture. This process is used to remove 
noncondensibles from the condenser while preventing loss of refrigerant to the atmos- 
paere. 

Flow of the coolant mixture from the auxiliary condenser cooling coils, storage 
tank, and cold trap are returned to the sump tank. Coolant flow from the test condenser 
tubes empties into a common header which then empties into the same coolant sump 
tank. 

In order to keep the coolant mixture temperature approximately constant, it 1s 
cooled by an eight ton refrigeration system located outside the building. This system 
continually recirculates the coolant using a 0.75 hp pump. Coolant temperature can be 
maintained at any temperature from ambient temperature (approximately 20 °C) down 
to -20 °C. 

4. Instrumentation 

Prior to entering each condenser tube, the coolant flows through a mixing 
chamber. The inlet temperature is measured using a single copper/constantan type-T 
thermocouple. Coolant flows from each condenser tube through an outlet mixing 
chamber. Outlet coolant temperatures are then measured using two thermocouples for 
each tube. An average outlet temperature is then taken for each tube. All thermocouples 
had been previously calibrated by Mabrey [Ref. 2]. 

| Refrigerant temperatures Were measured at several different locations in the 
evaporator/condenser test platform using similar single thermocouples. Vapor temper- 
ature was measured in two separate locations at the top of the condenser and one lo- 
cation just above the liquid level. Refrigerant liquid temperature was measured at both 
the top and bottom of the tube bundle. The thermocouples were inserted into stainless 
steel housings (with copper tips) which were installed in the condenser/evaporator test 


platform. 
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B. CLEANING AND DECONTAMINATION PROCEDURES 

As mentioned earlier, there was a contamination problem with the test platform 
during the work of Mabrey [Ref. 2]. Therefore, the experimental apparatus was initially 
disassembled and cleaned in accordance with the following procedure. 

Prior to removing the evaporator, all heaters and thermocouples were disconnected. 
The lagging around the mating flange was then removed. A Porta-lift was connected to 
the evaporator mating flange in order to support the weight of the evaporator prior to 
removing the bolts holding the evaporator and condenser together. After supporting the 
evaporator, the bolts were removed and the evaporator lowered from the condenser unit. 

Referring to Figure 2, the condenser unit 1s capped at each end by a steel end plate. 
The !eft end plate supports the auxiliary condenser units and is therefore removed after 
the right end plate. The test condenser tubes must be removed prior to the disassembly 
of the condenser. The assembly which holds the condenser tubes in place 1s composed 
of a nylon plate bolted to the condenser end plate. A gasket insures a tight seal between 
the nylon and condenser end plates. A stainless steel plate is bolted against the nylon 
plate with O-rings in between which ensure an air tight seal between the environment 
and the inside of the condenser. 

The plates holding the test condenser tubes in place must be removed prior to re- 
moval of the condenser end plates. The stainless steel plates at both ends are removed 
by unbolting the retaining nuts. The nuts holding the nylon plate in place are then re- 
moved. A rubber mallet is used to tap out the test condenser tubes from the left end of 
the platform. This causes the right end nylon plate to move away from the right 
condenser end plate. Screwdrivers can then be used to pry the nylon end plate away from 
the condenser. 

The condenser end plates can be removed after the test condenser tubes have been 
removed. As a word of caution, these end plates weigh in excess of 100 pounds. 
Therefore, they must be supported during the removal process. In these studies, a 
Porta-lift was used for this purpose. The right condenser end plate was removed first, 
followed by the left end plate. Prior to removal of the left end plate, the shroud assembly 
had to be removed. This was facilitated by first removing the thermocouple probes which 
penctrate the top of the condenser unit. The shroud was then pulled out of the 
condenser. 

If the auxiliary condenser tubes are to be removed, then the the left condenser end 
plate will have to be removed. This is because the five auxiliary condenser tubes are at- 


tached to the left end plate. The location of the test platform in relation to the coolant 


to 
to 


sump is such that this end can not be removed with the evaporator/condenser test plat- 
form in its present location. Rather, if the left head must be removed, the entire test 
platform must be rotated 180 degrees about a vertical axis. 

Initial inspection of the condenser unit showed only one small area of discoloration 
which surrounded the middle viewport. No other obvious signs of contamination Were 
evident. The condenser and evaporator were both cleanec using “3-D Supreme Cleaning 
Agent”, a commercially available product which leaves no residue. All components were 
rinsed with fresh water, dried, then given a final cleaning with acetone. The auxiliary 
condenser coils were cleaned and decontaminated in the same manner. 

The copper condenser tubes also appeared somewhat discolored. These were cleaned 
using “Shower Power”, a commercially available cleaning agent. Prior to installation into 
the condenser unit, they were given a final cleaning with acetone to ensure removal of 
any residue left over from the cleaning procedure. 

During the reassembly, all gaskets were replaced. Prior to replacement of the 
evaporator unit, a test blank was installed on the condenser and a vacuum test success- 
fully performed. After this test was completed, the test blank was removed and the 
evaporator reinstalled. All penetrations were reconnected, the boiling tube bundle 
reinstalled and another vacuum test performed. 

Prior to charging the system with R-113, the ethylene glycol/water coolant pumps 
were energized and flow started through the test condenser tubes. For this initial ‘test, 
smooth copper tubes were used. The temperature of the coolant was -15°C. Since 
Mabrey [Ref. 2] had noted that during condensation, a thickening between the fins of 
copper/nickel finned tubes was occurring which suggested ice formation, it was felt that 
water vapor might be present in the test platform since it had been open to the atmos- 
phere. Therefore, by running the coolant mixture at this temperature through the 
condenser tubes, if no condensation was seen, this would rule out the presence of water 
vapor. During this intitial run, no condensation was seen on any tubes. In addition, the 
orange peel appearance noted by Mabrey on the smooth copper tubes was not in evi- 
dence. 

After checking the evacuated system for possible contamination, the test platform 
was charged with fresh R-113. Coolant flow was again started. This time an orange peel, 
or dimpled appearance could be seen on the top condenser tube. This suggested some 


type of contamination laver. 
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Power (~ 20 kW) was applied to the heaters while the coolant temperature was 
raised to ~ 0°C. The orange peel appearance on the top tube was still present. Meas- 
urements of the vapor temperature at the top and bottom of the condenser showed ap- 
proximately a 20 °C differential, signifying the presence of significant amounts of 
noncondensible gases. These gases enter with the refrigerant and are released into the 
condenser when the liquid 1s boiled. 

A vacuum pump, connected to the top of the condenser was used to remove these 
noncondensible gases. When this was done, the temperature differential between the top 
and bottom of the condenser decreased to less than 0.5°C. The orange peel appearance 
quickly disappeared and all four condenser tubes demonstrated good filmwise 
condensation. Based on these observations, it was standard practice to degas the system 
for at least four hours after the system was open to the atmosphere or fresh R-113 was 
added to the evaporator. 


IV. EXPERIMENTAL PROCEDURES 


A. GENERAL PROCEDURES 

The experimental procedure was identical for all experiments. The coolant pumps 
were energized and flow started through one of the auxiliary condensers. Flow was also 
started through the four test condenser tubes at a nominal coolant velocity of 0.2 m/s. 
The auxiliary heaters were set at 3 KW, the instrumented tube bundle set at 3.5 kW, and 
the simulation heaters set at 1.5 kW. 

The vapor temperature at the bottom and the top of the condenser was allowed to 
rise to 47.5 °C which corresponds to the saturation temperature of R-113 at atmospheric 
pressure. The rate of temperature rise was controlled to approximately 0.3 degrees per 
minute by adjusting the rate of coolant flow through the auxiliary condenser. During the 
heat up period, a vacuum pump, connected to the top of the condenser via the cold trap, 
Was run continuously in order to remove any noncondensibles. The absence of 
noncondensibles was indicated when the vapor temperature difference between the top 
and bottom of the condenser unit was less than 0.5 °C. 

When the vapor temperature was about 45 °C, coolant velocity was adjusted to be- 
tween 1 and 1.2 m/s for the four condenser tubes. The heater power levels were then 
readjusted in order to maintain a relatively constant condenser pressure while leaving 
one auxiliary condenser coil on-line. At lower coolant velocities, condenser pressure was 
maintained using the auxiliary condenser coils. The power settings for the heaters were 
not altered during the course of the experiment. 

During all experimental runs, the vapor temperature at the top and bottom of the 
condenser unit was maintained between 46.5 and 48.5 °C which corresponds to a 
condenser pressure of approximately one atmosphere. Pressure control. was accom- 
plished by adjusting coolant flow through the auxiliary condenser tubes. 

I. Tube Bundle 
Once the vapor temperature was between 46.5 and 48.5 °C , coolant flow was 
adjusted to the desired coolant velocity for the four test condenser tubes. Here, the de- 
sired coolant velocity is an input to the computer program. The program then calculates 
the required flowmeter setting using the calibration curves (Appendix A). The initial 
coolant velocity used was 0.2 m/s. After setting the flow rate, the temperature difference 


between the inlet and outlet of all four condenser tubes was inonitored using the com- 
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puter program, DRPCON7. Approximately ten minutes was required for this temper- 
ature differential to stabilize. Once stable, a data set was taken. This procedure was 
repeated in increments of 0.1] m/s up to a maximum of 1.2 m/s. The procedure was then 
repeated by decreasing the flow in 0.1 m/s incrernents down to 0.2 m/s. Two data sets 
were obtained at each coolant velocity. It should be mentioned that for the 
copper/nickel finned tubes, the range of flow settings was changed. The minimal flow 
setting was 0.6 m/s while the maximum flow setting was 2.0 m/s. The reasons for this 
change are discussed in “3. Enhancement of Outside Heat Transfer Coefficient.” on 
page 30. 
2. Single Tube Experiments 

After completion of the tube bundle runs, the experiment was repeated for each 
individual tube starting with the lowest tube (tube D). Flow through the tubes above the 
test tube was stopped and flow through the tube of interest was set for a velocity of 0.2 

/sec. The vapor pressure was again maintained using an auxiliary condenser coil. 
During this time, the upper tubes were continuously observed for condensate formation. 
No single tube data was obtained until condensate formation on the upper tubes had 
ceased. This typically required twenty to thirty minutes. 

Once the coolant temperature rise had stabilized, a data set was obtained. Fol- 
lowing this, the coolant flow rate was increased in 0.1 m/s increments up to the maxi- 
mum flow velocity. A data set was obtained at each flow velocity. Throughout the 
experiment, tubes above the test tube were visually monitored to ensure that no 
condensate formation occurred. The flow velocity was then decreased in 0.1 m/s incre- 
ments down to the minimum coolant flow velocity, with a data set obtained at each flow 


velocity. This procedure was repeated for each test condenser tube in the bundle. 


B. SPECIFIC EXPERIMENTS 
1. Baseline Experiments 
The purpose of the baseline experiments was to validate the performance of the 
condenser/evaporator test platform by comparing data to that of Mabrey [Ref. 2]. These 
initial baseline experiments were performed using smooth copper tubes (inside diameter 
of 13.26 mm, outside diameter of 15.88 mm, thermal conductivity of 386 W/m K). The 


condensing length of the tubes was 1.219 m. The tubes were supplied by Wolverine Tube 
Co. 
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2. Enhancement of Inside Heat Transfer Coefficient 

The main purpose of this study was to obtain accurate values of the outside heat 
transfer coefficient, A,. In order to obtain accurate values for h,, the resistance to heat 
transfer on the inside of the tube = needs to be minimized. The reason that this re- 
sistance must be minimized is that if the inside resistance is a large part of the overall 
vapor-to-coolant resistance, then /, will be become inaccurate. In order to minimize the 
inside resistance, one must enhance the inside heat transfer coefficient. 

The inside heat transfer coefficient is influenced by the type of coolant flow 
through the tubes. In the test platform, the coolant 1s a 54/46 percent (by weight) 
ethylene glycol/water mixture. At temperatures below 0 °C, this fluid is extremely 
viscous with Reynolds numbers in the range of 300 to 2800, suggesting laminar flow. 
When the coolant flow is laminar, the inside heat transfer coefficient is low. However, 
if this laminar flow can be mixed, then the inside resistance can be decreased and the 
heat transfer coefficient increased. This can be accomplished through the use of radial 
mixing elements. Radial mixing elements have the effect of mixing the thermal boundary 
layer (formed on the inside of the tube) with the colder core flow thereby improving heat 
transfer. It should be noted that a penalty in . .ng such elements is an increased pres- 
sure drop which must always be kept in mind in practical applications. However, for the 
purposes of this study, increased pressure drop was not important. 

In this thesis, two types of mixing elements were investigated. The first was a 
twisted tape insert made of copper [Ref. 47,48]. This tape insert had a thickness of 0.559 
mm with a pitch for a 180 degree twist of three times the inner tube diameter. 

The second type of insert used in this study was a commercially available prod- 
uct referred to as HEATEX radial mixing elements.3 These elements are manufactured 
from stainless steel and consist of a number of small flexible loops (petals) of steel wire 
attached to a central wire core. These flexible petals are placed in the tube such that they 
face the oncoming flow. Their outside diameter is slightly larger than the inside diameter 
of the tube and are therefore kept in position with a “pinch” fit. The petals disturb the 
thermal boundary layer and promote fluid mixing. Figure 4 shows the two types of in- 
serts. [t should be noted that two different sizes of HEATEX were used since condenser 


tubes of two different inside diameters were examined. 


3 manufactured by CAL GAVIN Birmingham, England. 
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Figure 4. 
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As part of the baseline experiments, the smooth copper tubes were also tested 
with no insert installed. Experiments were repeated on the smooth copper tubes using 
the twisted tape inserts as well as the HEATEX elements. This data was then used as a 
baseline for comparison with other types of condenser tubes. It should be mentioned 
that when no insert is present, the inside tube resistance to heat transfer is very large 
relative to the outside resistance. In some cases (typically when the coolant flow rate is 
low), this can produce a negative value for the outside heat transfer coefficient. This is 
simply the inaccuracy in the calculation of A, when the inside resistance dominates and 
clearly demonstrates the need for an insert as mentioned above. 

3. Enhancement of Outside Heat Transfer Coefficient. 

The overall heat transfer coefficient, U,, 1s influenced by the outside as well as 
the inside heat transfer coefficient. The means by which the outside heat transfer coeffi- 
cient 1s enhanced 1s by altering the outside geometry of the condenser tube. 

Experiments to investigate enhancement of the outside heat transfer coefficient 
utilized three types of condenser tubes. The first type of tube tested was a copper/nickel 
finned tube (26 fins/inch). This tube has an inside diameter of 10.16 mm with a thermal 
conductivity of 43 W/m K. The tubes were manufactured by Wolverine Tube Co. These 
tubes were tested as a tube bundle as well as individually. Experiments were performed 
with no insert as well as with the twisted tape and HEATEX inserts. These tubes are 
shown in Figure 5. One point must be mentioned with regard to the finned tubes. - The 
use of fins on the outside of the tube enhances the heat transfer considerably. Conse- 
quently, the inside resistance again becomes a large fraction of the overall resistance, 
even when inserts are used. 

One additional point must be discussed with regard to the copper/nickel finned 
tubes. Because the size of the holes in the nylon end plates is fixed, the outside diameter 
of anv tubes to be tested must be the same. In order to have these finned tubes the same 
outside diameter, the inside diameter had to be reduced. Since the inside tube diameter 
was smaller for the finned tubes, higher coolant velocities were obtained. 

The second series of experiments used a copper/nickel KORODENSE tube also 
manufactured by Wolverine Tube Co. This tube, shown in Figure 5, has a shallow spiral 
groove along the entire condensing length with a pitch of 10 mm. As with the finned 
copper/nickel tubes, experiments were performed on these KORODENSE tubes with no 
insert as well as with the twisted tape and HEATEX inserts. The tubes were tested in a 


tube bundle as well as individually. 


Fujii and co-workers [Ref. 33] have previously shown that the condensation of 
refrigerants can be enhanced by wrapping smooth condenser tubes with a small diameter 
wire. Therefore, a third series of experiments were perforined in which copper/nickel 
KORODENSE tubes were tightly wound with wire as shown in Figure 5. Since these 
authors showed that there was an optimum wire size which gave the greatest enhance- 
ment of heat transfer (0.5 mm at a pitch of 3 mm), the KORODENSE tubes used in this 
study were wrapped with one of three different sizes of wire. The sizes used were 0.029” 
(0.736 mm) diameter stainless steel wire, 0.049” (1.244 mm) diameter stainless steel wire, 
and 0.0675” (1.715 mm) diameter copper wire. The wire lay in the shallow grooves of the 
KORODENSE tubes and was soldered to the tubes at the ends only. Although the wire 
was tightly wrapped around the tubes, it was not physically attached except at the ends. 
Experiments were carried out on each tube individually. Each tube was tested in all four 
locations within the bundle. This was done to see if position within the bundle (i.e. wake 
from tubes above) had any effect. Only the HEATEX inserts were used with these tubes. 

Once an optimal wire diameter was determined, a fourth series of experiments 
was conducted. This series used the same four KORODENSE tubes, only this time each 
was wrapped with the same optimal diameter wire. The tests were done with the 
HEATEX elements installed. Experiments were conducted on the tubes individually in 
order to compare with the previously obtained single tube data. In addition, these wire 
wrapped tubes were tested as a tube bundle in order to investigate whether or not the 


wire wrap had any effect on condensate innundation. 
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V. DATA REDUCTION 


A. DATA ACQUISITION 

The data acquisition system consisted of a Hewlett-Packard series 9000 computer 
coupled with a Hewlett-Packard model 3479A Data Acquisition system. The computer 
program used in this study was originally written by Mabrey [Ref. 2] but has been sub- 
stantially modified in the course of this thesis. The data acquisition thermocouple 


channel assignments are shown in Table 1. 


Table 1. THERMCOUPLE CHANNEL ASSIGNMENTS. 


0 | Vapor (condenser top) 
[1 [Vapor (condenser top) 
[3 tiguia (top of bundiey 
[a Tiquie (cop of bundtey 
[~s__[Liguid (bottom of poo 
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The computer program allows the user to acquire data, reprocess data, or plot data. 
During data acquisition, the user selects the fluid velocity desired. Based on coolant inlet 
temperature, the coolant fluid properties and mass flow rates are calculated (see Ap- 
pendix C). The program then uses equations derived from the flowmeter calibrations 
(Appendix A) in order to calculate the appropriate flowmeter settings. 

Once the flowmeters are set for the required coolant velocity, the computer system 


displays the vapor temperature at the top and bottom of the condenser as well as the 


a2 


coolant temperature rise for each condenser tube. When these temperature rises are 
stable (approximately 10 minutes), the user can terminate the display and acquire the 
data for that coolant flow velocity. While the temperature rises Were stabilizing, coolant 
flow through an auxiliary condenser tube was adjusted in order to maintain constant 


pressure. 


B. DATA REDUCTION 

The main function of the data reduction section of the computer program is to back 
out the outside heat transfer coefficient from an experimentally determined value of U, 
and an assumed inside heat transfer coefficient correlation. This is accomplished in the 
following manner. 

The data acquisition system measures the emf from the thermocouples listed in Ta- 
ble 1. These values of emf are converted to temperatures using correlations obtained 
from the temperature calibrations both contained in Appendix B. The volume flow rate 
of the ethylene glycol/water mixture is calculated based on the desired coolant velocity. 
The density, dynamic viscosity, thermal conductivity, heat capacity, and Prandt] number 
of the coolant are calculated using inlet temperatures (see Appendix C and Appendix 
E). From these quantities, the coolant mass flow rate and fluid velocity are calculated. 
The coolant Reynolds number for each condenser tube is also calculated. The heat 


transferred to the coolant is then found from 
q = mC, AT (5 —1) 


where m is the mass flow rate and AT is the coolant temperature rise. The outside heat 
flux, gq”, 1s calculated by dividing g by the surface area of the tube based on the outside 
tube diameter. 
The overall heat transfer coefficient, U,, is given by: 
Gee 


lees One 


where LMTD is the log mean temperature difference defined as: 
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The inside heat transfer coefficient (given as a dimensionless Nusselt number, Nu,) 
can then be calculated using one of three correlations. The choice of which correlation 
to use depends upon the type of insert installed inside the condenser tube. If no insert 
is installed, then Nu, can be calculated from a correlation provided by Hausen [Ref. 


49}4 for developing laminar flow in a tube at constant wall temperature. 
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If a twisted tape insert is installed, a correlation given by Hong and Bergles [Ref. 


Nap = 3.6541 OS) 


47; can be used to calculate Nu, where 
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Here, y is the twist ratio (pitch over diameter). Re, 1s the coolant Reynolds number, 
corrected for tape thickness. Re, 1s given by: 


4m 


Nes = ald — 45) 


(S — 6) 
where 6 is the tape thickness. This correlation is valid for Re, up to about 2300. 

A number of correlations for determining the inside heat transfer coefficient have 
been published when using a wire matrix insert [Ref: 50,51,52, 53 ]. However, these 
correlations are specific for the type of mixing element used, and are especially depend- 
ent upon the loop density of the wire matrix. 

The HEATEX mixing elements supplied by CAL GAVIN (nominal diameters of 
13.26 and 10.16 mm) did not match the loop densities used by previous investigators. 
However, CAL GAVIN did supply experimental data for the heat transfer coefficient 
versus fluid velocity at three different temperatures for both sizes of elements. To con- 


vert this data into a usable correlation, the following procedure was developed. 


4 Calculations were performed for coolant flow through the smooth copper tubes at coolant 
SOUS of 0.2 and 1.2 m/s. These calculations confirmed the coolant flow was not fully devel- 
oped. 
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The correlations of Oliver and Aldington [Ref. 51, 52 } took the form of 
Nu = cRe”" Pr" (5—7) 


where c, m, and n are arbitrary constants which are dependent upon the type of insert 
used. It was assumed that a similar type of equation could be used in the present study 
by fitting the constants to the data (once non-dimensionalized) provided by CAL 
GAVIN. 

The correlations of Oliver and Aldington (Ref. 51, 52 ] had Prandtl number raised 
to a power of 0.46. This was primarily to take into account temperature effects. There- 
force, the same value was assumed In the present study. Values of Nu/Pr°“ were plotted 
as functions of Re on a log-log basis. Figure 6 shows the data obtained for the 13.26 
mm elements while Figure 7shows the data obtained for the 10.16 mm elements. Based 
on the appearance of the non-dimensionalized data, a least squares linear regression was 
performed on each data set. The figures yield the values of the constants m and c (slope 
and intercept respectively) and provide the following correlations for the 13.26 and 10.16 
mm clements respectively: 


Nu, = 0.226Re”®? Pr? (5 38) 
and 
Nu, = 0.065Re”° pr?" (5 —9) 


Once Nu, was obtained from equations (5-4, 5-5 5-8 or 5-9), the inside heat transfer 


coefficient was calculated from: 
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Figure 6. HITRAN correlation for 13.26 mm elements 
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HEATEX CORRELATION (10.15 MM) 
NU = 0.067(RE**.76)(PR**.46) 
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Figure 7. HITRAN correlation for 10.16 mm elements 
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The overall heat transfer coefficient, U, 1s found from the summation of resistances 


to heat transfer across the whole tube: 


l I 
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h,A; 
where R, and R, are the fouling and wall resistances respectively. In order to find 4, 


equation (5-11) can be rearranged to give: 
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In all experiments, the tubes were meticulously cleaned and therefore, the fouling re- 
sistance, R,, was neglected. 

In order to better assess the calculated experimental value of the outside heat 
transfer coefficient, some form of comparison with prediction must be made. In this 
case, the experimentally obtained value of h, was compared with Nusselt theory. How- 
ever, because the wall temperatures are not known a priori, an iteration must be used 
to determine this predicted value of h,. An initial value of h, was estimated and based 
On this estimate, a wall temperature was calculated using the equation: 
qo 
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The condensate film properties are then calculated using a reference temperature given 
by: 
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Once these condensate fluid properties are known (Appendix D), the theoretical outside 


heat transfer coefficient can be calculated from Nusselt theory [Ref- 4]: 
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This calculated value is then compared to the initial guess. The error between the initial 


guess and the calculated value are compared. If the error is greater than 0.1%, the 
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“guessed” value for A, is updated and the iteration continued. This process is repeated 


until the error between the calculated and guessed values is less than 0.1%. 


C. MODIFIED WILSON PLOT 

The outside and inside heat transfer coefficients could be determined directly by 
measuring the vapor temperature, the tube wall temperature, and the coolant temper- 
ature rise. Of these, the most difficult temperature to measure directly is the tube wall 
temperature. This would require the use of instrumented tubes with thermocouples em- 
bedded in the tube wall. Unfortunately, the manufacture of such instrumented tubes is 
both costly and time consuming, especially when large numbers of tubes are to be tested. 

An alternative to the use of instrumented tubes is to employ the “Modified Wilson 
Plot” procedure. This procedure has been outlined in detail by Marto [Ref. 10]. The 
modified Wilson plot procedure relies upon the overall heat transfer coefficient being 
reliably measured from experimental data. A summation of heat transfer resistances 
from vapor-to-coolant is then expressed as: 
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As mentioned above, the fouling resistance is assumed to be negligible. This summation 


of resistances can be rearranged into a linear equation: 


I oe | a ae | Gea 
U Macy A; h; hy 


In looking at the correlations given above for h, and h, it can be seen that they all 
employ a leading coefficient. As an example, Nusselt theory, (equation 5-15) can be ex- 


pressed in the form: 
h, = oF (S — 18) 
where: 
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In a sinular manner, the equations for the inside heat transfer coefficients (equations 5-4, 


5-5, 5-8, and 5-9) can also be transformed into the form: 
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h, = CO (5 — 20) 


where the form of 2 would depend on which type of tube insert was being used. As an 
example, using equation 5-5, 22 would become: 


k R 125 10,5 
Q=— [ + s.asaxio™ Pro" Fe | (5-21) 


The linearized version of the summation equation (equation 5-17) can then be rear- 
ranged to give: 
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This is a linear equation with two unknowns, namely, C anda. Letting 








y=| : ~ RyAg F (5 — 23) 
U, 
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we can obtain a simplified linear equation: 


Y=mX+6 (5 — 25) 
where 
moos (5 — 26) 
es 
and 
l 
b= (5 — 27) 


Since 22 and F are both temperature dependent, an iterative solution procedure must be 
used. This is accomplished by applying a least squares fit to the linear equation. 

It should be pointed out that the accuracy of the modified Wilson plot procedure 
relics heavily on the number of data points taken as well as the range of coolant veloci- 
ties used. 
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VI. RESULTS AND DISCUSSION 


A summary of the experimental runs performed in this thesis is presented in 
Table 2. The abbreviation “SM” is for the smooth copper tubes, “KD” is for the 
KORODENSE tubes, and “CF” is for the copper/nickel finned tubes. 


Table 2. SUMMARY OF EXPERIMENTAL RUNS 


[—swios | Copper | HEATEX | —SS—SSCS 
[sor [Copper | Tape | SSCS 
[~KDo? | KORODENSE [None [SS 
[~KDos | KORODENSE | HEATEX |S 
[—KDor | KORODENSE | Tape [SSS 
Peres cwnifin [Tope] SSCS 
Pr Gros cwri fin” | HEATEX | SSS 
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A. DECONTAMINATION 

Prior to charging the svstem with R-113, the system was evacuated to approximately 
15 inHg and left for 24 hours. For this initial test, smooth copper tubes were used. The 
ethylene glycol coolant pumps were energized and flow started through the test 
condenser tubes. The temperature of the ethylene glycol was -15°C. No condensation 
was seen on any tubes. In addition, the orange peel appearance noted by Mabrey [Ref. 
2} was not observed. 

After checking the evacuated system for possible contamination, the test platform 
was charged with fresh R-113. Coolant flow was again started. No orange peel, or 
dimpled appearance could be seen on the top condenser tube. 

Power (~ 20 kW) was applied to the heaters while the coolant temperature was 
raised to ~ 0°C. An orange peel appearance was noted for the top tube. Measurements 
of the vapor temperature ai the top and bottom of the condenser showed approximately 
a 20 °C differential. 

A vacuum pump, connected to the top of the condenser, was used to remove any 
noncondensible gases. When this was done, the temperature differential decreased to less 


than 0.5°C. The orange peel appearance quickly disappeared from the top condenser 


tube. This suggests that at least part of the contamination problem reported bv Mabrev 


was due to the presence of noncondensibles. 
In a subsequent series of experiments conducted on a set of titanium finned tubes, 


a foreign substance appeared to build up on the tubes during the experimental run. This 
occurred despite continuous operation of the vacuum pump. The tubes were removed 
from the condenser and inspected. The foreign material had an oily feel to it and an odor 
simular to that of ethylene glycol. 

The titanium finned tubes had smooth extensions (6” in length) attached to their 
ends using epoxy. These extensions penetrated the nylon block and steel plate and pro- 
vided the smooth surface necessary for the O-rings to give an airtight seal between the 
condenser and the atmosphere. When these tubes were removed from the condenser, it 
was noted that the extensions on two of the tubes had come loose. The condenser was 
opened for inspection. Upon removal of the end plate and shroud assembly, pools of 
ethylene glycol were seen inside the condenser shell and on the auxiliary condenser tubes. 
In addition, ethylene glycol was found around the shroud assembly. The entire 
condenser and evaporator assembly were therefore dismantled and cleaned in accordance 


with the procedures outlined in “Experimental Apparatus”. 
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After reassembly of the condenser and evaporator, the smooth copper tubes were 
cleaned and reinstalled. Operation of the test platform failed to show any orange-peel 
or dimpled appearance on the tubes. In addition, all tubes demonstrated good filmwise 
condensation. 

It should be mentioned that in subsequent tests, copper/nickel finned tubes were 
installed in the condenser. However, unlike the titanium finned tubes, the smooth ex- 
tensions on these tubes were soldered. During system operation, no foreign material 
was noted between the fins. Condensate retention between the fins was negligible. The 


tubes demonstrated good film wise condensation. 


In summary, it is believed that the contamination problems noted by Mabrey [Ref. 
2] were due to two problems. The first problem was the presence of noncondensibles 
during operation of the condenser. This resulted in the relatively low values for overall 
heat transfer as well as the fact that the top tube, which would be most affected bv the 
noncondensibles, had the lowest heat transfer. Secondly, while it can not_be proved 


conclusively, the author believes that ethviene glycol mav also have been present in the 
svstem for certain tubes and that this is what was the material noted in between the fins 


bv Mabrev. 


B. OVERALL HEAT TRANSFER COEFFICIENT 

In the following sections, the results pertaining to the overall heat transfer coeffi- 
cient are presented. It is important to keep in mind that this value is experimentally de- 
termined and does not rely on any correlations other than for the fluid properties. In 
addition, it should be mentioned that in the following graphs, tube A is the top tube in 
the bundle, tube B the second, tube C the third, and tube D the bottom tube. 

1. Smooth Copper Tubes 

Figure § shows the overall heat transfer coefficient, U,, plotted as a function 

of coolant velocity. The tubes did not contain any inserts. Note the low values for U,. 
These values are almost constant for velocities over the range of 0.2 to 0.7 m/s. In this 
series of experiments, Reynolds numbers were typically of the order of 300 to 2300 in- 
dicating that flow is laminar. Without any insert to promote mixing, one would there- 
fore expect relatively low values for U, . This is because, in laminar flow, the inside 
resistance to heat transfer is relatively large and therefore leads to a lower value for heat 
transfer. Finally, there is essentially no difference between the tubes at low coolant ve- 
locities. This suggests that with these relatively low condensation rates, inundation is not 


a significant factor. At higher flow rate, it can be seen that the values for heat transfer 


start increasing and that differences between the tubes become more apparent. As 
coolant velocity is increasing, flow is beginning to become turbulent. This reduces the 
inside resistance to heat transfer which leads to an increase in the overall heat transfer. 
As heat transfer increases, the effects of inundation become more apparent. This ts pre- 
sumably due to the outside resistance becoming more important as the inside resistance 
becomes less important. 

The results obtained when these smooth copper tubes were tested individually 
is shown in Figure 9. Here, all the tubes behave similarly. In addition, comparison with 
Figure 8 shows that the values of U, for the tubes when tested individually are approx- 
imatelv the same as for the top tube during the bundle run as one would expect. 

The results for the same run with tape inserts installed are shown in Figure 10. 
It can be seen that U, increases as coolant velocity increases. The values are clearly 
greater than when no insert 1s used (approximately 50% greater) indicating the advan- 
tage gaimed when enhancing the inside heat transfer coefficient in situations where 
coolant flow is laminar. The effect of condensate inundation on the overall heat transfer 
coefficient is also clearly seen. Here, condensate inundation results in approximately a 
15% decrease in U, when the top tube 1s compared to the second tube. Interestingly, 
there appears to be little difference in the values for U, for tubes B, C, and D indicating 
that the condensate film is of similar thickness for these tubes. The data for a similar 
run performed by Mabrey [Ref. 2] is shown in Figure 11. The data obtained in the 
present study is almost twice the values reported by Mabrey. In addition, the present 
data shows that the top tube has the highest heat transfer which is to be expected as- 
suming noncondensibles are not present. This was not the case in Mabrey’s work. 

The results for these same tubes when operated individually are shown in 
Figure 12. All the curves are essentially the same with values for U, corresponding to 
that obtained for the top tube during the bundle run. The maximum variation is about 
12%! 

The next series of experiments were performed on the smooth copper tubes with 
the HEATEX elements installed. Figure 13 shows the results of these experiments when 
the tubes were operated as a bundle. As with the experiments performed with the 
twisted tape inserts, there is an increase in U, as coolant velocity is increased. In this 
case the increase is 100% in comparison to the case with no insert! Also, as in the pre- 
vious experiment, condensate inundation appears to result in about a 15% decrease in 
U, for the second tube compared to the top tube. There appears to be little difference 


in the heat transfer coefficients for tubes B, C, and D, although some effect may be ar- 
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Figure 8. Bundle run on smooth copper tubes without inserts. 
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Figure 9. Single tube run on smooth copper tubes without inserts. 
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Figure 10. Smooth copper tubes with tape insert operated as a bundle. 
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Figure 12. Smooth copper tubes with tape insert operated individually. 
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gued at the higher coolant velocities. This may reflect the fact that, at the higher coolant 
velocities, the inside resistance decreases while the outside resistance increases in im- 
portance. Anything which causes further increases in this resistance, such as condensate 
inundation, may cause noticeable effects in heat transfer. 

As before, these same tubes with the HEATEX inserts installed were tested in- 
dividually; the results are shown in Figure 14. The values of U, for all four tubes are 
virtually identical and match the values for the top tube obtained during the bundle run 
(Figure 13). The values appear to be about 10% greater than for the tape insert ex- 
periments indicating that the HEATEX provides better tube-side enhancement. 

While the preceding figures have shown the effect of inserts and also inundation 
on the heat transfer coefficient, two important questions must be addressed. First, what 
is the uncertainty associated with the calculated values for U,. Secondly, the repeat- 
ability of the data must also be examined. 

An uncertainty analysis was performed on U, as shown in Appendix F based 
on the procedures of Kline and McClintock [Ref. 54]. At the low coolant velocities, the 





uncertainty was 5.2%. This uncertainty was based primarily on the uncertainty in the 
mass flow measurement. Here, one must include a scale interpolation term as well asa 
term to account for the time-wise jitter for the flowmeter. Since the flow rate is low for 
this coolant velocity, these two terms become more significant than at the higher flow 
rates. This is especially true if the magnitude of the time-wise jitter is the same for all 
flow rates which was the case 11 this study. 

At high coolant velocities, the uncertainty in the value for U, was about 5.6%. 
Here, the major factor driving this value is the uncertainty in the determination of tem- 
perature. The reason for this is that the AJ between the tube inlet and outlet is smaller 
than for low coolant velocities. Since the relative uncertainty for the thermocouples re- 
mains the same regardless of coolant velocities, uncertainty in the teniperature meas- 
urement, specifically the LMTD, becomes the driving factor in determining the 
uncertainty in U,. 


The uncertainty analvsis performed here was onlv for the overall heat transfer 
coefficient. U,. In this analvsis, anv uncertaintv in the correlations used to determine the 





lvcol coolant or R-113 were neglected. Therefore 
the level of uncertainty reported here is slightly conservative. 

In order to examine the reproducibility of the data, repeat experiments were 
performed on the smooth copper tubes with the HEATEX elements installed. It is im- 


portant to mention that during the time between the two runs, the tubes were removed 
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Figure 13. Smooth copper tubes with HEATEX inserts. Bundle operation. 
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Figure 14. Smooth copper tubes with HEATEX installed. Individual operation. 
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from the condenser, cleaned and reinstalled. The results for the duplicate experiment are 
shown in Figure 15. When these values are compared to Figure 13, differences of the 
order of 3% 1s seen in the values of U, for the top tube. U, for the lower tubes are slightlv 
larger with a more pronounced difference between the tubes. However, these differences 
all lie within the 6% uncertainty band indicating that the system gives good 
reproducibility. 

In a similar manner, the results for a repeat experiment on the smooth copper 
tubes tested individually with the HEATEX installed is shown in Figure 16. When the 
values for U, shown here are compared with those shown in Figure 14, no significant 


differences are seen. Again, this demonstrates the good reproducibility obtained by this 


system. However, a word of caution 1s in order. In all the experiments reported here, 


great care has been taken to control as many variables as possible. Of particular im- 
portance is to ensure that the heater power levels are set at the same level for each ex- 


periment. This will tend to keep vapor velocity in the condenser at similar levels for each 
experiment which will result in better repeatability. 

The effects of having no insert as opposed to the twisted tape or HEATEX in- 
sert 1s shown in Figure 17. In this figure, the values for U, for the top tube only have 
been plotted as a function of coolant velocity for the smooth copper tubes. As previously 
noted, without any insert, coolant flow is laminar and JU, is relatively constant up to 
about 0.7 m/s. When either the twisted tape or HHEATEX inserts are installed, U, in- 
creases markedly as coolant velocity increases. Interestingly, the HEATEX results in 
about a 10 to 14% increase in U, when compared to the twisted tape inserts. Hence, 
from a heat transfer standpoint, the HEATEX element is the better insert. This is 
probably because the wire loops which make up the insert, touch much of the inside 
surface of the tube thereby “scrapping” the thermal boundary laver from the wall. With 
the twisted tape insert, mixing 1s provided only by inducing swirling which 1s probably 
not as an effective mechanism. 

2. Copper/Nickel KORODENSE Tubes 

The results obtained for the overall heat transfer eoellicient,. U. for 
KORODENSE tubes operated as a bundle are shown in Figure 18. This series of ex- 
periments was performed without any inserts installed. Note the similarity between these 
results and those obtained under similar conditions with the smooth copper tubes (Fig- 
ure 8). Again, values of U, appear relatively constant for coolant velocities up to 0.7 


m,'s. As coolant velocity increases further, values of U, increase dramatically. In addi- 
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Figure 15. Repeat run on smooth copper tubes with HEATEX insert. 
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Figure 17. Effect of tube inserts on U, 
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tion, an effect of condensate inundation begins to be apparent at these higher coolant 
velocities. 

Results obtained when these same tubes were operated individually are shown 
in Figure 19. Again, when there is no inundation, all tubes have similar behavior. In 
addition, the values obtained here match those obtained for the top tube during the 
bundle experiment (see Figure 18). Interestingly, the relative magnitude of U, is re- 
markably similar to that obtained for the smooth copper tubes (Figure 9). This 1s 
somewhat surprising. The KORODENSE tubes are corrugated on both the inside as well 
as the outside. Therefore, one would have expected somewhat higher values for heat 
transfer with the KORODENSE than with the smooth copper tube. 

Values of U, for the KORODENSE tubes with the twisted tape inserts installed 
are shown in [igure 20. This data was obtained during bundle operation. As coolant 
velocity increased, U, increased steadily. This was true for all four tubes. The effect of 
inundation reduced U, by about 15% (tube B relative to tube A). There appeared to be 
more separation in the values of U, for tubes B, C, and D than was the case for the 
smooth copper tubes (see Figure 10) suggesting a greater inundation effect. In compar- 
ison with the data obtained for the smooth copper tubes, the KORODENSE tubes re- 
sulted in a slight increase (about 6 to 8%) in U, This difference may be due to the fact 
that the insert reduces the inside heat transfer resistance so that the accuracy of the 
measurement is improved. 

The results obtained for these same tubes operated individually are shown in 
Figure 21. As expected, all four tubes behave similarly when no inundation is present. 
In addition, the values for U, match those obtained for the top tube during the bundle 
operation (Figure 20). In comparison to the smooth copper tube data (Figure 12), the 
data for the KORODENSE tubes appears somewhat higher (about 6 to 8%). However, 
given the known uncertainty, this difference is probably not significant. 

Results obtained for the KORODENSE tubes during bundle operation with the 
HEATEX elements installed are shown in Figure 22. Note that U, again increases 
markedly for all four tube as coolant velocity increases. An inundation effect can be seen 
for tubes B, C, and D which is similar to that seen in the twisted tape runs. Tube B 
shows approximately a 15% decrease in U, compared to tube A. This effect is magnified 
at higher coolant velocities (condensation rates) as one would expect. When this data 1s 
compared to that obtained for the smooth copper tubes with HEATEX installed (see 
Figure 13), it can be seen that U, is greater for the KORODENSE tubes. This suggests 


that roping of the tube provides some enhancement for the outside (and possibly the 
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Figure 18. KORODENSE tubes with no insert, bundle operation. 
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Figure 19. KORODENSE tubes operated individually, no inserts installed. 
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Figure 21. KORODENSE tubes with tape inserts, individual operation. 
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inside) heat transfer coefficients. While we expect some small difference, it must be 
pointed out that the differences observed here are at the bounds of the calculated un- 
certainty. 

Results for these tubes tested individually are shown in Figure 23. As before, 
the values of U, are similar for all four tubes and match the values obtained for the top 
tube during the bundle experiments (see Figure 22). In comparison to the data obtained 
for the smooth copper tubes, the values of U, are somewhat higher as expected, but are 
still at the bounds of the calculated uncertainty. 

The effect of using no insert, twisted tape and HEATEX on U, can be seen in 
Figure 24. The data are for the top tube only. Note that, as with the smooth copper 
tubes (Figure 17), the HEATEX results in about a 10 to 14% increase in U, in com- 
parison to the twisted tape insert. Both provide significant enhancement in U, when 
compared to no insert (almost 100% enhancement). As with the smooth copper tubes, 
when no insert is used, U, appears relatively constant for coolant velocities over the 
range of 0.2 to 0.7 m/s indicating laminar flow. 

3. Wire Wrapped KORODENSE Tubes 

This series of experiments was performed to determine an optimum wire diam- 
eter as discussed in “Experimental Procedure”. All tubes were tested individually and all 
had HEATEX mixing elements installed. Figure 25 shows U, as a function of coolant 
velocity for the case of an ordinary KORODENSE tube together with tubes wrapped 
with 0.029”, 0.049” and 0.0675” diameter wire, each tested in the top position (tube A). 
The 0.049” diameter wire gives the biggest enhancement in U, (about 25%). While the 
other two wire diameters also provide enhancement in U, (about 15%), it is below that 
of the 0.049” wire. The no wire case agrees with the previously obtained data (see 
Piguines2 2). 

Figure 26 shows similar data but in this case for tube B. Again, the 0.049” di- 
ameter wire provides the greatest enhancement in U, (again about 25%). 

Figure 27 shows what happens when the different diameter wires are used on 
tube C. In this case an interesting fact can be observed. Here, both the 0.049” and 0.029” 
wire provide virtually the same degree of enhancement for U, (about 20%) with both 
being better than the 0.0675” wire. Figure 28 shows the results obtained for tube D. 
Note that while the 0.049” wire again provides the greatest degree of enhancement for 
U, (about 16%), the differences between the three different diameter wires are not as 
great. 
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Figure 23. KORODENSE tubes with HEATENX inserts, individual operation. 
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Figure 24. Effect of tube inserts on U, 


65 


UZERO WIRE WRAPPED KORODENSE 
TUBE A, R-113, HEATEX 


UZERO (W/M**2 K) 
900.0 1000.0 1100.0 1200.0 


800.0 


LEGEND 
o= NO WIRE 
o= 0.029" WIRE 
4= 0.049" WIRE 
+ = 0.0675" WIRE 


700.0 





600.0 


0.7 0.8 0.8 1 i.1 1.2 1.3 


0.1 0.2 0.3 0.4 0.5 0.6 
VELOCITY (M/SEC) 


Figure 25. Effect of wire diameter on U, for tube A. 
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Figure 26. Effect of wire diameter on U, for tube B. 
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Figure 27. Effect of wire diameter on U, for tube C. 
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Figure 27 and Figure 28 suggest that the position of the wire wrapped tube 
within the bundle may affect the overall heat transfer coefficient. In order to gain further 
insight on this, the data obtained from these experiments was replotted so that U, for a 
given wire diameter was plotted as a function of coolant velocity for all four tube posi- 
tions. The data for the tube without any wire 1s shown in Figure 29. The data 1s all very 
similar and no distinct trends are evident. Note that comparison of this data with 
Figure 23 again demonstrates the reproducibility of the system. 

Figure 30 (0.029” wire diameter) shows that tubes A, B and C behave similarly. 
Tube D produces consistently lower values for U, than the other tubes. Figure 31 
(0.049” wire) demonstrates what appears to be a positional effect on U,. The data for 
tube A and B is virtually identical and lies above that of tubes C and D which are also 
virtually identical. Figure 32 suggests that there is also a positional effect for the 
0.0675” wire. The top tube seems to provide the greatest enhancement in U,. The level 
of enhancement for tubes C and D are approximately the same, both being lower than 
tubes A and B. 

Figure 33 shows U, plotted as a function of wire diameter for an arbitrarily 
chosen coolant velocity of 1.0 m/s. Note that the strongest positional affect appears to 
be for the 0.049” wire. Also, it is clearly evident that the 0.049” wire provides the greatest 
enhancement for heat transfer. 

The possible positional effect may be due to vapor velocity. Although the: ve- 
locity of the vapor was very small in this study (approximately 0.1 m/s) there may still 
be an effect. As the vapor penetrates the tube bundle, it slows down. As its velocity is 
reduced, it can no longer help to strip condensate from the tubes. Hence, heat transfer 
is decreased. This effect looks much like an inundation effect. An interesting point is 
why this shows up in these experiments and not with the plain KORODENSE or 
smooth copper tubes. With the wire wrapped tubes, the condensate layer is presumably 
thin. Therefore, anything which could help to thin it further would be of importance in 
improving heat transfer. The loss of this small amount of vapor shearing may be im- 
portant in this case. 

In this study, the 0.049” diameter wire appears to give the best enhancement in 
heat transfer (see Figure 33). This yields a pitch to wire diameter ratio of 7.2. Fujii et 
al. [Ref. 33] found that for condensation of R-12 on wire wrapped smooth copper tubes, 
the optimal pitch-to-wire diameter ratio was on the order of 2. This difference could be 
due to our use of the KORODENSE tubes. For a pitch-to-wire diameter ratio of 2, we 


would have had to use wire which was of the order of 5 mm in diameter. This would have 
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Figure 29. Effect of no wire on U, for all tubes. 
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Figure 30. Effect of 0.029” wire on U, for all tubes. 
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Figure 31. Effect of 0.049” wire on U, for all tubes. 
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Figure 32. Effect of 0.0675” wire on U, for all tubes. 
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Figure 33. U, as a function of wire diameter for the four tubes. 
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been unrealistic for our experiments. It may be that the optimal pitch-to-wire diameter 
may in fact be about 2. Additional experiments need to be performed with different 
pitth KORODENSE tubes. 

4. KORODENSE Tubes Wrapped with 0.049” Wire 

Once the optimal wire diameter had been determined, two identical experiments 
were performed on the KORODENSE copper/nickel tubes wrapped with the 0.049” di- 
ameter wire. The reason that two experiments were performed was to check the repeat- 
ability of the data. Figure 34 and Figure 35 show the results from these two experiments 
for bundle operation of the tubes. The two experiments show excellent repeatability. In 
comparing these results with the non-wire wrapped tubes (Figure 22), several interesting 
points can be made. The values of U, for all four tubes are greater for the wire wrapped 
tubes than for the plain tubes (approximately 16% for the top tube). In addition, the 
inundation effect (seen as a lowering of U,) is reduced for the wire wrapped tubes. 
Finally, note that the differences between the tubes are more apparent than for the 
non-wire wrapped tubes. This suggests that the wire wrapping may have an effect upon 
condensate inundation. This will be addressed in more detail later. 

Figure 36 and Figure 37 show the results obtained when the tubes wrapped 
with the 0.049” diameter wire were operated individually. All the tubes behave similarly 
and demonstrate good reproducibility between the two runs. 

5. Copper/Nickel Finned Tubes 

Figure 38 shows the results obtained when no insert was installed and the tubes 
were operated as a bundle. Comparison of these results with those shown in Figure 8 
for the smooth tubes shows that the finned tubes yield higher heat transfer coefficients 
as expected (approximately 12%). Also, as with the smooth copper tubes, the top tube 
again displays the highest value of U,, due to the fact that the lower tubes suffer from 
condensate inundation. The single tube data for the finned tubes is shown in 
Figure 39. All four tubes are virtually the same as the top tube during bundle operation. 

Figure 40 and Figure 41 show the results obtained for the finned tubes with the 
tape insert installed during bundle and individual tube operation respectively. In bundle 
operation, the values of U, are nearly double those obtained for both the smooth copper 
tubes (Figure 10) and the KORODENSE tubes (Figure 20) using similar inserts. As 
with the smooth and KORODENSE tubes, inundation again causes a reduction in U,. 
This reduction for the second tube appears to be about 10%. Comparison of the single 


tube operation (Figure 41) with the bundle operation (Figure 40) shows that all four 
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Figure 34. 0.049” wire wrapped KORODENSE tubes, bundle operation, expt. 1. 
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Figure 35. 0.049” wire wrapped KORODENSE tubes, bundle operation, expt. 2. 
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Figure 36. 0.049” wire wrapped KORODENSE tubes, individual operation. 
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Figure 37, 0.049” wire wrapped KORODENSE tubes, individual operation. 
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Figure 38. Copper/nickel finned tubes, no insert, bundle operation. 
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Figure 39. Copper/nickel finned tubes, no insert, individual operation. 
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tubes operated individually have approximately the same value of U, as the top tube 
during the bundle operation. 

Data obtained from experiments on the copper/nickel finned tubes with the 
HEATEX elements installed are shown in Figure 42 and Figure 43 for bundle and in- 
dividual tube operation respectively. Note that in comparison with the tape insert ex- 
periments (Figure 40 and Figure 41), use of the HEATEX inserts results in about a 
20% increase in U, which again indicates the improvement afforded by using this type 
of insert. Again the inundation effect results in a decrement of 10 to 14%. During in- 
dividual tube operation, all four tubes again behave like the top tube during the bundle 
Operation. 

As with the smooth copper tubes and the KORODENSE tubes, the effects of 
using different types of inserts can be summarized. This has been done in Figure 44. 
This data is for the top tube only. The similarity with data from the other tubes (see 
Figure 17 and Figure 24) is good with HEATEX yielding approximately a 20% increase 
in U, over the twisted tape insert and a remarkable 200% increase over the case with no 


insert! 


C. OUTSIDE HEAT TRANSFER COEFFICIENT 

The following figures deal with the outside heat transfer coefficient, A,. In all cases, 
except as specifically noted, the outside heat transfer coefficient was determined using 
the modified Wilson Plot procedure. As before, tube A is the top tube with tube D being 
the bottom tube in the bundle. In addition, no values were obtained for h, for exper- 
iments which had no insert. The reason for this is that, for the coolant velocities used, 
without inserts, the inside resistance was simply too large for any accurate results to be 
obtained for h,. 

1. Smooth Copper Tubes 

| Figure 45 and Figure 46 show the values obtained for A, for the smooth copper 
tubes with the tape insert installed during bundle and individual tube operation respec- 
tively. Several features deserve comment. First, values of h, for the top tube are below 
those predicted by Nusselt theory. These differences range from about 15% at low 
coolant velocities to about 7% at the higher coolant velocities. In addition, note the 
inundation effect (seen as a reduction in h, for tubes B, C, and D). Reductions in A, range 
from 28% to 19% as coolant velocity increases from 0.2 to 1.2 m/s respectively. As was 


the case for the overall heat transfer coefficient, when the tubes were operated individ- 
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Figure 40. Copper/nickel finned tubes, tape insert, bundle operation. 
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Figure 41. Copper/nickel finned tubes, tape insert, individual operation. 
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Figure 42. Copper/nickel finned tubes, HEATEX insert, bundle operation. 
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Figure 43. Copper/nickel finned tubes, HEATEX insert, individual operation. 


87 


KFFECT OF INSERTS 
CU/NI FINNED TUBES 


LEGEND 
So = NO INSERT 
Oo = TAPE INSERT 
4 = HEATEX INSERT 


UZERO (W/M**2 K) 


400.0 600.0 800.0 1000.0 1200.0 1400.0 1600.0 1800.0 2000.0 2200.0 





06 O7 O8 O88 10 11 #42 423 «214 «126 £426 «17 «180 «1.9 862.0 
VELOCITY (M/SEC) 


Figure 44. Summary of effect of tube inserts on U, 
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ually, all had values for A, which were virtually identical to those obtained for the top 
tube during bundle operation. 

Figure 47 and Figure 48 show values for f, obtained when the HEATEX ele- 
ments were installed in the copper tubes during bundle and individual tube operation 
respectively. As was the case for the tape insert, values of h, range from about 10 to 
14% below that predicted from Nusselt theory. Again, there is a substantial inundation 
effect which results in a 15 to 18% decrease in h, when comparing the top tube to the 
second tube. Note that the values for h, obtained for the individual tubes agree with 
those obtained for the top tube during bundle operation. Finally, comparison of the 
value for h, obtained when the tape inserts Were used as opposed to the HEATEX ele- 
ments reveal no effect on h, which is to be expected since we are examining the outside 
heat transfer coefficient. 

It is interesting to note that Goto et al. [Ref. 5] found that values of A, for 
R-113 condensing on long, smooth copper tubes were about 10% below that predicted 
by Nusselt theory. Marto et al. (Ref. 8 ] found that the vapor-side heat transfer coeffi- 
cient for R-113 condensing on short (0.133m) smooth copper tubes was in relatively 
good agreement with Nusselt theory. This result was subsequently confirmed by Michael 
et al. (Ref. 32]. However, it is important to point out that both the latter reports used 
very short condenser tubes. The tubes used by Goto et al. were about 0.6 m in length, 
while the tubes used in this study had a condensing length of 1.2 m. Hence tube length 
may be an important consideration. The reason for this is as follows. The coolant en- 
ters the condenser tube with a given inlet temperature. As the coolant flows through the 
tube, the coolant temperature increases. The wall termperature is also changing down 
the length of the tube. At the tube entrance, the temperature difference between the tube 
wal and the vapor saturation temperature causes the formation of a thick condensate 
layer. This layer is thinner at the tube outlet due to the smaller wall-vapor saturation 
temperature difference. This causes condensate to flow in the axial direction. These fac- 
tors taken together imply that the temperature difference is a function of the length of 
the tube. This factor is not taken into account in the Nusselt analysis. As coolant ve- 
locity increases the difference betwee the coolant inlet and outlet temperatures is re- 
duced. The wall temperature becomes more constant along the tube length so that values 
for h, are in better agreement with those determined experimentally. 

The effect of using the twisted tape and HEATEX insert on the inside and out- 
side heat transfer coefficients is demonstrated in Figure 49 and Figure 50 respectively. 


In both cases, the values for A, are virtually identical as expected. However, the 
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Figure 45. Smooth Copper tubes, tape insert, bundle operation. 
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Figure 46. Smeoth Copper tubes, tape insert, individual operation. 
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Figure 47. Smooth Copper tubes, HEATEX insert, bundle operation. 
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Figure 48. Smooth Copper tubes, HEATEX insert, individual operation. 
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HEATEX elements essentially double the inside heat transfer coefficient. This is what 
provides the improvement in overall heat transfer coefficient when the HEATEX ele- 
ments are used. 

The effects of inundation were examined by calculating an average and a local 
heat transfer coefficient. The local h, is found by dividing the value for h, for a given 
tube by the value for the top tube. Values obtained are then compared with Nusselt 
theory (Eq. 2-18) and the Kern model (Eq. 2-20). The average heat transfer coefficient 
is determined by taking an average value of h, for the tube of interest plus the tubes 
above it and dividing this value by 4, for the top tube. Values obtained are then com- 
pared to Nusselt theory (Eq. 2-17), the Kern model (Eq. 2-19), or the Eissenberg model 
(Eq. 2-12). For both the local and average values, the flatter the curves, the less the ef- 
fect of condensate inundation. 

The effects of inundation on A, for the smooth copper tubes are shown in 
Figure Sl, Figure 52, Figure 53 and Figure 54 at three arbitrarily chosen coolant ve- 
locities. As expected, there is virtually no difference on the local or average values of the 
heat transfer coefficient when the tape insert is used as opposed to the HEATEX ele- 
ment. In addition, the values for the average heat transfer coefficient lie in between the 
values predicted by the Kern and Eissenberg models (equations 2-19, and 2-21). Values 
for the local heat transfer coefficient lie slightly above those predicted by the Kern model 
(equation 2-20). These results suggest that Nusselt theory is too conservative with re- 
gard to the effect of condensate inundation. It 1s interesting to note that the main 
inundation effect occurs for the second tube with little additional effect for subsequent 
tubes in the bundle. This again supports the notion that the smooth copper tubes can 
only support a limited amount of condensate 

2. KORODENSE Tubes 

Figure 55 and Figure 56 show values for h, for the KORODENSE tubes with 
the tape inserts installed. The figures show data obtained when the tubes were operated 
as a bundle and individually respectively. The values of h, for the top tube are above that 
predicted from Nusselt theory which is due to the small enhancement provided by the 
roping of the tubes. This enhancement is approximately 15 to 20% in comparison to the 
sniooth copper tubes. As with the smooth copper tubes, there is again a substantial ef- 
fect of inundation (on the order of 20% for the second tube compared to the top tube). 
There is also more separation between the values of 4, for tubes B, C, and D than for 


the smooth copper tubes indicating that the lower tubes may not be “saturated” with 
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Figure 49. Effect of tape insert on A, and hf, 
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Figure 50. Effect of HEATEX insert on kh, and A, 
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Figure 51. Local A, for tape insert. 
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Figure 52. Local h, for HEATEX element. 
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Figure 53. Average A, for tape insert. 
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Figure 54. Average h, for HEATEX element. 
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condensate. Finally, note that the relative values for A, for the KORODENSE tubes are 
about 10 to 15% higher than the corresponding values for the smooth copper tubes. 

The data for the KORODENSE tubes operated individually show much more 
scatter than was seen for the smooth copper tubes especially at low coolant flow rates. 
Also, rather than being virtually identical with those obtained for the top tube during 
bundle operation, the values of h, lay scattered between those for the top tube (bundle 
operation) and the values predicted by Nusselt theory. Presumably, this is due to either 
inaccuracies or uncertainty in the inside correlation especially at low coolant flow rates. 
At higher coolant velocities, the scatter appears to be within about 10%. 

Values for h, for the KORODENSE tubes obtained when the HEATEX ele- 
ments Were installed are shown in Figure 57 and Figure 58 for bundle and individual 
tube operation respectively. As with the experiments conducted with the tape inserts 
installed, values of h, again he above those predicted by Nusselt theory. The enhance- 
ment appears to be approximately 15 to 20%. Again there ts a substantial inundation 
effect which is of the order of 20% for the second tube in comparison to the top tube. 
The data obtained for the individual tube operation shows less scatter. However, values 
for all the tubes again le between those predicted by Nusselt theory and those obtained 
for the top tube during bundle operation. 

As with the smooth copper tubes, the effects of changing inserts on the inside 
and outside heat transfer coefficients are shown in Figure 59 and Figure 60. The inside 
heat transfer coefficients are identical to those of the smooth copper tubes (Figure 49 
and Figure 50). It was felt that the enhancement provided by the roping would be much 
less than the enhancement provided by the inserts. Therefore, the same correlations 
were used to determine the inside heat transfer coefficient. The values for A, are slightly 
greater for the KORODENSE tubes than for the smooth copper tubes (approximately 
5 to 10%). 

Figure 61 and Figure 62 show the local heat transfer coefficient for the four 
tubes for the tape insert and the HEATEX elements. There 1s virtually no difference 
between the two. The values lie above the Kern estimate. Figure 63 and Figure 64 show 
the average heat transfer coefficient for the four tubes with both inserts. The differences 
for values using the tape inserts as opposed to the HEATEX elements are again negli- 
gible. The average values for h, also lie just above the Kern model. This is in contrast to 
the smooth copper tubes where the average values lay near the Nusselt estimate. This 
suggests that the effect of condensate inundation is indeed less for the KORODENSE 
tubes. 
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Figure 55. KORODENSE tubes, tape insert, bundle operation. 
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Figure 56. KORODENSE tubes, tape insert, individual operation. 
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Figure 87. KORODENSE tubes, HEATEX insert, bundle operation. 
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Figure 58. KORODENSE tubes, HEATEX insert, individual operation. 
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Figure 59. Effect of tape insert on , and hk, 
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Figure 60. Effect of HEATEX insert on A, and fh; 
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Figure 61. Local A, for tape insert. 
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Figure 62. Local 4, for HEATEX element. 
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Figure 63. Average A, for tape insert. 
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Figure 64. Average h, for HEATEX element. 
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3. Wire Wrapped KORODENSE Tubes 

The effect on A, of wrapping wire around the KORODENSE tubes is shown in 
Figure 65, Figure 66, Figure 67 and Figure 68. In viewing these figures, one should 
remember that the values for the unwrapped KORODENSE tubes were about 10% 
greater than predicted from Nusselt theory and about 15% greater than for the smooth 
copper tubes. From the figures, it appears that the 0.049” diameter wire results in the 
greatest enhancement in A, for all four tube positions. However, the relative magnitude 
of the enhancement decreases as lower tubes are examined. Since this data was obtained 
during individual tube operation, the results are not due to condensate inundation. 

Values for h, are shown for each wire diameter in all four tube positions in Fig- 
ure 69, Figure 70, Figure 7! and Figure 72. When no wire is used, the results are the 
same as those shown in Figure 58. For the 0.029” wire, positions A and B are essentially 
identical with tube C showing the greatest enhancement and tube D the worst. The 
magnitude of the deviation is about 20%. The 0.049” diameter wire results in the best 
enhancement in A, regardless of position. However, note that position D is not over- 
whelmingly greater than for the 0.029” wire. The data show that the level of enhance- 
ment depends upon bundle position. The top tube positions behave similarly while 
positions C and D show reduced levels of enhancement Finally, the 0.0675” diameter 
wire yields values for h, which are nearly the same although a positional effect appears 
evident. An interesting comparison are the values for h, for the 0.049” and 0.0675” wire 
in the fourth position. Here it appears that the 0.0675” wire gives roughly the same level 
of enhancement especially at high coolant velocities. 

4. KORODENSE Tubes Wrapped with 0.049” Wire 

The effects on h, of wrapping 0.049” diameter wire on KORODENSE tubes on 
are shown in Figure 73 and Figure 74. These two experiments are identical runs. In 
both cases, there is a substantial improvement in A, for all tubes when compared to val- 
ues predicted from Nusselt theory. It is interesting to note that even though tubes B, 
C, and D are all being inundated, the values for h, are still well above Nusselt theory 
(32% for tube A, 28% for tube B, 22% for tube C and 16% for tube D). In comparison 
with the non-wire wrapped KORODEWNSE tubes, wrapping with the 0.049” diameter 
Wire vields an increase of about 70% for tubes A and B and about a 30 to 50% increase 
for tubes C and D. Comparison with the top smooth tube yields almost a 90% increase 
in A,. 

The results for the single tube runs are shown in Figure 75 and Figure 76. Both 


runs demonstrate considerable scatter in the data of about + 22%. It is also interesting 
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Figure 65. Effect of wire diameter on A, for tube A. 
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Figure 66. Effect of wire diameter on hf, for tube B. 
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Figure 67. Effect of wire diameter on A, for tube C. 
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Figure 68. Effect of wire diameter on A, for tube D. 
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Ficure 71. Effect of 0.049” wire on A, for all tubes. 
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Figure 72. Effect of 0.0675” wire on A, for all tubes. 
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Figure 73. 0.049” wire wrapped KORODENSE tubes, bundle operation, expt. 1. 
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Figure 74. 0.049” wire wrapped KORODENSE tubes, bundle operation, expt. 2. 
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to note that the data for these single tube runs is somewhat higher than for the bundle 
runs (of the order of 11%). This may be due to uncertainty in the correlations used for 
h, which are reflected in A, 

The values for the average and local vapor-side heat transfer coefficients are 
shown in Figure 77, Figure 78, Figure 79 and Figure 80 for the two experiments con- 
ducted with the 0.049” wire wrapped KORODENSE tubes. The average vapor-side heat 
transfer coefficient is equal to that predicted by the Eissenberg model whereas the local 
vapor-side heat transfer coefficient is substantially above the Kern estimate. Recall that 
for the plain KORODENSE tubes, the average heat transfer coefficient was above the 
Kern estimate but below that predicted by the Eissenberg model. The local heat transfer 
coefficient was just above the Kern estimate. The differences between the plain and 
wire-wrapped KORODENSE tubes demonstrate that wrapping the tubes with wire 
substantially reduces the effects of condensate inundation. 

5. Copper/Nickel Finned Tubes 

The values for h, obtained from experiments on copper/nickel finned tubes using 
the tape inserts are shown in Figure 81. The values obtained for the top tube are ap- 
proximately 7 times greater than predicted by Nusselt theory. Condensate inundation 
caused about a 35% decrease in h, for the second tube. 

The values for h, obtained using the same tubes but with the HEATEX elements 
are shown in Figure 82. In comparison to the values obtained for the tape inserts (see 
Figure 81), this data appears very scattered with enhancement over Nusselt theory of 
the order of 10 to 15 fold. This is a large discrepancy. It is possible that the correlation 
for the small diameter HEATEX element may be suspect. The data supplied by CAL 
GAVIN was collected on smooth copper tubes with a diameter simular to that of our 
copper/nickel finned tubes. This data was collected over a lower coolant velocity range 
than that used in this study. Therefore, the correlation for this HEATEX element may 
not be valid for the finned tube over the range of coolant velocities used here. This area 
needs to be studied further, perhaps even to gathering new experimental data which can 
be used to develop a new correlation. 

The degree of inundation can be seen in Figure 83, Figure 84, Figure 85 and 
Figure 86. The data for the average vapor-side heat transfer coefficient for both 
HEATEX as well as tape inserts lie between the Nusselt and Kern models. This ts also 
true for the local vapor-side heat transfer coefficients. 

An interesting comparison is to look at the average and local heat transfer val- 


ues for all tubes tested. Regarding the average heat transfer coefficient, for the smooth 
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Figure 75. 0.049” wire wrapped KORODENSE, individual runs, expt. 1. 
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Figure 76. 0.049” wire wrapped KORODENSE, individual runs, expt. 2. 
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Figure 77. Average h, wire wrapped KORODENSE, expt. 
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Figure 79. Local h, wire wrapped KORODENSE, expt. 1. 
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Figure 80. Local h, wire wrapped KORODENSE, expt. 2. 
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Figure 81. A, for copper/nickel finned tubes, tape insert. 
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copper tubes, the values lay between the Kern and Eissenberg estimates. The plain 
KORODENSE tubes had values similar to those for the smooth copper tubes while the 
0.049” wire-wrapped KORODENSE had values which matched estimates based on 
Eissenberg’s model. The copper/nickel finned tubes, on the other hand, had values which 
were between the Kern and Nusselt estimates. Since the lower the values the more 
condensate effect present, this data would imply that the effect of inundation is worse 
on the finned tubes and least on the wire-wrapped tubes. If the data for the local heat 
transfer coefficients is analyzed in a similar manner, the same conclusions are reached. 

The wire wrapped KORODENSE tubes have the least effect of condensate 
inundation presumably because the wire serves to create a low pressure area at the base 
of the wire. This causes condensate to be pulled along the tube which enhances heat 
transfer. The values for the finned tubes are not in agreement with those of Honda et 
al. [Ref. 42]. These investigators found that the degree of inudation effect was substan- 
tially less than that reported here. There are two apparent reasons for this. The pitch- 
to-diameter ratio used by Honda et al. was 1.37 whereas in the present study it 1s 2.25. 
This difference may cause differences in the way condensate drains from one tube onto 
another. Alternatively, the correlations used to determine the heat transfer coefficient 
may be suspect. Hong and Bergles [Ref. 47] have reported that their correlation for the 
twisted tape has an uncertainty of about 15%. Therefore, it is possible that the large 
effect of condensate inundation is due more to uncertainty in the correlations than toa 
real phenomenon. 

In order to assess the validity of the data obtained on the copper/nickel finned 
tubes, a comparison was made with the model of Beatty and Katz [Ref. 13]. In order to 
do this comparison, the fin efficiency had to be determined using equations from 
Incropera and DeWitt [Ref. 55]. It was assumed that the fins could be approximated 
as rectangular fins with adiabatic tips. The fin efficiency was then approximated by: 

tanh(mL,) 


7 mL, ioe 


According to Incropera and DeWitt, the term mL, can be approximated by [Ref. 55}: 


| er, 
2h_\z,2 
mbe= (2 L, (6 — 2) 
Pp 
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Figure 83. Average h, for Cu/Ni finned tubes, tape insert. 
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Figure 85. Local A, for Cu/Ni finned tubes, tape insert. 
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where L, 1s defined as: 


L.=L+ aay 


== 
2 
Here, L is defined as the height of the fin and t is the fin thickness. The term 4, is de- 
fined as: 


A, = Lt (6 — 4) 


In order to calculate a value for the fin efficiency, 7,, an iterative procedure was 
used. First, a value for A, was guessed. Based on this value, 7, was calculated using 
equation (6-1). An effective area was then calculated from the following equation (Marto 
(Ref. 3]): 


where A, 1s the surface area of the fins sides, A, is the surface area of the fin tips and 
A, 1s the interfin surface area. Using these values, the effective diameter, D, was deter- 
mined from (Marto (Ref. 3]): 


me io2s a As 1 An | Ay ot 
Pal = 1.30n, Aep [025 tg ae pos Ar Ay pos (cr 0) 


where L is defined as: 


_ 2(D,~D/) 
L= 4D, (6 — 7) 
The heat transfer coefficient is then calculated from a modified version of the Beatty and 
Katz model (Eq. 2-2) which uses heat flux rather then the film temperature difference, 


iol ig 


kp’ ghye > 
h, = 6oss{ De 


ez Oo 


(6 — 8) 


The results of these calculations for the copper/nickel finned tubes with the tape 


inserts are shown in Figure 87. In this figure, the values of h, predicted from the Beatty 


and Katz model have been expressed in terms of the smooth tube area, that is zd@,L . 
Also note that values for h, have been calculated for copper, copper/nickel and titanium. 


It is important to note that the raw data lies below the predicted values for copper/nickel 
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finned tubes. The data obtained after reprocessing with the modified Wilson plot proce- 
dure lies somewhat above the prediction. However, both lie within + 20% which is 
within the range of validity for the Beatty and Katz model. Therefore, it would appear 
that the data obtained in the present study for the copper/nickel finned tubes with the 
tape inserts are reasonable. 

The results for the copper/nickel finned tubes with the HEATEX elements in- 
stalled are not as good. These results are shown in Figure 88. This figure shows the 
tremendous amount of scatter in the data irrespective of whether the raw data or the 
data obtained from the modified Wilson plot procedure is used. In either case, it is ap- 
parent that the data do not match values predicted from the model of Beatty and Katz. 
The only difference in these two experiments is the type of insert used. The values for 
h, are based on an experimentally determined value for U, and a value for h, which is 
obtained from a correlation. In this study, the correlations for the HEATEX elements 
were based on empirical data supplied from the manufacturer. This data was obtained 
on smooth copper tubes having the same diameter as the copper/nickel finned tubes used 
in this study. Another point which must be considered is that the range of velocity for 
which the empirical data was supplied was from 0.2 to 1.2 m/s. However, for the 
copper/nickel finned tubes used in this study, a coolant velocity of 0.7 to 2 m/s was used. 
These factors would suggest that the correlation used for the small diameter HEATEX 
element 1s in error. Further research will have to verify whether or not this is the case. 


D. IMPLICATIONS OF PRESENT WORK | 

This study used R-113 as the refngerant and a refrigerated ethylene glycol/water 
mixture as the coolant. However, Naval refrigeration systems use R-114 for the 
refrigerant and sea water as the condensing medium. Therefore, one must ask whether 
the results obtained from the present work can be extended to Naval applications. In 
order to perform the following analysis, the dimensions of the smooth copper tubes and 
copper/nickel finned tubes (26 fpi) used in this study were assumed. The inside heat 
transfer coefficient was calculated using the Dittus-Boelter correlation. For the smooth 
tubes, A, was calculated using the Nusselt equation, modified to use heat flux rather than 
the film temperature difference, AT . For the finned tube, h, was again calculated using 
the Dittus-Boelter correlation. However, /f, was calculated using the modified Beatty and 
Katz equation (Eq. 6-8). The fin efficiency was assumed to be 1.0. The condenser design 
parameters were modeled after the requirements for the DDG-51. Hence, the refriger- 
ation load was fixed at 242.7 tons. The number of condenser tubes was fixed at 340 
tubes. Coolant inlet temperature was assumed to be 65 F (18.3 C) while the saturation 
temperature of R-114 was 107.8 F (42.1 C). The coolant velocity was held constant at 
2m/s. 
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Figure 87. Experimental data (tape insert) vs. Beatty/Katz model. 
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Ratios of the surface areas for finned to smooth tube and a finned KORODENSE 
tube to smooth tube were calculated. In order to carry out this calculation for the 
KORODENSE tube, some type of correlation must be used to calculate h,. However, 
as shown by Cunningham and Ben Boudinar [Ref. 56] in their work on steam 
condensation on roped tubes, the enhancement due to the roping depends upon the type 
of groove cut into the tube. In this analysis, an enhancement of 1.6 for h, was assumed 
for the KORODENSE tube over the finned tube. Once the area ratio is obtained, a 
condenser tube length can be calculated. This is done as follows. The area ratio can be 
expressed in terms of length by: 


oy NayLy 
Ays  Npid,L, 





(6 — 9) 


where A,,1s the area of the finned tube, A,, 1s the area for the smooth tube, and N is the 
number of tubes. The term a, is given by: 


ap = (Ap + AN (6 — 10) 


where 4, and A, are the finned and unfinned areas respectively. 

Table 3 shows the results of this calculation. For the finned tube, the area ratio of 
1.75 leads to a decrease in the condenser length of about 50%. On the other hand, if a 
roped KORODENSE tube could be manufactured which had fins on the outside surface, 
the reduction in length of the condenser tubes would be 58%. This could lead to even 
greater Weight savings for refrigeration systems. 


Table 3. SUMMARY OF AREA RATIO CALCULATIONS 


Smooth Finned 
Aus Ae 
Ae 







ty 1.0 0.497 - 9.42 
iby 
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VII. CONCLUSIONS 
The following conclusions can be made based on the results obtained in this study. 


GENERAL 


. The condenser/evaporator test platform is fully operational and shows no signs of 


the contamination problems reported by Mabrey [Ref. 2 ]. 


Values for outside heat transfer coefficient obtained using smooth copper tubes 
appear to be reasonable and agree with other published data. 


OVERALL HEAT TRANSFER COEFFICIENT 


. KORODENSE tubes provide enhancements of 6 to 8% over smooth copper tubes 


regardless of the type of insert used. 


. For the KORODENSE tubes used in this study, wire wrapping results in an en- 


hancement in heat transfer. 


. The optimal pitch-to-wire diameter was about 7 for the KORODENSE tubes used 


in the present work. 


. Wrapping the KORODENSE tubes with 0.049” diameter wire resulted in an en- 


hancement of 16% over the non wire-wrapped KORODENSE tubes and a 22 to 
24% enhancement over the smooth copper tubes. 


. Copper/nickel finned tubes (26 fpi) yielded a doubling of the heat transfer in com- 


parison to smooth copper tubes. 


. The use of HEATEX radial mixing elements yielded a 10 to 15% increase in com- 


parison to the use of twisted tape inserts. 


OUTSIDE HEAT TRANSFER COEFFICIENT 


. The smooth tubes yielded values of k, which were 7 to 15% below those estimated 


from Nusselt theory. 


. The KORODENSE tubes gave values of A, which were about 15% greater than 


those estimated from Nusselt theory and approximately 25% greater than values 


obtained for the smooth copper tubes. 


. Wrapping the KORODENSE tube with 0.049” diameter wire resulted in a 90% 


increase in the outside heat transfer coefficient in comparison to the smooth copper 
tubes. 


. The copper/nickel finned tubes yielded a 7 fold increase in A, over the smooth 


copper tubes. 


EFFECTS OF CONDENSATE INUNDATION 


. For smooth copper tubes, condensate inundation resulted is a 15 to 18% decrease 


in h, for the second tube compared to the top tube. 


142 


2. KORODENSE tubes showed a 20% decrease in A, for this same comparison. 


3. Wrapping the KORODENSE tubes with 0.049” diameter wire resulted in a decrease 
of only 5 to 10% in hk, for the second tube compared to the top tube. 


4. For the copper/nickel finned tubes, condensate inundation resulted in almost a 
35% decrease in fA, for the second tube relative to the top tube. This large decrease 
is presumably due to uncertainty in the correlations used to determine the heat 
transfer coefficient. 
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VII. RECOMMENDATIONS 
Based on the results of this study, the following recommendations are made: 


A. GENERAL RECOMMENDATIONS 


1. The ball valves which control coolant flow to the flowmeters should be replaced 
with needle valves. This will provide better control of coolant flow. 


2. The flowmeters should be replaced with larger capacity meters. This will allow in- 
vestigation of heat transfer at higher coolant velocities and therefore higher 
condensation rates. This change will have to be accompanied by an increase in the 
diameter of the coolant return header in order to accommodate the increased 
coolant flow. 


3. Extra lines from the top ends of the condenser to the cold trap need to be added. 
This will ensure that no isolated pockets of noncondesibles exist in the condenser. 


B. EXPERIMENTAL RECOMMENDATIONS 


1. he correlation for the small diameter HEATEX element must be re-investigated. 
This may require obtaining additional data from CAL GAVIN or designing an 
experiment at the Naval Postgraduate School which can be used to obtain the re- 
quired data. 


2. The correlation used for the twisted tape insert also needs to be re-evaluated in 
light of the large inundation effects. As with the HEATEX, this may require that 
additional experiments be carried out at the Naval Postgraduate School. 


3. Once a reliable correlation is obtained, the experiments conducted on the 
copper/nickel finned tubes need to be repeated. 


4. Experiments need to be conducted using both titanium finned tubes as well as 
titanum KORODENSE tubes. This data can then be compared to that obtained 
with the copper/nickel finned tubes and the copper/nickel KORODENSE tubes to 
determine whether titanium tubes can deliver the required performance. 


5S. Smooth copper tubes should be wrapped with different wire diameters at different 
pitches in order to determine an optimal pitch-to-wire diameter ratio. 


6. Experiments need to be conducted on wire wrapping of KORODENSE tubes which 
' have a different pitch. These experiments should be based on results obtained from 
the smooth copper tube experiments. 


7. Experiments need to be repeated using R-114 and possible alternate refrigerants 
that the Navy is interested in such as R-24. 


8. A KORODENSE type tube which is roped on the inside and finned on the outside 
needs to be developed and tested since calculations have shown that this type of 
tube can lead to smaller condensers. 


144 


APPENDIX A. FLOWMETER CALIBRATION 


The flowmeters used to measure flow through the main condenser tubes were cali- 
brated in the following manner. The meters were disconnected from the condenser tubes. 
An empty 55 gallon drum was placed on a scale (resolution 0.5 Ibs). The condenser 
pump was started with flow through the bypass. 

Flowmeters were calibrated at 0 °C and 24 °C. Each flowmeter was initially opened 
to maximum flow. The flow was then throttled back to 10%. The weight of fluid in the 
drum was recorded and the timer started. The weight was then recorded at the end of the 
calibration time. This procedure was repeated in 10% increments up to 100% flow. 
Flowmeter number 2 was only calibrated up to 90% flow since a mechanical stop pre- 
vented it from going to 100%. The calibration runs were then repeated from the maxi- 
mum flow back to 10% flow in 10% increments. The mass flow rate was then calculated 
in kg/sec based on an average of the two readings. The data for the flowmeter cali- 


brations as a function of nominal flow rate are shown in Table 4. 


Table 4. MASS FLOW RATE (KG/SEC) AT << ea AND 24 C. 


Sl Flow Meter A Flow Meter B Flow | Flow Meter C | C. | Flow Meter D | Meter D 


Nominal 
Flow Rate 
(%) 


SS 
20 [oars | 0339 | ors | 0.295 | ono | oasie | oe | 01823 
[80 | oases] ove | 0.1533 | oars | 01699 | oie | _oa261_| 0.2409 
[so [oniees | 0.2028 | 0.1799 | 0912 | ones | 0197s | 0.2001 | 0.2667 
[10002019 | o2isn |_| | 0204] oars | oa7ee_| 02024 | 
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A least squares linear regression was performed for each flowmeter at each of the 
two temperatures for which the calibration runs were carried out. The linear regression 


equation takes the following form: 
m=a+bN 


where a is the y-intercept in kg/s, b is the slope in kg/s/% and N is the flowmeter setting 


in %. The results are shown in Table 5. 


Table 5. FLOWMETER CALIBRATION REGRESSION RESULTS. 


ea. Flow Meter A Flow Meter B Flow Meter C Flow Meter D 
Tre [me [ac | we | oc | we | oe [mc 


0.00206 | 0.00225 | 0.00204 | 0.00216 | 0.00208 | 0.0022 0.0028 0.0029 
-0.00329 | -0.00615 | -0.00138 | -0.0036 | -0.00211 | -0.00122 | 0.0048 





In order to obtain flowrates at any intermediate temperatures, a simplified interpo- 
lation procedure was used. The results of the calibrations for the four flowmeteres are 


shown in Figure 89, Figure 90, Figure 91, and Figure 92. 
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Figure 90. 
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Figure 91. Calibration for flowmeter C. 
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Figure 92. 
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APPENDIX B. EMF VERSUS TEMPERATURE 


The EMF read from the thermocouples must be converted into a temperature for 
pruposes of calculation. This is accomplished using a 7“ order polynomial which has the 


general form: 
T=C(O)+C(NV' 


where V is the EMF reading in volts, and I has values from | to 7. The leading coeffi- 
cients are given in Table 6. 


Figure 93 shows the relationship between EMF and temperature. 


Table 6. LEADING COEFFICIENTS FOR CONVERSION OF 
EMF TO TEMPERATURE. 
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Figure 93. EMF versus temperature. 


APPENDIX C. PROPERTIES OF ETHYLENE GLYCOL/WATER 
MIXTURE 


The physical and thermodynamic properties of the ethylene glycol/water mixture 
were determined using equations and figures from Cragoe [Ref. 56] and Gallant [Ref. 
oak 

The density of the ethylene glycol/water mixture is calculated in lines 2490 to 2530 
of the computer program. Density as a function of temperature for various concen- 
trations of ethylene glycol are shown in Figure 94. 

The kinematic viscositv of the ethylene glycol/water mixture is calculated in lines 
2370 to 2420. Figure 95 shows the kinematic viscosity plotted as a function of temper- 
ature for various ethylene glycol concentrations. 

The heat capacity is calculated in lines 2430 to 2480. Figure 96is a plot of heat ca- 
pacity as a function of temperature for different ethylene glycol concentrations. 

The thermal conductivity is calculated in lines 2560 to 2610 of the program. Values 
of the thermal conductivity as a function of temperature are plotted for various ethylene 


glycol concentrations in Figure 97. 
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Figure 94. Density of ethylene glycol/water mixture. 
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Figure 95. Kinematic viscosity of ethylene glycol/water mixture. 
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Figure 96. Thermal heat capacity of ethylene glycol/water mixture. 
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Figure 97. Thermal conductivity of ethylene glycol/water mixture. 
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APPENDIX D. PHYSICAL AND THERMODYNAMIC PROPERTIES OF 
R-113 


The physical and thermodynamic properties of R-113 were calculated based on 
equations supplied by Chapman [Ref. 58]. 


The heat capacity, C,, of R-113 was calculated from the equation: 


C, = (929.0 + 1.037) 


where T is in °C and C, is in kJ/kg K. Figure 98 shows C, as a function of temperature 
for R-113. 


The dynamic viscosity for R-113 is calculated from the equation: 


503 
a ane CEBAS, 
u=1.34x10 ~ x 10 


where T is in K and yw 1s given in kg/m s. The relationship between y and temperature 
is shown in Figure 99. 


The thermal conductivity of R-113 can be calculated from: 
k = 1.7308(4.846x107? — T x 6.57x107°) 


where T is in °F and K isin W/m K. Figure 100 shows the thermal conductivity as a 
function of temperature for R-113. 


The density of R-113 can be calculated from: 
o = 1.62x10° — T x (2.219 + T x 2.358x107°) 


where T is in °C and p is in kg/m. Density of R-113 is shown as a function of temper- 
ature in Figure 101. 


The latent heat of R-113 can be calculated from the equation: 


hy = 2326.1(70.55 — t x (0.04838 + 1.26x10~7)) 


where T isin °F and A, is in J/kg. Figure 102 shows the relationship between latent heat 


and temperature for R-113. 
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Figure 98. Heat capacity versus temperature for R-113. 
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Figure 99. Viscosity as a function of temperature for R-113. 
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Figure 100. Thermal Conductivity of R-113 as function of temperature. 
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Figure 101. Density as a function of temperature for R-113. 
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Figure 102. Latent heat versus temperature for R-113. 
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APPENDIX E. SAMPLE CALCULATION 


A sample calculation was conducted on run 1 of data set SM06. This run was carried 
out on smooth copper tubes with the HEATEX insert installed. The sample calculation 
was performed for the top tube only. 

Temperatures were recorded as emf’s using copper/constantan thermocouples. EMF 
values were converted to temperatures using the equation shown in Appendix B. 


Table 7 summarizes the results. 


Table 7. EMF AND TEMPERATURES FOR RUN 1 DATA SET SM06. 


713) 909 
(7 
(CC 





















For tube A, the temperature differential is calculated by: 


T(10) + T(11) 
AT =, — T(6) = 13.086 C 


The average temperature is then calculated from: 


Taye = T(6) + AT x 0.5 = 15.98 C 


Based on the ethylene glycol/water mixture physical properties (Appendix C), and using 
and average temperature of 15.98 C and an ethylene glycol concentration of 54%, the 


following fluid properties are obtained: 


p, = 1071.3 kg/m 

u, = 0.0051 kg/m s 
Cu= S196 29 ike K 
k, = 0.3742 Wm k 
ler 43.00 


Based on the calibration curves for flowmeter A, the mass flowrate of coolant cor- 
responding to a setting of 15.74 can be calculated. A simplified interpolation procedure 
is used to account for temperature differences. Based on this calculation, we obtain a 
mass flow rate of 0.2979 kg/s for this run. 

The coolant velocity can be calculated from: 


ae 
- pA; 





where the inside cross-sectional area is given by: 


32 
A;= 7 4 


A; = al (0.013259)? = 0.000138 m1” 


and 


0.02979 


Ve = 11071.3)(0.000138) 


= 0.201 m/s 
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The coolant Reynolds number, Re,, is calculated from: 


Vd, 


i 
Vo 


Re,.= 





— (0.201)(0.013259) 


ai6oxi0® 


Cc 


The heat transferred to the coolant is calculated from: 


q = mC, AT = (0.02979)(3196.99)(13.085) 


g = 1246.19 W 


and the heat flux is then calculated from: 





do= A 


O 
where 


A, = nd,L = r(0.015875)(1.2192) = 0.0608 m? 


where L is the condensing length of the tube. Based on this equation, we obtain: 


, _ 1246.19 


- os 
Too 0608 7 20494.9 Wim 


In order to calculate the inside and outside heat transfer coefficients, the physical 
and thermodynamic properties of R-113 must be calculated (Appendix D). These fluid 
properties are based on a reference or condensate film temperature. In order to obtain 
this film temperature, the wall temperature is needed. To obtain this, an iteration is 
carricd out with an initial value of h, assumed (in this case 1000 W/m? K). The wall 


temperature is then calculated from: 





qo 
Jans =] at 
il Sat h, 
T watt = 46.925 — == = 26.43°C 
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The film temperature is then calculated from: 


Tq 2 
[p= = = ea 





OT 


_ 46.925 


Ty 3 


+ = (26.43) = 33.26°C 


Based on the equations of fluid properties for R-113 (Appendix D), the following values 


were obtained for the condensate: 


h,, = 147036.5 J/kg 

k = 0.0734 W/m K 

p = 1544.34 kg/m 

u = 6.029x10- kg/ms 


Using these fluid properties together with equation (5-15), we obtain: 
h, = 1907.06 W)m’?K 


Using an iterative procedure, we adjust the initial guess of h, by taking the average of 
the calculated and guessed value of h,. The wall and film temperatures are then recalcu- 
lated and new fluid properties for R-113 are obtained. A new value of h, 1s then calcu- 
lated. This procedure continues until the newly calculated and previosuly calculated 


values are within 0.1%. For this run, we obtained: 
h, = 1251 Wim? K 
The next step is to find a value for the overall heat transfer coefficient, U,, where: 


” 
q o 


Uo=TMTD 


where, for this run, 
LMTD = 30.48°C 


and 
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U, = SS = 672.5 Wim? K 


The inside heat transfer coefficient, h, can be obtained from a correlation for the 
Nusselt number for single phase coolant flow in a tube. This Nusselt number is obtained 


from a correlation (in this case) for the HEATEX element where: 


Nu, = 0.226(Re,)°*(Pr,)°* 
. 65 46 
Nu, = 0.226(559.89)°?(43.75)*° = 78.41 


and 


, _ Nucke _ (78.41)(0.3742) 
id «0,013259 





h; = 2213.07 Wim? K 


— 


The outside heat transfer coefficient 1s calculated from equation (5-12) where: 


l 








Mo 1 d, a 
Us meal drain ae 
The wall resistance, R,, is defined as: 
a, 
in <2 | 
Rm = 97k 
ia 0.015875 
0.013259 i 
Rk, = See — ee eee 


PA Weep yl Giskoy 
Therefore, 


l 


1 __001se7s 3 
672.49 ~ (0.013259)(2213.07) _ 6h 71*10 0.0608) 


h, = 1070.2 WJm? K 
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Figure 103 is a copy of the computer printout for this run. The small differences 
seen in the computer calculations versus these sample calculations are the result of 


roundoff error. 
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Data set number 


ame 


RESULTS FOR TUBE ia 


Mass flow rate 
Inlet temperature 
Saturation temp 


BEET 


temp 


Cais 


Log. Mean temp 
Heat flux 


Resistance of Tube Metal 


Prandtl number 
Reynolds number 
Reynolds number H&B § 
Inside h.t.c. 
Inside NUSSULT NO. 
GUERRA tc 


TUGe 
# 
} 


Figure 103. 


FM 
(2%) 
foes 


(Deg C) 


ache 


Uo) 


VEG 
(m/s) 
@.20 


ha Ons 


OraQNnNnowipreryaAe & oN 


DET 


(CK) 


.980E-@2 
.445E+00 
.693E+01 
.307E+O1 
.O48E+O1 
.@48E+04 
.913E-@4 
,559e+0) 
.6@8E+@2 
SE ye 
.216E+95 
Sse 57 6) 
~(20E +02 


Uo 


6.720E+@2 


HO 
CW7 moe eke) 
1.@68E+03 


Computer printout of results for run I data set SM06. 
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HNUS 


1. 2SJenGe 


APPENDIX F. UNCERTAINTY ANALYSIS FOR U5 


Run number | of data set SM06 was also chosen for the uncertainty analysis. The 
analysis was only conducted on U,. The measured and calculated parameters found in 
the sample calculation (Appendix E) were therefore used in this section. In order to as- 
sess the effects of coolant velocity on the uncertainty in the measurement of U, , an 
uncertainty analysis was also performed on run number 11 of data set SM06. This run 
had a coolant velocity of 1.2 m/s as opposed to 0.2 m/s in run number 1. 

The uncertainty analysis utilizes the procedure suggested by Kline and McClintock 
(Ref. 54]. This procedure states that if 


then the uncertainty in R, namely, 6 is given by the equation: 


2 2 ~ 2 ~ 2 0.5 
R= | ($4) +( OR ix.) Ser eancst oes + (22 ox,) 
Ox, Ox Ox, 


where x, is the measured variable and 6x, is the uncertainty in the measured variable. 





The uncertainty in the temperature measurements is due to uncertainty in the volt- 
age measured by the thermocouples. For the thermocouples used corresponds to 0.1 
°C. The uncertainty in in this thesis, this value was assumed to be 4 microvolts or 0.1 
C. The uncertainty in AT, T,,, and coolant inlet and outlet temperatures, 7,, and T7,, , 


can be calculated as follows: 


 T(10) + 71) 
com, 2 


Taking the partial derivatives we obtain: 





likewise, 





al 


The uncertainty in the measurement of 7,, 1s thus: 


ST OTe sr)’ OTs sr) 
coo} \ S710) "\ 87d) 


Based on the thermocouple uncertainty, 6J = 0.1 C, we obtain: 





5T >) = 0.0778 C 


The coolant inlet temperature is calculated from a single thermocouple, number 6 for 


tube A. The uncertainty in this measurement 1s calculated as: 








aia § 
6T(6) 
and, 
Bae or sT 2 $0.5 
oem || | IG) 
On 
(6T,=0.1C 


Finally, the uncertainty in AT 1s calculated from: 
SAT = [(5T 9)? + (5T ei) | 
and, 
OAT = 0.1347 C 


Coolant velocity 1s calculated from: 


y tt 
: pA; 


In this run, m: was calculated to be 0.02979 kg/s. The uncertainty in the reading of the 
flowmeter is defined as the scale interpolation. This value corresponds to one-half the 


value of the smallest marked increment. In this case, the scale interpolation factor, 6,, 


2 


is 0.5% which corresponds to 0.00101 kg/s resolution. In addition, a time-wise jitter, 


6,, of 0.5% was noted during the runs. Therefore, the uncertainty in m is given by: 


Om _ Om s \24 6m s \2\> 
m noo mi oY 


Or 


om 


— = 0.0479 
m 


The uncertainty in the cross-sectional area of the tubes was estimated from: 
5A, 6d, |? 
ica E cee 


where 6d, is given as 0.1 mm based on tolerances supplied from the manufacturer. 


Therefore, 


5A; 
—~ = 0.0107 = 1.07 % 


i 


The uncertainty in the coolant velocity is now given by: 


OV, | 6m \2 OA, i 


hn 





Or 


SV, 
= 0.049 = 4.9 % 


“ 
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The uncertainty in the calculation of coolant Reynolds number 1s given by: 
bRe bV. Sah i ie 
7m rae oa) 
Kee 4. qi 


or substituting appropriate numbers we find 





6Re, 
—— = 0.0497 = 4.97 % 
Re, 


The uncertainty in the calculation of the heat transferred to the coolant is related 
to the uncertainty in the measurements of coolant mass flow rate and coolant temper- 


ature rise where: 


bq es AONE 
rane wklovee 


or 


- 
os = 0.0489 = 4.89 % 


The uncertainty in the heat flux is then: 


Og 6g \2 |, 9A ll 
= — )°+(— 
fe, [ ; Pie vy ) 


where the uncertainty in A, is given by: 
5A bd a 
A, d, L 


Here, the uncertainty in the length of the tube, dL , is assumed to be 5 mm while 6d, 1s 





assumed to be 0.1 mm (again from given tolerances). Substitution of the appropriate 


numbers vields: 


a 
£40 _ 9.9496 
q 


Oo 
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Finally, since U, is given by: 


q o 


Vo = TTD 





the uncertainty in the LMTD must be calculated. The LMTD is given by: 


AT 
as a (ay 
ln —_—_ 
eae _ lee 


In order to calculate the uncertainty in the measurement of the LMTD, the following 


LMTD= 


equations were used: 


@LMTD _ 
OAT ee = ie 
In | as ae 
sat C,0 
CLMTD = NJ << (iS a Sat T.;)] 
ier © Lim 2 
of = aes fs al ie oye Sat Te) 
Sal C,0 
@LMTD __ AT 
ee Pe T 2 
“ In) | (Trt — Teo) 
oe = ea 
@LMTD _ AT 


OT, 5 Tea i ea : 
ola ap (le (a) 


The uncertainty in the LMTD 1s therefore given by: 


! 
OoLMTD 2 oLMID 2 oLMID 2 oLMID A |e 
ITD =| ( = 24 ( Coe on ST, C= 5T,, 


For the sainple run with a coolant velocity of 0.2 m/s, the uncertainty in the measure- 
ment of the LMTD is: 


ag 
SLMTD = [[(2.329)(0.1347)}° + [( —1.0157)(0.0778)]? + [( —2.9097)(0.0778)]° + [(1.894)(0.1)]"]” 
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or 
OLMTD = 0.446 K 
and 


6LMTD _ 0.446 
LMTD 30.48 





= 0.0146 = 1.46 % 


The uncertainty in the overall heat transfer coefficient, U,, is then given by: 


: ” Ble 
OU, _| f ©To + ( 5LMTD ) 2 
U, Gia LMTD 


Or 
6U , ant 
sii, = [(0.0496)° + (0.0146)*]2 
and 
OU, 7 
Pls = 02.052 =5.2 % 


Table 7 summarizes the results of the uncertainty analysis for the smooth copper tubes 
with the HEATEX insert for coolant velocities of 0.2 and 1.2 m/s respectively. 

The uncertainty analysis neglects any uncertainty in the physical properties of the 
coolant or the R-113. For the most part, this is probably a reasonable assumption. The 
exception to this may be the viscosity due to its sensitivity to temperature. Hence, the 
uncertainty reported here is an underestimate. 

It is curious that the relative magnitudes of the uncertainty for the low and high 
coolant velocities are essentially the same. Two different mechanisms dominate in the 
calculated uncertainty at low and high coolant veiocities. In terms of the mass flow rate, 
the time-wise jitter and scale interpolation terms are the same for both coolant velocities. 
However, the actual mass flow rate is less at the lower flow rates. Hence, the uncertainty 
in the mass flow rate calculation is more significant at the lower coolant velocities. At 
higher flow rates the uncertainty in the calculation of the LMTD becomes the dominant 
term. This is because the coolant temperature rise is smaller at high coolant velocities 


but the uncertainty in the temperature measurements remains the same. The mass flow 
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Table 8. SUMMARY OR RESULTS FOR SAMPLE UNCERTAINTY ANALY- 
SIS. 


0.049 0.013 
i 

SLMID. 0.015 0.054 
~LMTD~ 
U, 


rate goes up so that the uncertainty in the coollant mass flow rate calculation decreases. 









3 5 






s| 2/3 | 










Gn 
rn 








Thus, at the higher coolant velocities, the uncertainty in the calculation of the LMTD 


is the dominant term. 
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APPENDIX G. PROGRAM LISTING 
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10@@ ! FILE: DRPCON7 

1Q9@@S5 ! PURPOSE: This program collects and processes condensation data for the t 
ube bundle. 

1015 ! CREATED: NOVEMBER 2, 1988 

bee?! UPDATE1: APRIL 25 1981 (S. MEMORY) 

Meg! : UPDATEZ: May 17 1391 (S. MEMORY) 

m@ee! UPDATES: JULY 26, 1991 (R.W. MAZZONE) 

1023 ! UPDATE4: 13 AUGUST 1991(R.W. MAZZONE ) 

1024 ! UPDATES: 1@ OCTOBER 1991 (R.W. MAZZONE) 

1@26 ! CHANGE HOD=4 AND USING CALCULATED CI FROM WILSON 
1Q3@ REEP 

1035 PRINTER IS 1 

1Q4@ PRuUNiewouNo “4x, ° "SELECT OPTION"”” 

1@4S PRINT USING “BX,""@ TAKING DATA OR REPROCESING PREVIOUS DATA" "" 
1@5@ PRINT USING “6X,°°! PLOTTING H VS DELTA-T""" 

1055 ProiewehNia GA) "2 PLOTTING HRAT VS N*** 

106@ CRiNiee@stNe “6X,° 3 PLOTTING WILSON""" 

1@65 PRINT USING "B6X,""4 PURGE FILES"*"" 

1@7@ BeniieetotNe 6X%,'"S XYREAD”"” 

1Q@75 PRINT USING “6X,""6B NUSSELT ESTIMATE” "" 

1Qa8@ PRINTER IS 7@1 

1985 INPUT Icall 

199@ IF Icall=@ THEN CALL Main 

1995 IF Icall=1 THEN CALL Plot2 

11a@ IF Icall=2 THEN CALL Plot] 

11@5 TF Icall=3 THEN CALL Plot3 

111@ IF Icall=4 THEN CALL Purge 

aS IF Icall=S THEN CALL Xyread 

tues IF Icall=6 THEN CALL Nusselt 

at25 END 

rsd SUB Main 

55 COM eC cvaet 7 > 

114@ COM /Baed/ Ifa 

nas COM /Nus/ Tin,Tsat,Qdpl,Hnus,Kf ,Rhof,Hfg,Muf ,Do,Itube 
1159 COM /Will/ Doa(4),D1ia(4),Kma(4),Tact ,Droot(4) 

mSS COM /Wil2/ Delta,Isat ,Nsets,Hod,Cia(3,3),Alpaa(3) 
116@ Oiveemr ely ots ),1(17),Ho(4),Qde(4),U0(4),Pc( 4) 
moe, DATA 1.0,1.0,1.0,1.8@ 

W216 Pewsey amore ee, 5.1 %2.5.172 

e217 Pain 8.226 ,0.226,0.226,0.226 

HZ 18 DATA @,063,@.063,0.@63,@.@63 

1229 READ Ciat+#) 

225 Bere OO. 1O0S6O91 25727 .94469,-7674545. 829578025595 .381 
123@ Olicmetcovesoecd.6,.c7b>o5E11,—-2.66192E13,35.94078E14 
255 READ CC) 

1240 DATA @.@15875,@.01400@ ,0.@15798,0.0159,0.0 

1245 DATA @.@13259,@.@1@160,0.@12776,0.@1346,0.@ 
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125@ DATA 386.0.42.975 422975, 2023-040 

E25) DATA @.@,@.014,0.@,.@141,@.0 

Less DATA @.0@05588, 3 

1750 READ Doal*#),Dial*),Kma(*#),Droot(*#),Delta,Hod |! Hod=H/Di1 

Pes L=1.2192 ! Condensing length 

127@ Jset=@ 

i271 Okaccpt=1 

LZ BEER 

128@ INPUT "ENTER MONTH, DATE AND TIME (MM: 0D:HH:MM:SS)”" ,Otg$ 

i2oc OUTPUT 709; 1B :Bies 

eae BEEP 

ies INPUT "SELECT OPTION (@=DAQ, 1=FILE)",Im 

13a@ Thard=1 

13@5 BEEP 

1312 INPUT "WANT A HARDCOPY PRINTOUT (1=DEF=YES,@=N)",Ihard 

Pas BEEP 

1320 INPUT “SELECT (@=#R-114,1=STEAM,2="R-113,5=EG)" ,1ft 

1325 Tin=1 

13350 Isat=2 

EAS REEP 

1340 INPUT "SELECT SAT TEMP MODE (@=LIQ, 1=VAP,2=(LIQ+VAP )/2=DEF )" Isat 

PaaS IF IThard=1 THEN PRINTER IS 7@1 

135@ ~ IF Im=@ THEN 

ei si= Beet 

136@ INPUT "GIVE A NAME FOR THE NEW OATA FILE", Files 

i255 CREATE BDAT Files, 20 

137@ ASSIGN @File TO File$ 

275 BEEP 

138@ INPUT "ENTER TUBE CODE (@=SMOOTH CU,1=FINNED CU,2=KORODENSE,3=TIT FIN 
yi ube 

lset sielei> 

1382 INPUT "@=NO INSERT, 1=TWISTED TAPE, 2=HITRANI, 3=HITRAN2",Insert 

1385 BEEF 

1392 INPUT "ENTER EG CONGENTRATION <WT FPERCENT ) = bean: 

aia ENTER 7@9;,0tgGs 

1402 OUTPUT 6File;:Dtg$ 

14@5 OUTPUT @File;Itube,Egrat,O0d1,0d2,0d3,0d4,0d5 

141@ Tact=1@ 

415 el S32 

142@ INPUT “SELECT (@=TOP, 1=SECOND,...., 1O-BUNDEE-DEFAVET) =~ lace 

ae oo PRINT 

1435@ PRINT USING “JOX Fille Neate ( IZA Gell es 

Tass PRINT 

144@ alias 

1445 BEEF 

1452 INPUT “ENTER NAME OF EXISTING FIRE @ewes 

1455 BEEP 

146@ INPUT "ENTER NUMBER OF DATA SETS STORED" ,Nsets 

1465 Twil=1 

147@ REEF 


mele) 


ms INPUT "WANT TO CALL WILSON? (1=DEFAULT=YES,@=NO)",Iwil 

148@ IF Iwil=1 THEN 

1485 BEE 

149@ INPUT "WHICH TUBE (@=TOP,1=SECOND,..,1@=BUNODLE)",Iact 

ibeshal INPUT “@=NO INSERT, 1L=TWISTED TAPE, 2=HITRANI, S=HITRAN2", Insert 
1493 CALL Wilson 

1494 IF Tact=1@ THEN 

1495 FOR Icount=@ TO 3 

1496 Cia(Insert,Icount )=Cia(Insert ,@) 

1497 Alpaa(Icount )=Alpaa(@) 

1498 NEXT Icount 

1499 NOE 

150@ IF Ihard=1 THEN PRINTER IS 701 

15@5 PRINT 

151@ PRINT USING “10X," "FILE NAME: o4 2A eles 

1a15 PRINT 

1520 ieedicet 23 THEN 

1S2 1 FOR IT1=@ TO &S 

Le22 PaiNimewounG  lOXe | RESULTS FOR TUBE NO."" sD-; 1 tl 

1526 PRINT 

1525 PR sUotNG 1Oke INSIDE COEFFICIENTS: ~~ ,00.50";Cialingert ,11 
) 

153@ PRU imUSTNGmet@x , ~ ALPHAS: Oe M2. SUE eee Alo 
aa({I1) 

Leo NEXT Ii 

P45 EiSE 

1540 SR eect Noman, INSIDE COEFFICIENT SEOR TUBE “",0," 9205). 00.3 
O"sTact+1,Cia(Insert,Iact) 

1Ses Geis tNG “10%, "ALPHA FOR TUBE ~" 852": ceerD. o 
OD"; lIact+1,Alpaa(lact ) 

1550 END IF 

1655 ERSE 

1560 BEEP 

156) PNB] -@O=79P, I=SECGOND, ..... 1@=BUNDLE" , Tact 

1562 INPUT "@=SMOQOTH CU, 1=FIN CU, 2=KOROD, 3=FIN TITANIUM" ,Itube 
EoaS INPUT "@=NQ INSERT, 1=TWISTED TAPE, 2=HITRANI, 3=HITRAN2Z",Insert 
1565 INPUT “ENTER CI VALUE ”,Ciat( Insert ,11) 

157@ END IF 

ses) ASSIGN @File TO Files 

1580 ENTER @File;Dtg$ 

Sees | ENTER @File;Itube,Egrat ,Od!1,0d2,0d3,0d4,0d5 

1586 INPUT "@=SMOOTH CU, 1L=CUNTI FIN, 2=KOQRQDENSE, 3=FIN TITAN “,I tube 

159@ END IF 

Jags Tout=1 

16QQ BEEP 

16@5 INPUT "WANT TO CREATE AN QUTPUT FILE? (1=DEF=YES,@=N@)", Tout 

161@ IF Tout=l1 THEN 

Ie le sell 

162@ iE sence aeNAMeE FOR GCUTPURBEILE 9 (Fouts 

Heo CREATE SDAT Fout$s,5 


Lisl 


1630 
1635 
1640 
1645 
1646 
165@ 
less 
1660 
1665 
1672@ 
US fal 
1672 
ee oS 
1674 
1675 
1676 
1677 
is 93 
168Q 
1682 
1683 
1685 
169@ 
1691 
/s)? 
1694 
P6455 
1696 
1697 
1698 
1699 
17@@ 
1701 
17@2 
1703 
1704 
17@5 
17@6 
177 
17@8 
ied 
t7i 
teZ 
re les 
ie 
Lys 
1716 
Vy 
1718 
1719 


ASSIGN ®Fout TO Fouts 
ENO IF 


Do=Doa(Itube) 
Dol=Doa(I tube) 
D1=Dia(Itube) 
Km=Kma( I tube ) 


Ax=PI#90i*2/4 ! Cross-sectional area 


Ao=P1*#Do#Ll 
! = =Rm=Do*#LOG(Do/Di )/(2*Km) 
IF Itube=1 THEN 
Dol=Droot(Itube) 
END IF 
IF Itube=3 THEN 
Dol=Droot(I tube) 
END IF 
Ao=PT*#Do#Ll 
Rm=LO0G(Dol1/01 )/(2*L*Km ) 
Al=PT*#Dole#L 


IF Im=@ THEN 
PRUNTER@0Seet 
BEEP 


INPUT "WHAT IS REQUIRED NOMINAL COOLANT VELOCITY IN m/s (MAX. 


",Vegst 
OUTPUT 709;"AR AFG ALG VRS 
FOR IT=6 TO Q 
OL Up ylahsie ies a Ss) 
VUsum=@ 7 
FOR J=1 TQ 5 
ENTER 7@9;E 
Usum=VUsumteE 
NEXT J 
Emf( 1 )=Vsum/S 
TC IJ=FNTvsv(Emf¢(1)) 
NEXT I 
Tavgl=(T(6)4+T(7)4+T(8)4T(9) 


)/4 


IF Tavgi<@ OR Tavgol?>24 THEN 


PRINT "COOLANT TEMP. NOT IN RANGE @ TO 24 


GOTO 1691 
END IF 
Rhoegi=FNRhoeg( Tavgi,Egrat 
Mdotl=Rhoegl *Ax*Vegst 
Pca@=(Mdot1+.00414)/.Q002@06 
Poa24=(Mdot1+.@@329)/.0Q@22 
Pcob@=(Mdot1+.Q00615)/.00204 
Peb24=(Mdoti+.@0138)/.0@21 
Pec@=(Mdot1+.@@36 )/.90208 
Pec24=(Mdoti+.@@3@1)/.@@22 
Ped@=(Mdot1+.@@122 )/.@@28 


) 


5 


6 


2 


Ca 


\ 


4 


™ 


7 Ny 
2) 
re? 2 
ihegig' 
1724 
yes 
1726 
ie? 
1728 
7 Pa 
173@ 
b7si 
Pe 
pee 
Pras 
1736 
laa | 
i738 
Iga: 
1742 
1741 
1742 
iy ae 
1744 
1 Qe 
1746 
Pea 7 
1748 
1 evs 
1752 
1 7S Lh 
1752 
755 
1754 
17 as 
1756 
Liag 
i752 
el2 

P7259 
1762 
1761 
1762 
7 bo 
1764 
1765 
1766 
rye? 
1768 
1769 


Pod24=(Mdoti-.@048)/.@0239 

Po(@)=PcaQ@+( Tavgl*(Pca2z4-PcaQ) )/24 
Pe( 1 )=PchQ@+( Tavgl*(Pcb24-Pcb@) )/24 
Pe (2 )=PecOt+( Tavgl*#(Pec24-Pcc®) )/24 
Po(3)=PcdQ@+( Tavol*(Pced24-PcdQ@) )/24 


PRINT “SET FLOWMETER READINGS CORRESPONDING TO: " 
PRINT SUGINGe ax, ~§4 OF METER A = **,DDD.D";Pc(Q) 
PRNiwotNGe: ba, 4% OF METER B = “* DDO.D*:Pc(1? 
vivo lNGe bAy0e UF METER © = DD DeaR"sPc< 2) 
pemeesiNombn, 64°0F METER 0 = ~",0DD.D°;Pets) 
PRINT 

PRINT "HIT CONTINUE WHEN READY" 

wall)s)s 


Repeat: 
PRINTER IS 7@1 


Dt 1ld=47.5 | Desired temperature of liquid 
Ido=l 

ON KEY @,15 RECOVER 1734 

PRINT USING "4X,""SELECT OPTION """ 

PRINT USING "6X,""@=TAKE DATA """ 

PRINT USING "6X,""1=SET Tsat (DEFAULT) """ 
PRINT USING “4X,""NOTE: KEY @ = ESCAPE""” 
Ido = desired option 

BEEP 

INPUT Ido 


Set default value for input 
IF Ido>l THEN Ido=l 

Take data option 

IF Ido=@ THEN 1797 


Loop to check saturation temperature 

IF Ido=1 THEN 
fi eNien Ofoihedeisay (DEPAULT=47, 5 5@—— (R-1is)" Otlea 
Peete PESTREG Tisat (DEFAULT=Z2.2 C = R=-114)",0tld 


Nn= 


j 


Nes=Nn MOD 15 


Nn= 


Nnt+l 


IF Nrs=l THEN 


PRINT USING 
Dells 


Del4 


END IF 
Read thermocouple voltages 
DUR 2S, AR APO ALI? VRS” 
Nend=17 
FOR T=@ TO 5 

QIU PIU Fauhsts Sissy) Syaia 


Vsum= 
J=1 TO 5 


FOR 


wax. "oDTsat ileleeit Tvat 


Tvab 


Q 


ENTER 7@9:;E 
Vsum=Vsumt+E 


NEXT 


J 


63 


Dell 


Bp 


177@ Emf(T)=VYVsum/S 


ee aah TCT =FNTvsv(Emf(l)) 
Due NEXT I 

Lie 3 FOR I=6 TO Nend 

1774 OUTPUTS OS4 eon sr. 
Vio Vsum=@ 

1726 FOR J=1 TQ 20 

LY ENTER 7@9;E 

1778 Vsum=VsumtE 

yy is NEXT J 

178@ Emf( I )=Vsum/29 

bg siah TCI J=FNTvsv(Emf¢(l)) 
17e2 NEXT I 

1785 -! Compute average temperature of liquid 
1784 Tlig=(1 Ca )4i( 49972 

1785 Tvap=(T(@)4T¢C 194702) )/3 
1786 Tvat=(T(@)4+T(1))/2 

1787 Tligl=T¢(5) 

1788 Dell=((7(19)4+1(11))720- 7G» 
i7eo DelZ=¢(<T< 12 410135007 ee: 
1730 Dels=((7( 14947 (lage sie 
ia Del4=((7T(16)4T(617))/2)-TC(9) 
ieee PRINT USING "4X,8(MD0.D0,3x%)";Dtid,7(35),1@2)jivat Bell , Deleware 
14 

lyase WAIT 2 

1794 SUTOMI7s5 

vals ENG wie 

ieee | TAKE BATA IF Im=@ LOQP 

lity = We PRINTER sean 

ste: OUTPUT 709; "AR AF@ AL17 VRS" 
17939 Nend=17 | INCREASE TO 20 IF FIVE TUBRESMINSeUNiese 
18@@ FOR I=@ T0 5 

18901 OUTPUT =7 2S omc 

18@2 Usum=@ 

1803 FHR 3-1) 10-5 

18@4 ENTER 795; 

18@5 Vsum=Vsumt+e 

18@6 NEXT J 

18Q7 Emf( I )=VUsum/5 

1808 NEXT I 

1869 . FOR I=6 TO Nend 

181@ OUTPUT 7@9;:"AS SA" 

1811 Vsum=@ 

Leute FOR J=1 TO 20 

shite ENTeRe? Gane 

1814 Usum=VsumtE 

1815 NEXT J 

1816 Emf( I )=Vsum/2@ 

182@ NEXT I 


Lec ELSE 


1 84 


1830 
1245 
1840 
1845 


185@ 


1855 
186a 
1865 
187@ 
1875 
18a@ 
1881 
ee5 
1890 
1335 
1896 
190@ 
1905 
19@6 
1907 
19@8 
1325 
Usa 
193@ 
25 
1940 
1945 
195@ 
is ies 
196@ 
1965 
no 7@ 
P75 
198@ 
28 1 
se eden 
1990 
1 ape p 
L3s2 


1994 


ios 
1386 
2202 
2@@5 
2210 
2@15 
2A20 
2@25 
2030 
2@35 


ENTER @File;Pc(*),Emf(*#),Tp(*) 
END IF 


DATA ANALYSIS 


Nend=17 
FOR I=@ TQ Nend 
TCID=FNTvsv(Emf(l)) 
NEXT I 
Nvap=(1¢@)+7< 1)4762))7/3 
Mlve=(16sI41(4))72 
Tligi=T(5) 
IF Isat=@ THEN 
Tsat=Tliql 
END IF 
IF Isat=1 THEN 
Tsat=Tvap 
END IF 
IF Isat=2 THEN 
Tsat=(Tvapt+Tlig)/2 
END IF 
Jset=Jsett+l 
IF Qkaccpt=@ THEN Jset=Jset-l 
PRINT 
PRINT USING "10X%,""Data set number = "",DD";Jset 
PRINT 
Theg=@ 
Tend=3 
IF Tact<1@ THEN 
Ibeg=lact 
Tend=lIact 
END IP 
FOR I=Ibeg TO Iend 
Grad=FNGrad(Emf(I+6)) 
Mout=(7(2*14+1@)4+7(2*1+11))/2 
Delt=Tout-TCI+6) 
Tavg=T(I+6)+Delt+*.5 
IF I=Iend THEN 
Delt=T(1@)-TC(9) 
Tavg=1(1@)-T(9)/2 
ENG IF 
Rhoeg=FNRhoeg( Tavg,Egrat ) 
Nueg=FNNueg( Tavg,Egrat ) 
Mueg=Nueg*Rhoeg 
Cpeg=FNCpeg(Tavg,Egrat ) 
Keg=FNKeg(Tavg,Egrat ) 
Preg=Cpeg*Mueg/Keg 
Mdat=FNemeal<(1,1( 1+6),Rek Ip 
Veg=Mdot/(Rhoeg*A~ ) 
Reeg=Veg*#Di/Nueg 
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20402 
2245 
2250 
2@55 
2260 
2@65 
2070 
2275 
2280 
2081 
2285 
2886 
2@88 
2@89 
2890 
2091 
2032 
2094 
2@935 
2@96 
2297 
2298 
2@9q 
2180 
21@1 
ZUO2 
2104 
2105 
21@6 
Fea sh MMs 
Fad hes 
2120 
ZS 
Zan 
Zio) 
eee 
2134 
Zsa 
2 he5 
Felines a 


21a | 


214@ 
2145 
2146! 
2147 
2148 
2uS0 
Fo cai 
Zia! 


Res=4*Mdot/(PI*Mueg*#(Di-4*Delta)) 
Qdot=Mdot*#Cpeg*Delt 
Qdp( I )=Qdot/fAo 
IF I=@ OR I=Iact THEN 
Qdpl=Qde(1) 
CALL Nusselt 
END IF 
Lmtd=Delt/LOG((Tsat-T(1+6))/(Tsat-T¢(1+6 )-Delt)) 
Uol I =Qdp(I)/Lmtd 
Cfeg=l 
IF Reeg<4@0@ THEN 
IF Insert=@ THEN 
Bbh1=.@668+*(Di/L )*Reeg*Preg 
Bb2=1+.@4*((Di/L)*Reeg*Preg)*.6666 
Nueg=3.66+(Bb1/Bb2 ) 
ENDS Te 
IF Insert=1 THEN 
Nueg=1.86+*( Preg*Reeg )* .333333*(Di/L)*.33343*Cfeg 
Nueo=Cia(Insert,1)#(1+5.484E-3*Preo”.7*(Res/Hod? 1°25 ae 
END IF 
IF Insert=2 THEN 
Nueg=Cia(Insert,1I)*(Reeg*.65 )*#Preg” .46 
END 
IF Insert=3 THEN 
Nueg=Cia(Insert,1I)*(Reeg”.76 )*Preg*.46 
END IF 
ise 
BEEP 
PRINT USING “1@X,""“TURBULENT CORRELATION®*" 
Nueg=.@27*Reeg*.8*Preg®.3333*Cfeg 
a [8 108 
Hi=Nueg*Keg/Di 
Ho(T)=1/¢(1/Uo0( 1 )-Do/(Di*Hi )-Rm*Al ) 
IF I=@ OR [T=lIact THEN 


PRINT 

PRINT 

PRINT USING "10x, “"“RESUETS FOR TUBE. ste ret 

PRINT 

PRINT USING "1@X,""Mass flow rate = “",MZ.3DE° 7 Maan 
PRINT USING "1@X,"“"Inside Tube Dia. (m.) = "" MZ. 3DE";Dia(Itube 
PRINT USING “10X,""*18@ DEG OVER Oaa. GHOD) = %,MZ.3DE"; hee 
PRINT USING "10X,""Inlet temperature =n yMZ.30E" ; ieee 
PRINT USING “"1@X,""Saturation temp (Deg C) = "",MZ.2DE";Tsat 
PRINT USING "10X,""DELTA Tape Thickness = ""(MZ.3DE";ODelta 
PRINT USING "1@X,""QDELT temp Oif. = "" .MZ.30E”; Dem 
PRINT USING "10X,""Log. Mean temp ODif. = "! .MZ.4DE" ; Lowe 
PRINT USING “1@X," “Heat ein = = “"MZ.3DE" : Odea) 
PRINT USING "10X,° “Conductivity e.ce = 9° ,MZ.30E": Keg 
PRINT USING “10X,""Conductivity Tube Metal = "",MZ.3DE";Km 
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POBET: Rm 


SUE Pee 


. SDE" ;Reeg 
»~ADE";:Res 


Sole Hi 


. SDE"; Nueg 


Stee 7 UStie) 


Siew. kh 


»oWE > Riven 


ANS Sahie 


Powe > Mut 


Uo 


Civic 2 


sO ait 


2153 PRINT USING "1@X," "Resistance of Tube Metal= "",MZ 

2beaS PRINT USING “1@X,""Prandti number Ser rae 

21690 PRINT USING "“1@X,""“Reynolds number = "" MZ 

2151 PRINT USING "1@X,""Reynolds number H&B S Ses My, 

Ss PRINTOUSING “~10X,° “Inside h.t.c. = aie 

AUIS s) PRINT USING "10X,""Inside NUSSULT NQ. SP Miz 

ZASG PRINT USING "1@X,""QVERALL H.t.c. (Udo) = elles. 

Zn © PRINT 

eigi PRINT USING “10x, ° “KF SP ee 

MN GPa PRINTSUSING» “10K,—° RHOF ee OAM 

Zig S| Pei aiNae 10%.) HFG = ee Fag 

2174! PRINT USING "10X,""MUF = "MZ 

2175 PRINT 

Zid PRIN@eUoING 19x." TUBE FM VEG GET 
HO HNUS"" 

2180 PRUNT USING “10X,°" # Ch) (m/s) (K) 

oo 

ees ee musiNon  FOx, of 4X%-50.0D,354 ,2.00,4X.00.350,2XxX.M2, SDE, 3X Mz. SDE 

map? . SWE, SX"; ([+1,Pc(1),Ven,Delt,Uol(lI),Hol(I),Hnus 

2190 diss 

Zig SicotNon Oc our 4X. 50,00 ,4%,.2.D0,4% DRS. 2% (MZ .50E,3xX (MZ. SDE 

mere it! Pc( 1) Veg Delt ,Uo( I), Ho(l) 

2200 END IF 

22¢@5 NEXT I 

2210 IF Im=@ THEN 

2215 Qkaccpt=1 

2220 eelele 

2225 vier nOeAeGePl THIS@SeET (1=DEFAULT=YES, 9=N0Q907" ,Okaccot 

2235@ IF Qkaccpt=1 THEN QUTPUT @File:Pc(#),Emf(*#),Tp(#) 

2255 END IF 

2240 IF (Qkaccpt=1 QR Im=1) AND Tout=1 THEN 

2245 FOR I=0 TO 2 

2250 OUTPUT Seeut wiecl ),Oda< 1) 

2255 NEXT I 

2269 eR Ol SS 

22bs IF Im=@ THEN 

2210 Okrpt=l1 

2275 REEP 

2280 INPUT "WILL THERE BE ANOTHER RUN (1=YES=DEFAULT , @=NQ )" 

2285 IF OQkrpt=1 THEN 1683 

2290 el Sis 

2295 IF Jset<Nsets THEN 1683 

2380 ENO IF 

2305 ASSIGN ®File TO * 

Zon0 IF ITout=1 THEN ASSIGN @6File TQ * 

@o1> SUBEND 

232@ DEF FNGrad(T ) 

2429 Grad=-3.8778S7E-S-244.71428S57E-8*T 

2330 RETURN Grad 

23355 FNEND 


iow 


2340 
2345 
2350 
2656 
2360 
VASES 
2370 
2513 
fae 

385 
2390 
2595 
242 
2405 
2412@ 
2415 
2420 
2425 
2440 
2435 


244Q 


2445 
2450 
2455 
2462 
2465 
247@ 
2475 
2480 
2485 
249@ 
2495 
2582 
2005 
251@ 
Zale 
252@ 
205 
253 
2350 


2542 . 


2545 
2558 
agen os 
2564 
2565 


2570 | 


Z2o72 
258 
2565 


DEF FNKcu(T) 

QFHC COPPER 25@ TQ 3@@ K 
Tk=T+2735.15 
K=454-.112*Tk 
RETURN K 

FNEND 

DEF FNNueg(Tc,Egr ) 
RANGE OF VALIDITY: 
Th=tetz7 4.15 


720 [0 20 Hebe 


Nul=7,1196507E-3-Tk#*(7,.4863547E-S-Tk*#(2.6294945E-7-1k #35. 085332 9S eee 
Nu2=4.9257658E-3-Th#(4.9213912E-Sa)k*(1.6457534E-7-(k*1 83553535 leeeee 
Nu3=8 .6586293E-5-Tk#*(8.88579@2E-S-1Tk#*(3,0495032E-7-][k *3.499999b— eee 


A2=(Nu3—-2*Nu2+Nul )/20@ 
Al=(Nu2-Nul-94@0#AZ )/1@ 
AQ@=Nui-42#A1-1764*A2 
Nu=AQ@+Egr*(AlL+Eor #A?2 ) 
RETURN Nu 

FNENO 

DEF FNCpeog( Ic. Ear? 
RANGE OF VALIDITY: 
Dr=Te+274,5 
Col=1.6701S50E+5+ eS 

Co2=1,4746125E+4+Tk*6 25 

Cos=9.5800500E+2+TK*#7.3 

A2Z=(Cpo5—-2 *€a2+Cn 17200 

Al=(Cp2-Cp1-900+*AZ2 )/1@ 

AQ=Cp1-4@*A1-168Q#A2 

Co= min tere Sear, 

RETURN Cp 

FNEND 

DEF FNRhoeg(T,Egr ) 
Rol=1,@6@7095E 45-7 *( 5.705 12836-14[*4 O857 1 ese oo 
Ro2=1.0748272E+3—-1T#( 4. 42661 9S5E—1 41 #4 69347 @bee»D 

Ro3=1 ,G@885S54E+4—1 *( S575 5S6S5S9E—14) 46.12 tag oes 
A2=(Ro3-2#Ro2+Ral )/2a0 

Al=(Ra2-Rol-9@Q0#A2 )/1@ 

AQ=Ro1l-40#A1-160@*A2 

Ro=AQ+Egr*(Al+Egr #A2 ) 

RETURN Ro 

FNEND 

DEP VENPred 1 sean 
Pr=FNCopeg(T,Eor )*FNNueg(T,Eor )*FNRoeg(T,Eor )/FNKeg(T,Ear? 
RETURN Pr 

FNEND 

DEF FNReGak te bene 
RANGE QF VALIDITY: 
Th=le+27 ses 
Ki=2.2824708E-1+1k «(5S S9892e66-44 104s) 57 ee coe ee 
K2=2.5984661GE-it+DRe( 2° 5976S 7te 44 eee 7 


WO TOr 2espeGre 


=20) 10-20 DEGaEe 


~KS=3.2138952E-1-Te#( 4200428576 —4—1 ie eee ee 
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2590 
oes 
2600 
26@5 
261@ 
Zeiss 
2620 
HAS) fig 
2630 
2639 
2640 
2645 
2650 
2655 
2660 
2665 
2670 
2500 
2682 
2ee5 
2690 
2645 
2100 
2/@5 
2/1@ 
eels 
21720 
2725 
2730 
2235 
274@ 
2746 
2750 
Z1as 
276@ 
ed fs) 
2779 
ati 
2850 
2955 
2860 


2865 


28702 
ZedS 
Z288Q 
2885 
2890 
2835 
2900 
2205 


A2Z2=(K3-2*K2+K1 1/200 
AL=(K2-K1-90G+*Az2 )/10 
AQ@=K1-40+*A1-1600*A2 
K=AQ+Egr*#(ALt+Egr*#Az ) 
RETURN K 
FNEND 
DEF FNTanh(X ) 
P=EXP(X ) 
Q=1/P 
Tanh=(P-Q)/(P+Q) 
RETURN Tanh 
FNEND 
DEF FNTvsv(V) 
COMea7zee. C(7) 
T=C(Q) 
FOR I=1 TO 7 

T=T+C( 1 )*¥U%] 
NE Xia 
RETURN T 
FNEND 
DEF FNBeta(T) 
Rop=FNRho( T+.1) 
Rom=FNRho(T-.1) 
Beta=-2/(Raptkom )*(Rop-Rom)/.2 
RETURN Beta 
FNENQO 
DEF FNPsat(Tc) 
@ TQ 8@ deg F CURVE FIT OF Psat 
Pt leo et aZ 
Pa=S,945575+TF#( .159552082+TF*#(1.4840963E-3+TF*9.6150671E-6)) 
Pgo=Pa-14.7 
IF Poo@ THEN 

Psat=Pq 
Sle 

Psat=P9#29.92/14.7 
ENO IF 
RETURN Psat 
FNENO 
SUB Xyread 
BEEP 
INP ENTER FILE NAME’ Files 
REEP 
Ne ieee ner NUMBER OF X,Y PAIRS’ ,N 
ASSIGN @File TO File$ 
FOR I=1 TQ N 

ENTER @File;X,Y 

Roy 
NEXT I 
SUBEND 
SUB Purge 


Pet -PSiG,-—in Ho 


189 


2912 
fe Ae 
Zac 
2025 
2930 
Fikira fogre 
2940 
2945 
2959 
2354 
2962 
2965 
2970 
facta: 
2976 
2507 
2980 
fais! 
2990 
22495 
3990 
3QQ5 
3@1@ 
S@15 
3929 
3025 
3Q5@ 
4035 
3@4Q 
3@45 
3@46 
3059 
S@o1 
3@52 
S0a3 
3@54 
S055 
3056 
3057 
3958 
evao 
3@6@ 
3062 
4063 
3@65 
3072 
3@75 
3Q8@ 
4985 
3@9@ 


Bee 

INPUT “ENTER FILE NAME TO BE DELETED" ,FileS 

PURGE Files 

GOTO 291@ 

SUBEND 

SUB Wilson 

COM /Will/ Doa(4),Dia(4),Kma(4), Tact ,Droot(4) 

COM /Wil2/7 Delta, leat Neets Hoceeracs 3) Aleaaca: 

DIM Emf(l7), 1117). hae vac 10g eee Siena) 

BEER 

INPUT "PLEASE RE-ENTER NAME OF FILE" ,File$ 

ASSIGN @File TO Files 

INPUT "@=SMOOTH COPPER TUBE, 1=CUNI FIN, 2=KQRODENSE, 3=FIN TITAN", Itube 
Beer 

INPUT "“@=NO INSERT, 1 
BEEF 

INPUT “GIVE 4 NAME EGR ea - lees 
CREATE BDAT Xy$,20 

ASSIGN @Xy TQ Xy$ 

EEte2 a2 

Do=Doal( Ttube?) 

Di=Dia(Itube) 

Km=Kma( Itube ) 

Ax=PI*Di*2/4 | Cross-sectional area 
Ao=PI*DoetLl 

Rm=Do*#LQG(Do/Di )/(2*Km) 


TWISTED TAPE INSERT, 2=HITRANI, 3=HITRAN2",Insert 


Initial values 
Tf=Tsat 
Alpa=.655 

IF Insert=@ THEN 
Ci=1.9 

END If 

IF Insert=1 THEN 
Ci=Se 172 

ENO Se 

IF Insert=2 THEN 
Cir= 2226 

END IF 

IF Insert=3 THEN 
Ci=.@63 

EN le 

6=9.81 

Theg=@ 

Tend=@ '!CHANGE TO 4, IF FIVE TUBES IN BUNDLE 
IF Tact<1@ THEN 
Ibeg=lact 
Tend=lact 

END IF 
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NM PNM mM WN ry 


FOR I=Ibeg TQ lIend 


S~=@ 
Sy=@ 
Sxs=@ 
Sxy=@ 
Iset=@ 
ASSIGN @File TO File$ 
ENTER @File;Dto$,Itube,Eorat ,Odi ,Od2,0d3,0d4,0d5 
ENTER @File:Pc(*#),Emf(#),Tp(*) 
FOR J=@ TO 17 
TC J =FNTvsvieEmf( J )) 
NEXT J 
Tvap=(7T(@)+T( 1)4+T(2))/3 
Tlig=(7T(3)+T(4))/2 
Tliqgl=T(5) 
IF Isat=@ THEN 
Tsat=Tliagli 
BileSiE 
Tsat=Tvap 
ENO SIF 
Grad=FNGrad(TCI+6)) 
Tout=(T(2*14+1@)4+T(2*T411)9)/2 
Delt=Tout-T(I+6) 
Tavg=1(1+6 )+Delt+#.S 
IF I=Iend THEN 
Delt=T(1@)-T(Q9) 
Tavg=(T(10)-T(9))/2 
END IF 


Water/Ethylene Glycol Mixture Froperties 
Rhoeg=FNRhoeg( Tavg,Egrat ) 

Nueg=FNNueg( Tavg,Egrat ) 

Mueg=Nueg*Rhoeg 

Cpeg=FNCpeqg( Tavg,Egorat ) 

Keg=FNKeg( Tavg,Egrat ) 

Preg=Cpeg*Mueg/Keg 


Mdot=FNFmcal(1,T(1+6),Pc(1)) 
Veg=Mdot/(Rhoaeg*Ax ) 
Reeg=Veg*Di/Nueg 
Res=4*Mdot/(Mueg*#(PI*Di-4*Delta)) 
Qdot=Mdot+*Cpeg*Delt 

Qdp=Qdot/fAo 
Lmtd=Delt/LOG((Tsat-T(1+6))/(Tsat-T(1+6)-Delt )) 
Uo=Qdp/Lmtd 

IF Insert=@ THEN 

Bbl=.@668#(0i/L )*Reeg*Preg 
Bb2=1+.@4«((Di/L)*Reeg*Preg)*.5556 
Qmega=3.66+(Bb1i/Bb2 ) 

END IF 


OL 


N 
] 
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f. 


WW WNW bo! OF tl OW &G) DW Ww bl t! OF 6b) 
Io) I! tf) OF bt! OF IN Ot! Oy OF fb) OJ 
Mmsems NSsmsemngsennaeimn 


3482 
54@5 
341@ 
3415 
3420 
3425 
3430 
3435 
3440 
3445 
3450 
3455 
346@ 
3465 
3472 
2475 
348Q 
3485 
3490 
3495 


IF Insert=1 THEN 


Omega=(1+5.484E-3*Preg*.7*(Res/Hod)*1.25)%.5 


END IF 

IF Insert=2 THEN 
Omega=(Reeg*.65)*(Preg*.46) 
END IF 

IF Insert=3 THEN 
Omega=(Reeg*.76)*(Preg”.46 ) 
END IF 


R-114/R113 Properties 
Hf g=FNHfg( Tsat ) 
KFEFNKC TF) 
Rhof=FNRha( Tf ) 
Muf=FNMu( Tf) 


F=(Kf°S*Rhof *2*G*Hfo/( Muf *D0#Qdp ) )*,. 33333 
Ho=Alpa*F 
Two=Tsat-Qdp/Ho 
Tf=Tsat/3+2*Tuwo/3 
Y=(1/Uo-Rm )#F 
X=Do#F/( Keg*OQmega ) 
PRINT “OMEGA=";Qmega;: "F=";3F;"X="iX3 "Y="; VY 
Xa( Jset =x [ INEFFICIENT (MODIFY LATER® 
Ya(Jset )=Y 
S™=Ox+X 
Sy=Syty 
Sxs=SmuqtX#X 
Sxy=SxytX#y 
Jset=Jsettl 
IF Jset<Nsets THEN 3135 
ASSIGN @File TO * 
Slonpe=(Nsets*#Sxy-Sx#Sy )/(Nsets*#Sx5-Sx*2 ) 
Intcpt=(Sy-Slope*#Sx )/Nsets 
Cic=1/Slope 
Alpac=l/Intcpt 
Cere=neowce1-Ccic)7Cic?) 
Aerr=ABS((Alpac-Alpa)/Alpac ) 
IF Cerr>.@@1 OR Aerr>.@@1 THEN 
Alpa=(AlnpatAlpac)*.5 
er=(Cirseice)*_5 
PRINT "CIC=";:Cic; "ALPA=";Alpa 
6CQTO 31a@ 
END IF 
BEEP 
BEEP 
PROUNTER. 05 ei 
PRINT “CIC=";Cic; “ALPA="; Alpe 
Cia(Insert,I)=Cic 
Alpaa(I)=Alpac 


Te 


4500 
35@5 
351@ 
S515 
352@ 
gacs 
5530 
seo5 
3540 
3545 
3550 
S555 
356@ 
3565 
4570 
e275 
3580 
ese5 
359@ 
o595 
36@@ 
S6@S 
361@ 
Sa15 
3620 
Sa25 
363@ 
Sa05 
364@ 
3645 
365@ 
S655 
366@ 
S665 
367@ 
Se 75 
3680 
3685 
369 
3695 


——_ - be 


3720@ 


3705 


371@ 
aS 
e720 
e272 
373@ 
B7S5 
374@ 
S745 


PRINTER IS 701 
FOR J=@ TO Nsets-l 
OUTPUT @Xy;Xa(J),YatJ) 


NEXT J 
PRINTER IS 1 
NEXT I 
ASSIGN @Xy TO * 
SUBEND 


SUB Nusselt 
COM /Nus/ ITin,Tsat,Qdp,Hoc,Kf,Rhof ,Hfg,Muf,Do,Itube 
Do=.0159 
Ho=100@ 
IF Tin=@ THEN 
BEEP 
INPUT "ENTER TSAT AND HEAT FLUX" ,Tsat ,Qdp 
ENO IF 
Hfo=FNHfg(Tsat ) 
Two=Tsat-Qdp/Ho 
Tf=Tsat/3+2*Two/3 
KfF=FNKC TF) 
Rhof=FNRho( Tf ) 
Muf=FNMu( Tf ) 
iee=eeeoe(t hf o*knes 2*4 8iefa/(Muf*00o*Qdp))”.3333353 
IF ABS({(CHo~-Hoc )/Hoc )>.@@1 THEN 
Ho=(HotHac )*.5 
GOTO 3585 
END IF 
IF ITin=@ THEN PRINT “HO=";Hoc 
SUBEND 
SUB Plotl 
DIM Yaa(4) 
PRINTER IS 70S 
Tdv=1 
BEEP 
INPUT "OK TO USE DEFAULT VALUES (1=DEF=Y,Q@=N)" ,Idv 
IF Idv=1 THEN 
Ttn=2 
Xmin=1 
Xmax=5 
Xstep=l1 
Ymin=@ 
Ymax=2.@ 
Ystep=.5 
ESE 
BEEP 
INPUT "ENTER MINIMUM AND MAXIMUM X-VALUES" ,Xmin, Xmax 
Beer 
INPUT "ENTER MINIMUM AND MAXIMUM Y-VALUES",Ymin, Ymax 
BEE 
INPUT “ENTER STEP SIZE FOR X-AXIS" ,Xstep 


18, 


475@ 
yaya 
376@ 
S755 
377Q@ 
aS 
378@ 
e735 
379@ 
eteeke 
38QQ 
38@5 
3810 
3815 
5820 
3825 
3830 
S455 
3840 
3845 
385@ 
S255 
S386Q 
3865 
38702 
Ba/75 
388 
3885 
3B9A 
Seco 
39@@ 
3905 
eal@ 
6915 
eae 

S425 
ga50 
e355 
394@ 
3945 


Saou 6 


S355 
3962 
B205 
397@ 
S272 
3989 
Bess 
S300 
So 26 


BEEP 
INPUT "ENTER STEP SIZE FOR Sasa) omer 
ENO IF 
PRINT "IN: SP1;IP 23@0,1800,830@,680@; " 
PRINT "SC @,100,0,1@@;TL 2,@:" 
Sfx=1Q0/( Xmax-Xmin? 
Sfy=10Q/0Ymax-Ymin) 
PRINT "PU @,@ PD" 
FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin)#Sfx 
PRINT "PA";sX,",@3 X73" 
NEXT Xa 
PRINT "PA 10@0,@;PU; “ 
PRINT “PU PA @,9 PD" 
FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)#Sfy 
PRINT PALO 
NEXT Ya 
PRINT = PAnG,IOGe Te a 2” 
FOR Xa=Xmin TQ Xmax STEP Xstep 
X=(Xa-Xmin)#*Sfx 
PRINT PA 7eeeeelOO. XT” 
NEXT Xa 
PRINT “PA 100,100 PUNFA 100.0 Fae 
FOR Ya=Ymin TQ Ymax STEP Ystep 
Y=(Ya-Ymin)#Sfy 
SIS UE eG) leit Wale eee IE 
NEXT Ya 
PRINT "PA 10@,1@@ PU" 
PRU Nie eGo eee 
FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin)#*Sfx 
PRINT "PA";X,",@3° 
PRINT Ce =22 =e Bee ee 
NEXT Xa 
PR ee ee eee 
FQR Ya=Ymin TQ Ymax STEP Ysten 
Y=(Ya-Ymin)#Sfy 
SUSE ete We ei 
PRIN) “GP = 2S eee ya 
NEXT Ya 
sels |= 
INPUT "SELECT MODE (@=HN/H1,1=HN(avg)/H1)",Ism 
Ism=Ismt+l 
IF Idv=1 THEN 
IF Ism=1 THEN Ylabel$S="HN/H1 " 
IF Ism=2 THEN Ylabel$S="HN(avg)/H1" 
Xlabel$="Tube Number" 
Eas 
GEEP 
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4QQ@ 
4@@5 
4Q1@ 
4@15 
4Q2Q 
4025 
4230 
4@35 
4Q4@ 
4@45 
A@5@ 
4@55 
4262 
4@65 
4272 
4@7S 
4@8@ 
4085 
4@9@ 
4@95 
4122 
41@5 
4112 
4115 
412@ 
4125 
4132 
4135 
4142 
414% 
415@ 
4155 
416@ 
4165 
4172 
4175 
4182 
4185 
419@ 
4195 
4222 
42@5 
421@ 
4215 
4220 
4225 
425@ 
4235 
4242 
4245 


INPUT "ENTER X-LABEL" ,Xlabel$ 
BEEP 
PNPU) eee Niliekey-LObeL” ,YlabelS 
END IF 
Cees eee e eer rn S610 CP 7 =-LEN(XlabelS)/2; "0:LB™xXlabel$:"" 
Brel eee retour  e, 3 -LCEN( Ylabel$)/2+*S/6; "DI @,1;LB"; Ylabels:;°* 
PRINT “CP @.@ DI" 
Okp=1 
SGEEP 
INPUT "WANT TO PLOT DATA FROM A FILE (1=DEF=Y,@=N)?7" ,0Okp 
IF Okp=1 THEN 
BEEP 
eve oeNen ste NAME OF THE DATA FILE” ,Dfiles 
ASSIGN @File TO Ofile$ 
BEEP 
INPUT “ENTER THE BEGINNING RUN NUMBER" ,Md 
BEEY 
INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED" ,Nsets 
BEEP 
mee SBEECT A SYMBOL FOR THE PLOTTER (1=*, 2=+, dec ,4=0,S5=*)", Sy 
PRINT “PU OI" 
IF Sy=1 THEN PRINT "SM#" 
eo —2  IheEN PRIN] “SM+~ 
IF Sy=3 THEN PRINT "SMc" 
IF Sy=4 THEN PRINT "SMo" 
tess =o THeEN PRINT "Sh" ” 
FOR T=1 TO Nsets 
FOR J=O0510 4 
ENTER @File: Yaa(J),D 
IF J=Q@ THEN Ytop=Yaa(Q) 
Yaa(J)=Yaa(J)/Ytop 
NEXT J 
FOR J=@ TQ 3 
X=(J+1-Xmin)#Sfx 
Y=(Yaal( J )-Ymin)*Sfy 
Re Nee rule ak. Yc 
NEXT J 
NEXT I 
BEEP 
ASSIGN @File TO * 
GOTO 4042 
END IF 
BEF . 
INPUT "LIKE TO PLOT THE NUSSELT RELATION (1=Y,@=N)7" ,Qknus 
PRN Us oh 
TF Oknus=1 THEN 
FOR Xa=Xmin TO Xmax STEP Xstep/5@ 
X=(Xa-Xmin )#*9Of x 
IF Ism=1 AND Xa>Xmin THEN Ya=Xa*.75-(Xa-1)*. 75 
IF Ism=2 AND Xa>Xmin THEN Ya=Xa*(-.25) 


SS 


4258 IF Xa=Xmin THEN Ya=1 


4255 Y=(Ya-Ymin)*Sfy 

426@ PRINT Ge reer aml 

4265 NEXT Xa 

427@ BEEP 

4275 PRINT SEG 

4280 INPUT "MOVE THE FEN TO LABEL THE NUSSER IEE wer 
4285 PRINT “LBNusselt”" 


4290 an {@) LI 
4295 IF Ism=2 THEN 


430@ slelele 

43@5 INPUT "LIKE TQ PLOT EXPTL CURVE (1=Y,Q@=N)" ,Qkex 
431@ Nq=@ 

4315 IF OQkex=1 THEN 

4320 BEGI 

4325 INPUT “ENTERD THEVEXPONENT. Je™ 
4335@ FOR Xa=Xmin TO Xmax STEP Xstep/1@ 
4335 Nq=Nq+l 

434@ Ya=Xa*(-Ex ) 

4345 X=(Xa-Xmin)#*Sfx 

435@ Y=(Ya-Ymin)#Sfy 

4355 IF Na MOD 2=@ THEN 

4360@ UNM IRA pong iis a 

4365 EESE 

437@ Cree ty eee el 

A475 SNe te 

438Q NEXT Xa 

43385 PRINT "PU" 

439@ BEEP 

43595 INPUT “MOVE PEN 70 LABEL AND HIT ENTER =o: 
440Q PRINT "LBs=@" 

4405 PRINT “PR -1 @" 

441@ PRN ee ene 

4415 GOTO 4300 

447Q ENG” TF 


4425 END IF 

4442 GOTQ 4510 

4435 REEP 

4440 INPUT “LIKE 70 PLOT KERN RELATIONSHIP (l=) (O=Nog ay ce 
4445 We ieee) Ula 


4450 | FOR Xa=Xmin TQ Xmax STEP Xstep/2@ 

4455 Ya=Xa*(-1/6) 

446@ X=(Xa-Xmin )#*Sfx 

4465 Y=(Ya-Ymin)«#Sfy 

4470 PRINTS | At. aye 

4475 NEXT Xa 

4480 PR tee Ge 

4485 slsieie 

4492 INPUT *MOVE THE PEN 20 LABEL KERN REG RON Sree or 
4495 PRINT "LBKern;PU" 
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4522 
45@5 
451@ 
4515 
452@ 
525 
453@ 
4535 
4540 
4545 
4550 
4555 
4562@ 
4565 
457@ 
4575 
458@ 
4585 
4590 
4595 
462@ 
46@5 
4610 
4615 
462@ 
4625 
4620 
4635 
4642 
4645 
465@ 
4655 
4662 
4665 
467@ 
4675 
468Q@ 
4685 
469@ 
4695 


47@@ | 


4725S 
471@ 
4715 
472@ 
4725 
473@ 
4735 
4742 
4745 


END IF 
PRINT "PU PA @,0" 
Beer 


PPT Sie Oe PLOT SEISSENBERG RELATION (1=Y,@=N)?", Oke: 


IF Okei=l THEN 

FOR Xa=Xmin TO Xmax STEP Xstep/1@ 
Ya=.6+.42*Xa%(-.25) 
X=(Xa-Xmin )#Sfx 
Y=(Ya-Ymind«#Sfy 
CRuNieeetrm .X,y, PR" 

NEXT Xa 

PRINT “PU" 

Beer 


pie noe Tae PEN T@ LABEL THE EISSENBERG LINE” , Ok 


PRINT “LBEissenberg:; PU" 
(0) hs 
PRINT "PU SPQ" 
SUBEND 
DEF FNPvst(Tc) 
wome/Fid/ Ift 
DIM K(&) 
IF Ift=@ THEN 
BEEP 
PRINT "PUST CORRELATION NOT AVAILABLE FOR R-114" 
Sor 
fe 5) Ihe 
TE eiit=! THEN 


Wages? se 3'254554 ,-26.08025696 ,-168.17@6546,64.23285504 ,-118. 9646225 


Pole dete iy oe. 2, 4/506 765,1.69,6 

READ K(*) 

ieee +2 /5.15)7647.5 

Sum=@ 

FOR N=@ TO 4 
Sum=oum+tkK ( Noee¢ 1-7 °C N+ i) 

NEXT N 


Smesum ie Glth (5 )*( l-l)thCa etl 7 92) = ai (RG) #el—T )“2+KC 8) 


Pe=cAr (Br) 
P=22120@@@*Pr 

FNO IF 

IF Ift=2 THEN 
Mr=Te*! .8+52+459.6 


Megan 555—4350.98/1F-9.263S5*LGT( Tf +2 .@539E-S3amt ) 


P=P*1913525/14.696 

END IF 

IF Ift=3 THEN 
G@adeocebos-506G.1/(ict273,15) 
P=133.38241@°A 

ENG 

RETURN P 

F.NEND 


Md 


475@ 
4755 
476@ 
47765 
477@ 
4775 
4782 
4785 
4793@ 
4795 
48@@ 
48@5 
481@ 
4815 
4816! 
4820 
4825 
483@ 
4855 
484@ 
4845 
485@ 
4855 
4860 
4865 
487@ 
4875 
488 
4885 
489@ 
4895 
49@a 
49@5 
491@ 
4911 
4913 
4ai5 
4920 
4322 
4932 


A555 . 


494@ 
4946 
4952 
4955 
496@ 
4965 
497@ 
AS 
AQ98@ 


DEF FNHfg(T) 

COM /Fld/ Ift 

IF Ift=@ THEN 
TF=T#1, +32 
Hfg=6.1451SS8E+1-Tf#(6.9S1Q79E-2+TF*(1.39886R8E-44+1 . 9GQ7R843E-7*TF)) 
Hf g=Hfg*2326 

END IF 

IF Ift=1 THEN 
Hf 9=2477200-2450*(T-10) 


ENO IF 
IF Ift=2Z THEN 
TF=T#1.8+32 


Hfg=7.0557857E+1-Tf*(4.848@052E-2+1,2619048E-4*T Ff ) 
Hf g=Hfg*2326.1 
Hfe=(iebli-.@031*1 er Gets 

ENO IF 

IF Ift=3 THEN 
Tk=T4+2732e5 
Hfg=1.35264E+6-Tk#(6.38263E+2+Th*.747462 ) 

ENG site 

RETURN Hfg 

FNENO 


~ DEF FNMu(T? 


COM 7Fid7 att 
IF Ift=@ THEN 
Th=T+273.15 
MuSEXPS=>4 46364101147 /0 wis 
END Sita 
IF Ilft=1 THEN 
A=247.8/( 1 tts 
Mu=2.4E-5*10°A 
END IF 
IF Ift=2Z THEN 
Mu=8.9629819E-4-T*(1.1094609E-S-71*5,. SBESZ9E se» 
Tete. lS 
Mu=1].34E-5*10.0°(5035.0/7( Th ee)? 
ENO IF 
IF Ift=3 THEN 
ely lees. toe) 
Mu=EXPC-11.0179+7TK*( 1. 74dE+S3-TR*CZ: 803558 tS aie een 
ENO IF : 
RETURN Mu 
FNEND 
DEF FNUvst(Tt) 
COM /Flid/ Ift 
IF Ift=@ THEN 
BEEP 
PRINT "VYUST CORRELATION NOT AVAILABLE FOR R-114" 
50g 
ENO IF 
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4985 TE Litt HEN 


4990 P=FNPvst(Tt) 

a5 W=ttF275.15 

5@aa@ X=15Q0Q/T 

5@@5 ply Gliese tl c=4) 

5010 ele Gk) 2 Dee XP(X )/X*.5S 
5@15 B=.Q0@15*F1-.@@@0942+*F2-.@0@4882+*x 
5020 K=2*P/( 461.527 ) 

5825 VEC 14+(1+2*Rek 0°. 5)/K 


5030 Ss 
5@35 IF Ift=2 THEN 


5042 te =Tt*1. 8+52 
5@45 Seems oScov athe, (G12 /7262—1F*5. 172619@E—-4 >} 
5@5@ V=V/16.@18 


5@55 ENDIF 
S2@6a IF Ift=3 TREN 


5@65 Pe=1t+276.15 
527 P=FNPvst( Tt) 
sy75 Wao, oo*1k/P 


528 END IF 

5@e85 RETURN V 

5@3a FNEND 

5095 Den eENGp<{T > 

51@@ Ghetto, Ift 

51@5 IF Ift=@ THEN 

511@ emt ee (S215 

SNS Cpo=.4@118+Th#(1.65@0Q7E-S3+Tk*#(1.51494E-6-Tk*6.67853E-1@) ) 
512@ SN)" 0s 

ali IF Ift=1 THEN 

513 Pomme oesa | = 2 eeecoeoo71*( 4, A4cS3b1E-S+2 ria cce-7*T ) ) 
e355 ENDIF 

514@ IF Ift=2 THEN 

51451 Pg@ececata 2 Seal tl * (9, 6400455E-4+1.720779ZE-B*T ) 

5146 Cp=(929.0+1.03*T )/1000. 

5159 END IF 

Bioo IF Ift=5 THEN 

516 PrsT+2735.15 

S165 Comin l6oo*( 1), 6884E—-2+7hs(3.35083E-S-Tk#(7.224E-6-Tk sie G1 748E-Seme) 
517@ Nels 

5175 RETURN Cr*#1Q@QQ 

518@ . FNEND 

ees DEPSENRhe tT) 

519@ COME ldy a1 ft 

S135 TF lft =O HEN 


520@ Tk=T+273,15 

S2@S R= eh 755.95) 

5210 Ro=36.32+61.146414*X*%( 1/3 )+16.418015*X4+17.476858#X* .S+1.119828#x*2 
2 1s Ro=Ro/.@62428 


5220 ENE 
eee IF Ift=#l THEN 


te, 


5230 
5255 
5248 
5245 


5246 |! 


5258 
5225 
5269 
5765 
5270 
S22 
5280 
S205 
5290 
5/4 ote 
5380 
BSUS 
531@ 
Zol9 
=320 
2525 
5330 
B55 
5340 
S45 
535@! 
B45 1! 
Bao | 
555 
5354 
SSe5 
Daa / 
536@ 
Se )shs 
5370 
S575 
538@ 
Seas 
5390 
podS 
5420 
54@5 
5410 
S415 
542@ 
9425 
5430 
5445 
5442 
5445 


Ro=999 .S2946+1T*(  O@1269-14(5. 4825.5 -s-(*1 > 740 eeoee 
END ke 
IF Ift=2 THEN 
Rom! 62074796 +5—=14(2 (715654040222 oa7e coe —o 
Ro=1./((251/+. 00054 741° 101) ste a 
END IF 
IF Ift#3 THEN 
PREt+t2 73 spo Sees 
Vf=G,.2484BE-447k*(5 . 2795E-74+1* *< 3244s eo oe eee) 
Ro=1/VF 
END ir 
RETURN Ro 
FNEND 
DEP VENEG Gls) 
Pr=FNCp(T )*FNMu( T D/FNKCT ) 
RETURN Pr 
FNEND 
DEF FNK(T) 
COM el at 
IF Tft=@ THEN K=.071-.Q@0261*T 
IF Ift#l THEN 
REC T+*273- 15) 275715 
Ke-,922474+X%"(2,.8395-X4+(1  80@7—xX*( 5257 (= 1754s oD 
END IF 
LP Pi is =2 en 
K=8.2095238E-2-T#( 2.22 14286E-4+T*2.3809524E-8 ) 
K=,@802- .@00203*7T 
K=.071-.Q000261*7 
ef aT 21 ease 7 
K=4,84619@5E-2-TF*6.5714286E-S 
K=K*1.73@8 
END IF 
IF Ift=3 THEN 
TReT+Z275.i5 
K=4,1868E-4*(519.442+.3520920*Tk ) 
EIN ite 
RETURN k 
FNEND 
DEF FNHf(T) 
CONS Eras iit 
IF Ift=@ THEN 
BEER 
PRINT “HF CORRELATION NOT FOR R-114" 
SIGE 
ENG ate 
IF Ift=1 THEN 
Hf=t*(4.203849-T#(S.88132E—-4-1*4, 5516051 7E=8)) 
ENG Ie 
TPetit=z THEN 
Lea ga Ks 
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5452 
5455 
S46@ 
S465 
5472 
$475 
548@ 
S485 
5499 
5495 
5500 
fe0S 
5510 
BaS 
5520 
Boe. 
Sel 
SS 
554@ 
6545 
5550 
S555 
5562 
Sa5 5 
557@ 
BS/5 
5589 
Saco 
5599 
= haga 
5620 
56@5 
5610 
Sbis 
5620 
BoZ0 
563@ 
BbSo 
564 
S645 
565@ 


BESS 


556@ 
5665 
567@ 
5575 
S68@ 
S685 
569@ 
5695 


HfF=8.2078S71+TF*(. 19467857+T£*#l .5214286E-4 ) 
HfF=HF #2. 326 
ENO IF 
Pe. it —s) een 
Ht —2o0 en eee e VERIFIED 
Sina) Us 
RETURN Hf*#10@Q@ 
FNENO 
S sla (2 lake Z 
COME Orel  otar, sym, icon 
ihe? eae 1 ft 
DIM €(9),Xya(7),Doa(4) 
DATA @.0@158,0.0158,0.0158,@.0158 
READ Doa(*) 
Fu=l 
PRINTER IS 1! 
SlelSir 
PRINT USING "4X,""Select Option X-Y Limits:""" 
PRINT USING “6X,""@ Use default values*"" 
PRINT USING “6X,""1 Use new values""" 
INPUT Qkd 
BEE 
INPUT "ENTER TUBE CODE",Icode 
Do=DoalItube) 
Tht=2 
REEP 
PRINT USING "4X,""Select option:""" 
PRINT USING "6X,""@ h versus g""" 
Ee UNieUSINGS@=Sx, 1° o versus Delta-T*"" 


PRINT USING "6X,""2 h versus Delta-T (default)""" 


INPUT Iht 
PRINTER IS 705 
IF OQkd=0 THEN !AXIS DEFAULT VALUES 
IF Iht=@ THEN !¢(h vs q) 
Ymin=@ 
Ymax=6Q 
Ystep=1@ 
Xmin=.2 
Xmax=1.4 
Xstep=.2 
= Gs 
IF Iht=1 THEN !(q vs t) 
Xmin=Q 
Ymin=@ 
Ymax=.5 
Xmar=15 
Xstep=3 
Ystep=.1 
ENOQ IF 
Poeent—=calheEN '{h vs t) 
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ro 


nonin own 
DO wo 


ID 
> p HD! I 


ng sy 


Xmin=@ 
Ymin=@ 
Xmax=5@ 
Ymax=6 
Xstep=12 
Ystep=] 
END IF 
ENG. te 
IF Okd=1 THEN 
BEEP 
INPUT “ENTER MINIMUM AND MAXIMUM X-VALUES",Xmin,Xmax 
BEE 
INPUT “ENTER MINIMUM AND MAXIMUM Y-VALUES" ,Ymin, Ymax 
BEEF 
INPUT “ENTER STEP SIZE FOR }-Axto. 3 xcuen 
BEEF 
INPUT “ENTER STEP Si2E FORMY Ax ts cen 
END IF 
BEEP 
PRINT “IN; SP1;51P 235007 RGU easUG eau 
PRINT “SE 0.100 G2 100, jieec. es 
Sfx=100/( Xmax-Xmin) 
Sfy=100/(Ymax-Ymin) 
BEEP 
Tcg=@ 
INPUT "LIKE 70 BY-PASS CAGE f=) {0=\-DERAUI ee etca 
TF leo=)l THeNeal7S 
PRINT ae OU Orig 
FOR Xa=Xmin TO Xma» STEP: Xstep 
X=(Xa-Xmin)*#Sfx 
PRUNt Seeesn | Os kee 
NEXT Xa 
PRINT “PA 100,0:FU;~ 
PR UN Oe pee cie tay 
FOR Ya=Ymin TQ Ymax STEP Ystep 
Y=(Ya-Ymin)#Sfy 
PRINT “PA @,":Y,"YT" 
NEXT Ya 
PRINT ~FA @,10@ TL @ 2: 
FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmind#Sfx 
BRUNT “PA>; Xo | CO exe 
NEXT Xa 
PRINT “PA 100,100 PU PA 100,@ PD" 
FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)«#Sfy 
minha a cHO ems alates CIES 
NEXT Ya 
PRINT "PA 10@,1@@ PU" 
PRINT "PA @,-2 SR 1.5,2" 
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5952 
Sie) 
5962 
S255 
oo 40 
Bo 72 
598 
Bus 
§939@ 
Seeks! 
B2aa 
6205 
6212 
69015 
6922 
6025 
623 
6035 
624 
6045 
6250 
6955 
6262 
6@65 
6272 
6275 
6@8a 
6285 
6292 
6@95 
6122 
6105 


kW/m" 


611@ 
ee15 
612@ 


0,0. 


6125 
61322 
Bio 


wee 2s PR 


6142 
6145 


2,@ PU: 


615@ 
bos 
6162 
6165 
6179 
S175 


FOR Xa=Xmin TO Xmax STEP Xstep 
K=(Xa-Xmin)#*Sf 
PRIN @@rob sk Os 
Teethers OOTHEN PRINT "CP -1.5,-1:.80:PR -1,0:LB";Xa:"" 
IF Xe=2> TREN PRINT "CP -.5,-1;:L80" . 
Xin=@ 
IF Xa MOD 1=@ THEN Xin=1 
IF Xa>=1@ THEN PRINT “CP -2,-1;LB";Xa;"" 
IF Xa>i AND Xa<i@ AND Xin=1 THEN PRINT “CP -1.25,-1;LB";Xa;"" 
IF Xa>l AND Xin=@ THEN PRINT “CP -2,-1;LB";Xa;“" 
[pe aSteineNeeR IN’ “€P -1,-l:Leh 0° 
NEXT Xa 
PRINT "PU PA @,@" 
FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)#*Sfy 
Reh eee ieee Yous 
IF Iht=@ AND Ya>@ THEN 
IF Ya<l AND Ya<>o@ THEN 


PRINT "PR 2,Q" 
PRINT "CP -4,-.25;LBQ;PR -2,@;LB";Ya:*" 


IF Ya=@ THEN PRINT "CP -2,-.25;LBQ" 
eect Oi nite erehee RINT 9 CP —S,—,25;LB°3¥a;”’* 
Peaveas teat oelint=ZetroN PRIN) “CP —4,—-,25;LB";Ya;"* 
teen Cero yas. l0 AND Iht=2 THEN PRINT “CP =3,.—.25:LB":Ya:"* 
teed ie Pee werR INT = GP —4.—.25;LB1.0° 
NEXT Ya 
IF Okd=1 THEN 
Beer 
INPUT “ENTER X-LABEL” ,Xlabel$ 
Beer 
NEU SeeENTER yY-EOBEL”  Ylabels 
END IF 


Peat eT heN 
Cees nee rere -12,59;01 O,1;LBRPR 140.5; EBo,PR -1,@.S 1887 ( 


aoe eed Soe lS; LB27SR 1,5,2:PR 25. .5;L8.2PR 
ELSE 

Cee eee eee Ss Oi OB gy (MW Ms oR belle ower —1..0,5;5L62;SR 12a.2;P 
Sve iLjs) ve 
END IF 
IF Iht=@ THEN 

PRINT "SR 1.5,2;PU PA 40,-12;D1:LBqg/(MW/m;PR @.5,1;SR 1,1.5;LB2;5R 

SS ei seh 

Bao 

Poiana — | Oiee( 1 sPRoe 5-1, bbs; PR 
ere 

PR UN eeleeWwios FR 
END IF 
eel N ieee ewe | © 
Xlg=1.E+6 
AiG =— tere 
Xal=50 


-©,0; CBR 


Sa Cham 2.4, 5 eee 


So jee 
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618@ 
SIRS 
619@ 
eV ESE 
B2a@ 
62@5 
6212 
S215 
6220 
ae ae) 
6230 
6235 
624@ 
6245 
625@ 
S255 
626@ 
EZoS 
6270 
6275 
6280 
6285 
6290 
S232 
63a0 
63@5 
631@ 
5) a) es 
6328 
SSZo 
6350 
Be55 
634@ 
6345 
635@ 
Beso 
6356@ 
6365 
6370 
eos 
B58 
6385 
6392 
So25 
6420@ 
64@5 
641@ 
6415 
6420 
6425 


Yal=9%5 


Nrun=@ 
BEEP 
INPUT "WANT FO PEON DATA FROM A Pi Be Gl ee)? eee 
X1ll=1.E+6 
Xul=-1,.E+6 
Okp=2 
IF Ok=1 THEN 
BEEP 


INPUT "ENTER THE NAME OF THE SPEGISDATh er ee Seat eues 
ASSIGN @File TO DO_file$s 

IF Icomb<>@ THEN 6265 

Sx=@ 

Sy=@ 

Sx2=0 

Sxy=@ 

Md=1 

BEEP 

INPUT "ENTER THE BEGINNING RUN NUMBER (DEF=1)" ,Md 
Npairs=9 

BeBe 

INPUT "ENTER THRE NUMBER OF X—Y PAIRS STORER SCMEE-—4)) Sheena. 
Nrun=Nrun+Npairs 

re nS RE Sinema lls | 

BEEP 

PRINT USING "4X,""Select a symbol:""" 


PRINT USING "4X,"" 1 Star 2 Eee mns dec) 
PRINT USING "4X,*" 3 Circle 4 Square"”” 
PRINT USING "4X,"" S Rombus""" 

PRINT USING "4X,"" 6 Right-side-up triangle""" 
PRINT USING "4X,"" 7 Up-side-down triangle""" 
INPUT Sym 

BEEP 


INPUT “ENTER TUBE NUMBER FOR PLOTTING (Q=TOP,1=SECOND,...) (ime 
PRINTER IS 7@S 
IF Sym=1 THEN PRINT "SM#" 
IF Sym=2 THEN PRINT "SM+" 
IF Sym=3 THEN PRINT "SMo" 
LPM te TREN 
PQPsst=t 700 Mcgee 
ENTER @File:Xya(*) 
NEXT I 
ENGr 
FOR I=1 TO Npairs 
ENTER ®File;:Xya(+*) 
Ya=Xya(Itube*2z ) 
Xa=Xya(Itube*2+1 ) 
Yo=LQG( Xa) 
Xc=LOG(Xa/Ya) 
Sx=Sx+Xe 
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(=o> 


643@ SY=SV+ iG 

64355 SxZmoret nc 2 

6440 Sxy2=SxytXc#Ye 

6445 IF Ihnt=@ THEN 

6450 Xt=Ya 

6455 Ya=Ya/Xa 

646@ Xa=Xt 

6465 IF Xa¥1.6€+6>Xu1 THEN Xwl=Xa/1.E+6 

6472 IF Xa/1.E+6<X11 THEN X11=Xa/1.E+6 

6475 ea Mer Ie 

6482 IF Iht=0 THEN 

6485 X=(Xa*il.E-6-Xmin)*Sfx 

6490 Y=( Yael] .E-3-Ymind«Sfy 

6495 END IF 

6500 IF Iht=1 THEN 

6505 X=(Xa-Xmin )#*Sfx 

651@ f= Ya*!.E-6G-Ymin ?*Sfy 

6511S ENO IF 

652@ IF Iht=2 THEN 

eves! X=(Xa/Ya-Xmin d#Sf 

653@ V=( vane oyun) * Sty 

Sogo ENO IF 

6542 IF Y>100 OR Y<@ THEN 6585 

6545 If Symos TREN PRINT °SM” 

655@ ie Sym<4 TREN PRINT °SR 1,4,2.4” 

Sys! PRUNT SPAT Ex. y¥ © 

656@ ie Sm See NPeR INT “SR t.2,1./6- 

6565 IF Sym=4 THEN PRINT eee es) Ue eae 
6570 IF Sym=S THEN PRINT pe als (Sila 4 Sige 

eoy > IF Sym=6 THEN PRINT “UC@,5.3 138, ae ale, 4 56 ; 
6580 IF Sym=7 THEN PRINT "UC@,-S. oa. Poor: — Seo, 
Bou NEXT I 

6590 Bee 

6595 TNPUR WANT TQ LABEL ¢(1=Y ,@=N)7" fiero! 

66aG IF Tlbl=i1 THEN 

6605 tee syne o THEN PRINT Sh 

6610 i sym 4.7THEN PRINT “SR 1 J7a2 .4° 

Bes PRIN anXal.Yal, is 

662@ IF Sym>3 THEN PRINT "SR 1.2,1.6" 

6625 TE Sym=4 THEN PRINT “UCG2,4,99,0.-2,-4,9,0,8,4,@:;' 
6632 | IF Sym=S THEN PRINT aes See eee, 5, 6: 
6635 IF Sym=6 THEN PRINT "UC@,S.3,99,3,-8,-6,8,3,8:;" 
664@ IF Sym=7 THEN PRINT "UCQ, -5, 5. ae — Seen ome.-o.-c> 
6645 PRoNi» of” 

665@ [eee mcs HEN FPRUNT "PR 2,0 

S655 av llieeatre 2 )-1.0;5R PoOeiyes te re files: 

666@ Yal=Yal-5 

6665 REEP 

667@ INPUT “WANT TO ADD ANOTHER STRING (1=Y,@=N)?",Tas 
6675 IF ITas=l1 THEN 
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668@ 
6685 
669 
6695 
6722 
67@5 
671@ 
S715 
672@ 
6725 
673@ 
5735 
674@ 
6745 
5750 
S/o. 
676@ 
B7ao 
6772 
Bio 
5780 
6785 
679@ 
2) (fale 
6800 
68@5 
6819 
6815 
5820 
BBLS 
683@ 
6835 
6842 
6845 
685@ 
Bes 
5860 
6865 
687@ 
6875 
688@ 
6885 
6890 
yeaa: 
592@ 
69@5 
691@ 
Scio 
6920 
Stel 


BEEP 
INPUT “ENTER THE STRING" ,LabelS 
PRINT “PR 2,0;S5R 1.0,1.93:308 j babeleres 
GOTO 6665 
ENG Tr 
END IF 
BEER 
INPUT “WANT TQ COMBINE ANQTHER FILE? (1=Y,@=N)",Icomb 
ASSIGN @6File TQ * 
X11=5 
Xul=45 
IF Icomb<>@ THEN 62290 
BEEP 
INPUT "WANT TQ PLOT A LEAST-SQUARES LINE (1=Y,90=N)",I1s 
IF Tls=1 THEN 
BEEP 
INPUT "SELECT EXPONENT: @=COMPUTE, 1=@.75)",Iexp 
BEEP 
INPUT “SELECT CURVE TYPE (@=SOLID,1=DASHED)" ,I1t 
T1t=lTit+i 
PRINT. or. 
IF ITexp=@ THEN 
Bh=(Nrun*Sxy-Sy*#S~« )/(Nrun#Sx2-S*%Z ) 
Bese 
Bb=.75 
END IF 
Aa=(Sy-Bb*Sx )/Nrun 
Aa=EXP( fia) 
CR UNG eel el 
PRINT USING "10X,""a = “",2Z.4DE"; Aa 
PRINT USING “10X °°n = “32. 40G2E5 
PRINTER IS 7985 
Tn=Q 
IF Iht=@ THEN Xxstep=Xstep/4@ 
IF Iht>@ THEN Xxstep=Xstep/19 
FOR Xa=X11 TQ Xul STEP Xxsten 
IF Xa>.99*Xmax THEN 6995 
IF Iht=i1 THEN Ya=Aa*Xa*Bb 
IF ITht=@ THEN Ya=Aa“(1/Bb)*(Xa*l1.6+6)°((Bb-1)/Bb) 
IF Iht=2 THEN Ya=Aa*Xa*(Bb-1 ) 
IF Iht=@ THEN 
Y=(Ya#l.E-3-Ymin )#*Sfy 
X=(Ka-Xmin)#*Sfx 
ENO IF 
IF Iht=1 THEN 
Y=( Yael. E-6-Ymin)#*Sfy 
X=(Xa-Xmin)*#*Sfx 
END IF 
IF Iht=2 THEN 
Y=( Yar lvb-o-1m1 mo eony 
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693@ 
)is) 2s 
6942 
6945 
695@ 
eel =| 
6960 
goes 
6970 
6975 
6980 
B9@5 
6992 
eho 
7002 
7@Q@S 
701 
7@15 
7220 
7825 
7230 
7035 
7042 
7045 
7250 
7055 
726Q@ 
7@65 
7272 
7@75 
7282 
7085 
709@ 
7@95 
71@@ 
71@5 
7110 
mils 
712@ 
fae 5 


lei 


isS 
714@ 
7145 
715@ 
7155 
716@ 
ABS 
717@ 
E175 


X=(Xa-Xmin)#*Sfx 
END IF 
IF Y<@ THEN Y=@ 
IF Y>1@@ THEN 6990 
PP it=1> THEN 
Baie fay ek. Ya PO" 
ei Lote 
In=Intl 
Ir=In MOD Il1lt 
Peete ne Nee RINT "PAR 3k, Y, “PD” 
foete- Oe teen PRINT “PA .X,Y, "PU" 
eiMED = 
NEXT Xa 
RN Pu 
ee) iF 
Icomb=@ 
Saige olss 
ENG re 
Penner U Sie 
BEE 
INPUT "WANT TQ PLOT NUSSELT LINE (1=Y,@=N)?",Inp 
IF Inp=@ THEN 7125 
Beer 
PNG ewe NER eons “SEFAULT=18 DEG ©)” ,Tsat 
Hf g=FNHfo(Tsat ) 
Xi]=5 
Xul=45 
FOR Xa=X11 TO Xul STEP Xstep/5@ 
Tfilm=Tsat-Xa*.5 
Kf=FNK(Tfilm) 
Rhof=FNRho( Tfilm) 
Muf=FNMu( Tfilm) 
Ya=.728#(Kf*S3*Rhof*2*9,.81*Hfo/(Muf*Do#Xa))*.25 
X=(Xa-Xmin )#Sfx 
P= yeti ce SayYmin tof y 
CRU eee Amaia, PD 
NEXT Xa 
Eran “PU PA @,0" 
PRINT "PU PA @,@ SPQ" 
SUBEND 
SHBSP lots 
COM /Dri/ Star,Sym,Icon 
COlmaaas 1 ft 
OiM Gc4> 
Fw=l 
mallee abe al 
BEEP 
PRINT USING "4X,""Select Qption X-Y Limits: """ 
PRINT USING "6X,""@ Use default values""" 
PRINT USING "6X,""1l Use new values""" 
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7180 
(Gers 
719@ 
7 leks 
7220 
72@5 
7218 
WAS 
7220 
Uj (aflage 
7230 
W255 
7240 
7245 
7250 
{are 
726@ 
W255 
ti10 
G25 
7280 
7285 
7290 
Geos 
7580 
73@5 
731@ 
(215 
73520 
(25 
7330 
w355 
7542 
7345 
72350 
mos5 
7368 
7365 
7370 
Tat 
7380 
vecS 
Goa? 
Hoos 
7420 
7405 
741Q 
HAs 
7420 
7425 


INPUT Okd 
PRINTER IS 705 
IF Okd=@ THEN 
Xmin=2 
Ymin=@ 
Xmax=15 
Ymax=15 
Xstep=3 
Ystep=3 
ELSE 
BEEP 
INPUT "ENTER MINIMUM AND MAXIMUM X-VALUES",Xmin,Xmax 
BEEP 
INPUT "ENTER MINIMUM AND MAXIMUM Y-VALUES" ,Ymin, Ymax 
BEEP 
INPUT “ENTER STEP S072 POR R Ase) oer 
BEE 
INPUT “ENTER SJ ERPes te POR Ss eet eD 
END IF 
BEBE 
PRINT "IN:SP1;IP 2300,1800,8300,68Q@@: " 
PRINT "SC 0,100,@,1@00;TL 2,@;° 
Sfx=10Q0/( Xmax-Xmin) 
Sfyv=1@Q/(Ymax-Ymin) 
REE 
Icg=0 
INPUT “LIKE 10 BY-PASS CAGE (l=yY¥ ,O=N=NEFAUEI D7 eles 
IF Icg=1 THEN 762@ 
PRINT SFU OG Cera 
FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin)#*Sfx 
Pi Nimes Giau eres XA 
NEXT Xa 
PRINT "PA 10@,@;PU; ° 
PRING@ RS Para Gs Oeri. 
FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)«Sfy 
PR UNG fs epee) eae ale 
NEXT Ya 
PRINT =PA O00 ieee. 
FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin)#S fx 
PRINT| “RA ole ale: 
NEXT Xa 
PRINT "PA 10@,10@ PU PA 10@,@ POD" 
FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)#Sfy 
PRINT ©] PO Fa ug 27. yi 
NEXT Ya 
PRINT "PA 10@@,1@@ PU" 
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7430 
7435 
7440 
7445 
7450 
7455 
7460 
7465 
7470 
7475 
7480 
7485 
7490 
7495 
7500 
7505 
7510 
7515 
7520 
7525 
7530 
(oss 
7540 
7545 
7550 


foo5 


7560 
7565 
7570 
go (> 
758O 
(Ss 
7590 
G45 
7608 
7605 
7612 
7615 
762@ 
G25 
763@ 


7a 


7640 
7645 
7650 
fo55 
766@ 
wo55 
7672 
feo 


PRINT “PAO -2eern, 1.5,2° 
FOR Xa=Xmin TQ Xmax STEP Xstep 
X=(Xa-Xmin)#Sf» 
PRIN Teepe 0, 
[eea ee Vincee@mimeN PRIN? €P -1.5,-1;LB@;PR -1,0;LB";:Xa;:"" 
IF Xa=@ THEN PRINT “CP -.5,-1;L80" 
Xin=Q 
IF Xa MOD 1=@ THEN Xin=1 
IF Xa>=1@ THEN PRINT “CP -2,-1;LB"; Xa; "" 
ae aoe ea teen Kim=) FHEN PRINT “CP =-1.25,-1:LB";: Xa; °" 
IF Xao!l AND Xin=@ THEN PRINT "CP -2,-1;LB";: Xa; "" 
eee aal THENSRRINT "CP —1,-1;LB1.0" 


NEXT Xa 
SCN SNE felt Waa fa) = (fal 
Tht =2 PeMowlr y 


FOR Ya=Ymin TO Yma» STEP Ystep 
Y=(Ya-Ymind#Sfy 
PRINT “RA @aaeyY jee 
IF Iht=@ AND Yao@ THEN PRINT "PR 2,0" 
Reece ae yan eC THEN PRINT “€P —4\—. 25maeBO;>PR —-2,0;L6";Ya:"" 
Nee ae? THEN PRINT “CP -—2,-.25;LBQ" 
teva! AND lnt<Z2 THEN PRINT “@P -Sh=.25;3L8": Ya: 
MeeeyacQ AND Imt=2Z THEN PRINI “GP -492925;°L8"; Ya: "” 
teeraeOeaNe Yaciu AND Iht=Z THEN PRINT “CP -3,-.25:L8":Ya;"" 
[erase IneNerR ING €P —4,-.25;LB1.2° 


NEXT Ya 
TF Okd=@ THEN 
XlabelS="X" 
bo ee 
ELSE 
EEE 
PCW eNPeR A-LCABSEL” ,Xlabels 
BEEP 
INPUT “ENTER Y-LABEL" ,Yiabel$ 
END IF 


Aone ee Porn SO -10 CP", —-LEN( KlabelS)/2:°@:LB";XlabelS;~ 
Cr eee oer O92 -LEN(Ylabel®)/2*5/5; “D1 @,1;LB6°:Ylabel$; a 
PRONG eae O08 D1” 
Nerun=@ 
BEEr 
INPUT "WANT TQ PLOT DATA FROM A FILE (1=Y,Q=N)7" ,Ok 
Ok p=@ 
IF Ok=1 THEN 
BEER 
INPUT “ENTER THE NAME OF THE PLOT DATA FILE”,D_files 
ASSIGN @File TO D_ files 
IF Icomb<?@ THEN 7690 
Sx=0 
a 


Sy 
é 


hl 


WN 


"4 


‘ 
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768Q 
7685 
769@ 
fede 
77@Q@ 
77@5 
ius 
oral) 
7728 
yeti kas 
773@ 
(ss 
774@ 
7745 
775@ 
Woe 
7760 
idsks: 
7779 
tia 
778@ 
tts 
773@ 
7fbs\= 
7892 
7885 
7810 
Tea 
7820 
7825 
783@ 
FilciGa= 
784Q 
7845 
7850 
7856 
7860 
7865 
787Q 
(ots 
788 
72385 
789@ 
(ese 
79@@ 
79@5 
791@ 
fianks 
Tac 
(2es 


Sxy=@ 


Md=1 

BEEP 

INPUT “ENTER THE BEGINNING RUN NUMBER (DEF=1)" ,Md 
Noairs=9 

BEER 


INPUT "ENTER THE NUMBER OF X-~Y PAIRS STORED (DEF=9)",Npairs 
Bee 
INPUT “SELECT TUBE NUMBER (@=(0P  1=SECONDS ee es ltube 
Nrun=NruntNpairs 
PRINTER som 
slelele 
PRINT USING "4X,""Select a symbol!:""" 
PRINT USING "4x,"" Star 2 Plus sign 
PRINT USING “4x,"* Circle 4 Square""" 
PRINT USING "4xX,"" Rombus""" 
PRINT USING "4x,"" Right-side-up triangle 
PRINT USING "4x,"" Up-side-down triangle""" 
INPUT Sym 
PRINTER IS 7@5 
IF Sym=1 THEN PRINT "SM#" 
IF Sym=2 THEN PRINT "SM+" 
IF Sym=3 THEN PRINT "“SMo" 
Md=Itube*#9 
IF Md>1i THEN 
FOR IT=@ TO (Md-1) 
ENTER @File:ixX,yY 
NEXT I 
END IF 
FOR J=1i TO Npairs 
ENTER @File:Xa,Ya 
Sx=Sx+Xa 
Sy=SytVa 
Sx2=Sx2t+Xa%2 
Sxy=Sxy+Xat#vVa 
X=(Xa-Xmin)#Sfx 
Y=(Ya-Ymin)#Sfy 
IF Y>1@@ OR Y<@ THEN 791@ 
IF Sym73o THEN PRINT “She 
IF Syme4 THENSERIN] “Shed teen 
PRINT “PeMax.Y. > 
[Ee syme3 THENBPRINT {Shela oe 
TF Sym=4 THEN PRINT “UC2,4, 9950-35440 eo oe 
IF Sym=5 THEN PRINT *UC3,@,99,-5.-6,-3,6.4,692)-0>" 
IF Sym=6 THEN PRINT "UC9.5.5,9975,-2,-6,0,3,8:- 
IF Sym=7 THEN PRINT “UC, -S.6, 99 (See yee eee 
NEXT I 
BEEP 
INPUT "WANT TO LABEL (1=Y,@=N)?",I1bl 
IF Tibl=1 THEN 


J) Om NN 
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793@ 
2355 
7340 
7945 
795@ 
7359 
796@ 
7965 
797@ 
Gave 
798@ 
(ahete 
739@ 
Zass 
RAAe 
8@05 
R212 
8@15 
RAL 
8825 
RA50 
82355 
R242 
R@45 
BA5@ 
8O55 
B26 
B@6E5 
8272 
R75 
R282 
B@es 
8292 
Rags 
8122 
81@5 
811 
e115 
8122 
Bt2o 
81352 


a. age 


8140 
8145 
8150 
efoS 
B16@ 
alleys 
B17@ 
ei75 


IF Sym>3 THEN PRINT "SM" 
IF Sym<4 THEN PRINT "SR 1.4,2.4° 

PRINT “PA",Xal,Yal,"" 

MS Syme URISN) ep meee se 

IF Sym=4 THEN PRINT “UC2,4,99,0,-8,-4,0,0,8,4,@:" 
IF Sym=S THEN PRINT "UCZ,@,99,-3,-6,-3,6,3,5,3,-6:" 
IF Sym=6 THEN PRINT "UC@,5.3,99,3,-8,-6,0,3,8:" 

IF Sym=7 THEN PRINT “UCQ,-S.3,99,-3,8,6,@,-3,-8:3" 
PRINT "SM" 


IF Sym<4 THEN PRINT ee Zye 

Pit heen shel Oot 6s b8 sO files: “” 
Yal=Yal-5 

Been 


INPUT “WANT TO ADD ANOTHER STRING (1=Y,@0=N)?",Ias 
IF Ias=l1 THEN 
Beer 
INPUT “ENTER THE STRING" ,Label$ 
PeiVieeneerCnorn lael. 8,LB 3 Labels; ~” 
Soc) 7390 
END IF 
ena 1F 
BEEP 
INPUT “WANT TO COMBINE ANOTHER FILE? (1=Y,Q@=N)",Icomb 
ASSIGN @6File TO * 
IF Icomb<>@ THEN 7645 
Ils=l 
BEEP 
PNEUT = WANT TCO PLOT A VEAST-SQUARES LINE (1=DEF=YES,@=NQ)",Iis 
IF Ils=1 THEN 
2) Bale 
INPUT “SELECT CURVE TYPE (@=SOLID,1=DASHED)”,11t 
Lit=Lite+l 
PRINT SM" 
IF Jexp=@ THEN 
Bo=(Nrun* sxy-oy*Sx )7 CNrUn*SxZ2—5x "2 ) 
ERSE 
Bb=.75 
END IF 
Aa=(Sy-Bb#S~ )/Nrun 
Pie etl Sa 


CReN wUSINGw 10x, "a = ~*,MZ.35DE°":Aa 
ee wiUciNee hk Sy =") M7, SDE"; 8b 
PRENTER 15 705 

In=Q 


FOR Xa=Xmin TO Xmax STEP (Xmax~-Xmin) 
Ya=Aat+Xa*Bb 
Y=(Ya-Ymin )#Sfy 
X=(Xa-Xmin )*Sfx 
IF Y<@ THEN Y=@ 
IF Y>10@ THEN GOTO 822@ 
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8180 
8185 
8192 
ela 
8280 
8205 
821@ 
eZis 
8222 
e222 
8230 
e205 
8240 
8245 
825 
B2s5 
8262 
8278 
828@ 
B82 9@ 
8302 
R32@ 
8330 
8342 
8350 
B37 
838@ 
8390 
8400 
842@ 
8430 
844 
8450 
846@ 
847@ 
8480 
8490 


IF Tit=1 THEN 
PRINT SR Ar a so ee 
BSG 
In=Int+l 
Ir=in MOBrIIt 
TF fr=1SiHEN PRINT “PAs eee eee 
IF Ir=@ THEN PRINT "PA" ,X,Y, “PU” 
ENO IF 
NEXT Xa 
PRINT "PU" 
ENO IF 
Icomb=@ 
GOTO 762e 

END IF 

PRIN@ = eUSepeo. 

PRINT “PU PA @,@ SPQ" 

SUBEND 

DEF FNFmcal(1I,T,Pc) 

IF IT=@ THEN 
Mdot@=.@@0206+*+Pc-.@@414 
Mdot24=.@@225*Pc-.@@329 

END IF 

IF T=1 THEN 
Mdot@=.@02@4+*Pc-.@@615 
Mdot24=.@@216#Pc-.@@138 

ENO IF 

IF JT=2 THEN 
Mdot@=.@@2@8+*#Pc-.@Q36 
Mdot24=.@0@223*Pc-.@@3@01 

ENO IF 

IF IT=3 THEN 
Mdot@=.@@028*Pc-.@@122 
Mdot24=,0@@29#*Pc+.Q@@48 

END IF 

Mdt=(T#*(Mdot24-Mdot@)/24 )+Mdot@ 

RETURN Mdt 

FNEND 
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